Springer Tracts in Modern Physics 
174 



T.Schapers 

Superconductor/ 

Semiconductor 

Junctions 




Springer Tracts in Modern Physics 
Volume 174 



Managing Editor: G. Hohler, Karlsruhe 

Editors: J. Kiihn, Karlsruhe 

Th. Muller, Karlsruhe 
A. Ruckenstein, New Jersey 
F. Steiner, Ulm 
J. Triimper, Garching 
R Wolfle, Karlsruhe 

Honorary Editor: E. A. Niekisch, Jiilich 




Now also Available Online 

Starting with Volume 165, Springer Tracts in Modern Physics is part of the Springer LINK service. 
For all customers with standing orders for Springer Tracts in Modern Physics we offer the full 
text in electronic form via LINK free of charge. Please contact your librarian who can receive a 
password for free access to the full articles by registration at: 

http://link.springer.de/series/stmp/reg-form.htm 

If you do not have a standing order you can nevertheless browse through the table of contents of 
the volumes and the abstracts of each article at: 

http://link.springer.de/series/stmp/ 

There you will also find more information about the series. 



Springer 

Berlin 

Heidelberg 

New York 

Barcelona 

Hong Kong 

London 

Milan 

Paris 

Singapore 

Tokyo 



c 

Physics and Astronomy 



ONLINE LIBRARY 



http://www.springer.de/phys/ 




Springer Tracts in Modern Physics 



Springer Tracts in Modern Physics provides comprehensive and critical reviews of topics of current 
interest in physics. The following fields are emphasized: elementary particle physics, solid-state 
physics, complex systems, and fundamental astrophysics. 

Suitable reviews of other fields can also be accepted. The editors encourage prospective authors to 
correspond with them in advance of submitting an article. For reviews of topics belonging to the 
above mentioned fields, they should address the responsible editor, otherwise the managing editor. 
See also http://www.springer.de/phys/books/stmp.html 



Managing Editor 
Gerhard Hohler 

Institut fur Theoretische Teilchenphysik 

Universitat Karlsruhe 

Postfach 69 80 

76128 Karlsruhe, Germany 

Phone: +49 (7 21) 6 08 33 75 

Fax: +49 (7 21) 37 07 26 

Email: gerhardhoehler^physik.uni-karlsruhe.de 
http://www-ttp.physik.uni'karlsruhe.de/ 

Elementary Particle Physics, Editors 
Johann H. Kiihn 

Institut fur Theoretische Teilchenphysik 

Universitat Karlsruhe 

Postfach 69 80 

76128 Karlsruhe, Germany 

Phone: +49 (7 21) 6 08 33 72 

Fax: +49 (7 21) 37 07 26 

Email: johann.kuehn@physik.uni-karlsruhe.de 
http://www-ttp.physik.uni-karlsruhe.de/'^jk 

Thomas Muller 

Institut fiir Experimentelle Kernphysik 

Fakultat fur Physik 

Universitat Karlsruhe 

Postfach 69 80 

76128 Karlsruhe, Germany 

Phone: +49 (7 21) 6 08 35 24 

Fax: +49 (7 21) 6 07 26 21 

Email: thomas.muller@physik.uni-karlsruhe.de 

http://www-ekp.physik.uni-karlsruhe.de 

Fundamental Astrophysics, Editor 
Joachim Trtimper 

Max-Planck-Institut fur Extraterrestrische Physik 

Postfach 16 03 

85740 Garching, Germany 

Phone: +49 (89) 32 99 35 59 

Fax: +49 (89) 32 99 35 69 

Email: jtrumper@mpe-garching.mpg.de 

http://www.mpe-garching.mpg.de/index.html 



Solid-State Physics, Editors 

Andrei Ruckenstein 
Editor for The Americas 

Department of Physics and Astronomy 

Rutgers, The State University of New Jersey 

136 Frelinghuysen Road 

Piscataway, NJ 08854-8019, USA 

Phone: +1 (732) 445 43 29 

Fax: +1 (732) 445-43 43 

Email: andreir@physics.rutgers.edu 

http://www.physics.rutgers.edu/people/pips/ 

Ruckenstein.html 

Peter Wolfle 

Institut fiir Theorie der Kondensierten Materie 

Universitat Karlsruhe 

Postfach 69 80 

76128 Karlsruhe, Germany 

Phone: +49 (7 21) 6 08 35 90 

Fax: +49 (7 21) 69 81 50 

Email: woelfle@tkm.physik.uni-karlsruhe.de 

http://www-tkm.physik.uni-karlsruhe.de 

Complex Systems, Editor 
Frank Steiner 

Abteilung Theoretische Physik 

Universitat Ulm 

Albert-Einstein-AUee 11 

89069 Ulm, Germany 

Phone: +49 (7 31) 5 02 29 10 

Fax; +49 {7 31) 5 02 29 24 

Email: steiner@physik.uni-ulm.de 

http://www.physik.uni-ulm.de/theo/theophys.html 





Thomas Schapers 



Superconductor/ 

Semiconductor 

Junctions 

With 91 Figures 




Springer 




Dr. Thomas Schapers 

Forschungszentrum Jiilich 

Institut fur Schichten und Grenzflachen (ISG) 

52425 JUlich, GERMANY 

E-mail: th.schaepers@fz-juelich.de 



Library of Congress Cataloging-in-Publication Data applied for. 

Die Deutsche Bibliothek - ClP-Einheitsaafnahme 
Schapers, Thomas: 

Superconductor, semiconductor junctions/Thomas Schapers. - 
Berlin; Heidelberg; New York; Barcelona; Hong Kong; London; Milan; Paris; 
Singapore; Tokyo: Springer, 2001 
(Springer tracts in modern physics; Vol. 174) 

(Physics and astronomy online library) 

ISBN 3-540-42220-X 



Physics and Astronomy Classification Scheme (PACS): 74.50. +r, 74.80.Fp, 74.40.+IC, 
73.23.-b, 73-6i.Ey 



ISSN print edition: 0081-3869 
ISSN electronic edition; 1615-0430 

ISBN 3-540-42220-X Springer- Verlag Berlin Heidelberg New York 

This work is subject to copyright. All rights are reserved, whether the whole or part of the material is concerned, 
specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, reproduction on 
microfilm or in any other way, and storage in data banks. Duplication of this publication or parts thereof is 
permitted only under the provisions of the German Copyright Law of September 9, 1965, in its current version, and 
permission for use must dways be obtained from Springer- Verlag. Violations are liable for prosecution under the 
German Copyright Law. 

Springer- Verlag Berlin Heidelberg New York 
a member of BertelsmannSpringer Science+Business Media GmbH 

http://www.springer.de 

© Springer- Verlag Berlin Heidelberg 2001 
Printed in Germany 

The use of general descriptive names, registered names, trademarks, etc. in this publication does not imply, even in 
the absence of a specific statement, that such names are exempt from the relevant protective laws and regulations 
and therefore free for general use. 

Typesetting: Camera-ready copy from the author using a Springer KT^ macro package 
Cover design: design & production GmbH, Heidelberg 

Printed on acid-free paper SPIN: 10791132 57/3141/tr 543210 




Preface 



The present book covers the transport properties of superconductor/two- 
dimensional-electron-gas Josephson junctions. Starting with the basic el- 
ement, a superconductor/two-dimensional-electron-gas interface, phase-co- 
herent Andreev reflection in hybrid Josephson junctions is introduced and 
further on, multiterminal structures are discussed. Special care is taken on 
explaining the underlying theoretical concepts related to the transport mech- 
anisms, Employing a two-dimensional electron gas in a semiconductor as a 
normal conductor opens up the possibility to observe effects not found in 
purely metallic junctions. One example is the light sensitivity of the semi- 
conductor, which has a direct impact on the supercurrent in the Josephson 
junction. 

Many of the effects reported here rely on the fast technological progress in 
the epitaxial growth of III-V semiconductor heterostructures. By using these 
layer systems, fascinating quantum effects have been found. Two examples 
out of many are the quantized conductance in a point contact and electron 
optics using ballistic electron beams. By combining heterostructures with su- 
perconductors, many of the effects found in purely semiconductor systems can 
in a sense, be transferred to the superconducting state. A prominent example 
is the quantization of the critical current in a superconducting quantum point 
contact. 

The present book is based on the Habilitationsschrift that the author 
submitted to the Rheinisch-Westfalische Technische Hochschule (RWTH) 
Aachen. A large portion of the experimental results reported here are from 
measurements performed in Jiilich. By including results of other groups ac- 
tive in this research field, the author attempted to obtain a more or less com- 
plete coverage of the present status of the physics of superconductor/two- 
dimensional electron gas structures. Since most of the phenomena are de- 
scribed in a detailed fashion, this book should suit not only specialists active 
in the field but also interested students and scientists not familiar with these 
structures. 

The author wishes to express his gratitude to former and present col- 
leagues in Jiilich who were involved in the superconductor/semiconductor 
hybrid project; J. Appenzeller, A. I. Braginski, G. Crecelius, A. Forster, 
V.A. Guzenko, H. Hardtdegen, A. Kaluza, J. Kaufmann, T. Klocke, S.G. 
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Preface 



Lachenmann, J. Malindretos, K. Neurohr, D. Uhlisch and A.V. Ustinov. Valu- 
able discussions on the theory of superconductor/semiconductor structures 
with A, A. Golubov (Twente University, The Netherlands), P. Sarauelsson 
(Chalmers University of Technology) and M.Yu. Kupriyanov (Moscow State 
University) as well as the hospitality of the NTT Basic Research Laboratory 
during the final stage of the preparation of the manuscript are gratefully ac- 
knowledged. I am indebted to H. Liith for his support of the hybrid activities 
in Jiilich. Above all, I would like to thank my wife Gitta for her patience and 
understanding during the time when this manuscript was being written. 



Jiilich, May 2001 Thomas Schdpers 
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1. Introduction 



Introducing a semiconductor as a normal conducting layer in a weak-link 
structure [1] has induced a vast variety of new potential applications and 
unique physical effects [2, 3, 4, 5]. One motivation behind this approach is 
of course the aim of combining superconducting electronics, e.g. rapid single 
flux logic (RSFQ) [6], with state-of-the-art semiconductor circuits. The key 
element is the so-called Josephson field effect transistor (JoFET) [7], where 
the Josephson supercurrent can be controlled by a gate. JoFET operation 
has been demonstrated in the past by using several different semiconducting 
materials [8, 9, 10, 11]. However, since the voltage necessary to control the 
supercurrent is relatively large compared with the output voltage, voltage 
gain is difficult to obtain. A possible way towards better performance might 
be to use interference effects, i.e. Fabry-Perot resonances [12, 13, 14]. 

As a promising alternative to a JoFET, nonequilibrium carriers injected 
into the normal-conducting area can also be used to control the supercurrent 
in a weak-link structure [15, 16, 17, 18, 19, 20, 21, 22, 23]. Here, the underlying 
mechanism is that the occupation of the Andreev levels, which determine 
the supercurrent flow, is altered by the injected carriers. It was shown that 
relatively small injection currents are sufficient for an efficient supercurrent 
control [22]. 

The application of superconducting three-terminal structures such as 
JoFETs or injection-type devices in real electronic circuits is still an open 
question. Many technological obstacles have to be overcome before this is 
done. Nevertheless, superconductor/two-dimensional-electron-gas structures 
can serve as an almost ideal system when fundamental questions about su- 
percurrent flow are to be addressed. Novel quantum effects can be observed, 
which rely on the large Fermi wavelength found in semiconductors compared 
with purely metallic structures. A prominent example is the superconduct- 
ing quantum point contact, where a quantized supercurrent was predicted 
[24, 25] and experimentally observed [26]. 

Superconductor/two-dimensional-electron-gas structures can be classified 
into two groups. First, an inversion layer, e.g. on the surface of InAs, can 
be employed as a two-dimensional electron gas (2DEG). In this case the 
transport through the structure usually has to be described in the so-called 
dirty limit, where the coherence length ^ is larger than the elastic mean 
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free path. The carrier transport is diffusive in this case. In the alternative 
regime, the clean limit, the coherence length is smaller than the elastic mean 
free path l^i. In order to have access to the latter regime, two-dimensional 
electron gases in semiconductor heterostructures are used. Here, the spatial 
separation of the donors from the electrons in the quantum well leads to an 
efficient suppression of the mobility-reducing ionized-impurity scattering. As 
a result, an elastic mean free path of several microns can be obtained easily at 
low temperatures. Since the relevant dimensions in a superconductor/2DEG 
structure are mostly smaller than lei, the carrier transport can be considered 
as ballistic. In this regime, interesting quantum effects have been predicted 
theoretically [24, 25] and confirmed experimentally [26, 27, 28]. 

Concerning the choice of an appropriate semiconductor, the height and 
thickness of the Schottky barrier at the interface is one of the most important 
parameters. In order to obtain a good coupling between the superconductor 
and the semiconductor, the transparency of the interface barrier should be 
as large as possible. In this respect InAs is an ideal candidate, since the 
Fermi level is pinned above the conduction band edge. In addition, the large 
electron concentration of the 2DEG at the surface and the low effective elec- 
tron mass lead to a large coherence length in the normal-conducting material 
[ 8 , 29, 30, 31]. A drawback is that the conductance of the 2DEG at the sur- 
face can result in insufficient isolation between adjacent junctions. However, 
if the InAs surface is covered with an additional InAlAs layer, this surface 
conductance can be suppressed, as shown by Lachenmann et al. [32]. In ad- 
dition to InAs, Si [9] and InGaAs [33] have been used as a semiconducting 
material in hybrid weak-link structures. 

In order to contact a superconductor to a 2DEG embedded in a semi- 
conductor heterostructure, two different approaches can be pursued. If, ul- 
timately, a high mobility is needed, AlGaAs/GaAs heterostructures are the 
material of choice [34, 35]. However, a drawback of this system is that owing 
to the diffusion process necessary for the formation of the superconductive Sn 
contact, the junction geometry is not precisely defined. This problem can be 
circumvented if the heterostructure is etched down to or below the 2DEG so 
that, subsequently, the superconductor can be contacted from the top or from 
the side to the channel layer [14, 36, 37, 38]. The latter method is only ap- 
plicable if the Schottky barrier is sufficiently low, e.g. in InAs-based channel 
layers. 

For many of the effects discussed here, results on structures based on InAs 
surface layers as well as on Ino. 57 Gao. 43 As/Ino. 77 Gao. 23 As/InP heterostruc- 
tures have been reported [39, 40]. For the latter, the 2DEG is located in 
a strained Ino. 77 Gao. 23 As channel layer. Owing to the high mobility, even 
ballistic transport through a quantum point contact could be demonstrated 
[41, 42, 43]. Because of the high indium content, the Fermi level pinning oc- 
curs approximately at the conduction band edge. This property leads to a 
low-resistance coupling to the superconductor [44, 45]. 
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This book is organized as follows. The basic principles of carrier transport 
through a single superconductor/normal-conductor interface are discussed in 
Chap. 2. The properties of two-dimensional electron gases and the prepara- 
tion steps necessary to fabricate a junction are described in Chap. 3. The fol- 
lowing chapter deals with the characterization of single superconductor/two- 
dimensional-electron-gas interfaces. Special care was taken here in optimiz- 
ing the superconductor/normal-conductor interface properties by employing 
several alternative methods to clean the semiconductor surface before the 
deposition of the superconducting material. In Chap. 5 the voltage-carrying 
state of superconductor/two-dimensional-electron-gas junctions structures is 
investigated, whereas the next chapter is dedicated to the superconductive 
state. As outlined above, gate control of the carrier concentration is one of the 
unique features of a semiconducting channel. In Chap. 7, first, the Josephson 
field effect transistor is introduced, while in the second part more funda- 
mental effects, e.g. quantization or fluctuations of the critical current, are 
described. The effect of nonequilibrium carriers on the Josephson supercur- 
rent is presented in Chap. 8. 



2. Carrier Transport Through a 
Superconductor /Normal-Conductor Interface 



In this chapter the carrier transport through the simplest structure, a super- 
conductor/semiconductor bilayer is, discussed. A detailed understanding of 
this key structure is essential for the study of more sophisticated structures, 
i.e. Josephson junctions or three-terminal structures, as presented in the fol- 
lowing chapters. First, we shall introduce Andreev reflection as a fundamen- 
tal mechanism related to transport through a superconductor/semiconductor 
interface. Next, the Blonder-Tinkham-Klapwijk (BTK) model [16] is ex- 
plained, which treats the carrier transport through a non-ideal interface with 
a barrier in between. Because of its simplicity, the BTK model is referred 
to frequently to describe transport through a superconductor/semiconductor 
interface. However, the experimental results obtained on single interfaces sug- 
gest that some extensions, e.g. the inclusion of the proximity effect, are nec- 
essary in order to obtain a better theoretical description of the experimental 
findings. 

For the discussion of the basic transport phenomena it is often sufficient 
to use an ideal normal conductor (N), which means that this normal conduc- 
tor consists of the same material as the superconductor (S). This simplifies 
the theoretical description considerably, since all material parameters, e.g. 
electron concentration, are identical. However, as we shall see later, the cou- 
pling of a semiconductor to a superconductor can differ significantly from the 
above-mentioned idealized case, i.e. the effective electron mass or the elec- 
tron concentration can be much smaller than the corresponding values for 
the superconductor. Therefore, in some cases it is necessary to extend the 
description of an S/N interface in order to take the specific properties of the 
superconductor/semiconductor interface into account. 



2.1 The Bogoliubov Equation 

In 1957, Bardeen, Cooper and Schrieffer (BCS) [47] succeeded in explaining 
the appearance of the superconductive state at an atomistic level. They found 
that two electrons with opposite wave vector and spin form a Cooper pair 
(fc t) i)- electron pairing is mediated by a weak electron-phonon cou- 
pling, the Frohlich interaction. Owing to this attractive interaction. Cooper 
pairs, which are bosonic particles, condense to an energetically lower ground 
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state. Despite the lower ground state energy, the kinetic energy of the elec- 
trons is increased compared with a noninteracting Fermi gas. However, this is 
overcompensated by the pairing energy, which is gained owing to the attrac- 
tive electron-phonon interaction. Because of the increased kinetic energy, the 
probability Vq to find a Cooper pair (fe — fc J.) at T = 0 is not step-like as 
for a noninteracting Fermi gas, but rather is smeared out around the Fermi 
wave vector kp, as shown in Fig. 2.1. 




Fig. 2.1. Probability Vq that a pair state (k f, —k f ) is occupied along a particular 
k axis. Two single-particle exitations are shown, with electrons placed at fci and k2- 
Since these electrons do not form a Cooper pair, the states at —ki and —k2 must 
be empty 



The vanishing resistance of a superconductor is closely related to the 
opening of an energy gap Aq between the ground state of the Cooper pair 
condensate and the spectrum of excited single-electron states. Owing to this 
energy gap, scattering processes, which are responsible for the resistance of a 
conductor, are suppressed. In order for an electron to occupy a single-particle 
state, the electron energy must be larger than the gap Aq. 

Let us have a closer look at what happens if a state is occupied by a 
single electron. In Fig. 2.1 this situation is depicted for an electron with a 
wave number k\ below fcp and for a wave number k 2 above fcp- In this case, the 
corresponding states are occupied with probability one. Since these electrons 
must not form a Cooper pair, the corresponding states at the opposite side 
of k space —ki and — ^2 must be empty. As can be seen in Fig. 2.1, at k\ the 
system gains a small portion of additional electron character, which can be 
quantified by 1 — Vq. Owing to the empty state at — fci, the additional hole 
character, expressed by Ug, is larger. Thus, as a net result, an excitation with 
a wave vector below k-p can be regarded as hole-like. The opposite situation 
is found for a wave vector located above kp. As shown in Fig. 2.1, the system 
gains more electron character at /c 2 than hole character at — ^ 2 , thus this 
excitation can be regarded as an electron- like quasiparticle. 

The quasiparticle states introduced above can be described by the Bo- 
goliubov equation [18], 
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f H{r) A{r) \ f Uk{r) \ Uk{r) \ 

A*{r) -H{r) ) Vk{r) ) ~ \ Vk{r) ) ' 



( 2 . 1 ) 



The solutions of the Bogoliubov equation are electron-like or hole-like quasi- 
particles represented by the vector (ufc,Ufe). Here, H{r) is the one-electron 
Hamiltonian, defined as 

= + ( 2 . 2 ) 

where ^ is the electrochemical potential, m* is the effective electron mass 
and U{r) is a scalar potential. The coupling between the two components of 
the vector (uk,Vk) is provided by the superconducting pair potential A{r). 
In a metallic superconductor m* is identical to the free electron mass rUe- 
Following the discussion above, the component Vk represents the probability 
amplitude for a hole-like state, while Uk represents the probability of an 
electron- like state. The two are coupled by the pair potential A{r). The 
character of the quasiparticle, electron-like {\uk\^> |wfcp) or hole-like (|Mfcp < 
|ufep), is determined by the predominant component of the vector (uk,Vk)- 
In many cases, the common spatial component g(r) can be split off from 
(Uk,Vk)'- 



e.g., for a homogeneous superconductor with A(r) = Aq, plane wave solutions 
exp(ifc • r) can be assumed for g{r). In this case uo and vq are given by 



o to 

II 




E J ’ 




(2.4a) 


II 

(N O 


1 — Mq , 






(2.4b) 


where the energy eigenvalues are 












1/2 


(2.5) 


E = ± 


V 2me 


-g] +A*qAq 





The quantity Vq in (2.4b) is identical to the above-discussed probability of 
finding a Cooper pair in the BCS theory. Thus, for the case where |fc| is larger 
than the Fermi wave vector, the quasiparticle is electron-like, since > 
|uoP, while for |fe| < kp the character of the particle is hole-like (|uoP < IrioP)- 
As can be seen in (2.5), a gap of Aq exists in the excitation spectrum. Within 
this gap no single-particle states exist in the superconductor. 

For a normal conductor, the superconducting pair potential is zero. This 
implies that the two components of (2.1) are decoupled and pure electrons 
(1,0) and holes (0,1) exist. For a hole state (0,1), the energy eigenvalues are 
positive for |fe| < kp and negative for |fc| > kp. This can easily be seen if the 
Bogoliubov equation is solved for a hole state with Aq = 0. As an important 
consequence, the direction of the group velocity 
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Vk = iVfcE' (2.6) 

n 

is, for a hole, opposite to its wave vector k. This is in contrast to the case of 
an electron, where Vk and k are oriented in the same direction. 

In addition to its application to homogeneous systems, as discussed above, 
the Bogoliubov equation can also be used as a very powerful tool to describe 
inhomogeneous systems like superconductor/normal conductor interfaces. 



2.2 Quasiparticle Transport Through 
a Superconductor/Normal-Conductor Interface 

Here, the models describing carrier transport through a superconductor/ 
normal-conductor interface are introduced. The Andreev reflection process 
is the basic ingredient of the theoretical description. However, in order to be 
able to model nonideal interfaces, normal scattering and the proximity effect 
have to be included also. 

2.2.1 Andreev Reflection 

If a normal conductor is coupled to a superconductor a unique reflection pro- 
cess, namely Andreev reflection, can be observed [49]. This reflection mech- 
anism is depicted in Fig. 2.2. Let us assume an electron with an energy 
E < Aq slightly above the Fermi level fj, moving towards the normal con- 
ductor/superconductor interface. Since no quasiparticle states are provided 
in the superconductor, transmission is excluded. Furthermore, no normal re- 
flection is allowed, because there is no barrier present at the interface, which 
can absorb the momentum difference. However, a Cooper pair can be formed 
in the superconductor, with the consequence that an additional electron is 
taken from the completely filled Fermi sea of the normal conductor. In a 
degenerate electron gas this electron vacancy can be interpreted as a hole 
carrying a positive charge. 

For the generation of a Cooper pair in the superconductor, the wave 
vector of the electron taken from the Fermi sea needs to have a direction 
opposite to the wave vector of the incident electron. Since in a completely 
filled Fermi sea the total wave vector is zero, the remaining wave vector of the 
system after one electron has been removed for the Cooper pair has the same 
direction as the incident electron. The wave vector of the system is identical 
to the wave vector of the hole. Thus, the incident electron and the created 
hole have a wave vector in the same direction. However, as mentioned above, 
the wave vector and the group velocity of a hole are in opposite directions. 
Consequently, the hole takes the same path as the incident electron in the 
reverse direction; therefore this process is called retroreflection. Owing to 
the opposite charge and opposite group velocity of electrons and holes, this 
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also implies that the conductance is twice as large as for an ideal normal 
transmission through the interface. 



electron 

E j O ► 






,'^0 


O’- — 




. ■ ■ M' 


hole 




► X 




Fig. 2.2. (a) Energy diagram of the Andreev reflection process: an incident electron 
from the normal-conductor side is retroreflected as a hole, while a Cooper pair is 
formed in the superconductor, (b) In real space: in contrast to the normal reflection 
process, the retroreflected hole takes the same path as the incident electron in 
reverse 



For an incident electron with \E\ > Aq, quasiparticles, electron- or hole- 
like, are excited in the superconductor. In this case not only Andreev reflec- 
tion but also normal specular electron reflection take place with a certain 
probability. An Andreev reflection probability less than one is found even for 
I if I < Aq if a potential barrier is introduced at the superconductor /normal- 
conductor interface, because now normal reflection at the barrier is allowed. 

2.2.2 Blonder Tinkham Klapwijk Model 

A model which describes Andreev and normal reflection at a superconduc- 
tor/ normal-conductor (SN) interface, as well as the transmission of particles 
through this interface, was first developed by Blonder, Tinkham and Klapwijk 
(BTK) [16]. A rigorous microscopic derivation of this model was presented a 
few years later by Zaitsev [50] and Arnold [51]. In the BTK model a poten- 
tial barrier at the SN interface, e.g. an oxide layer or a Schottky barrier, is 
approximated by a ^-shaped barrier as depicted in Fig. 2.3a. In the original 
BTK model no difference of the Fermi energy between the superconductor 
and the normal conductor was assumed. However, in the case of a semicon- 
ductor in contact with a superconductor, the Fermi energies can differ sig- 
nificantly, since the electron concentration in the semiconducting material is 
usually much lower than in the superconductor. Therefore, if these materials 
are brought together, the alignment of the Fermi energy leads to a poten- 
tial step of height Uq at the bottom of the conduction band at the interface 
(Fig. 2.3a). Similarly to the case of a ^-shaped barrier, this also leads to en- 
hanced normal reflection at the normal-conductor/superconductor interface. 
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b) 




Fig. 2.3. (a) Model of a normal-conductor/superconductor interface. A step-like 
increase of the pair potential A{x) is assumed at the NS boundary. The potential 
step between the superconductor and the normal conductor is f/o- In addition, a 5- 
shaped barrier is assumed at the interface, (b) Andreev reflection process where an 
electron (e) enters the normal-conductor/superconductor interface. The orientation 
of the arrows represent the direction of the group velocity 



For simplicity, a step-like increase of the superconducting pair potential 
is assumed at the interface: A{x) = 0{x)Aq. The potential U{x) contains 
the step Uq as well as the J-shaped barrier at the interface: 

U{x) = Uo9{-x) + ^-^ZS{x) , (2.7) 

me 

where fcps = 2mefJ,/h^ is the Fermi wave number in the superconductor. 
The height H of the ^-shaped barrier is expressed by the dimensionless Z 
parameter 

Z = Hme/h^kps ■ (2.8) 

The scattering process of an incident electron wave 




from the normal-conductor side N is depicted in Fig. 2.3b. As a result, 
electron-like (uq,vo) and hole-like {vq,uo) quasiparticles are transmitted into 
the superconductor: 

V^t,ans(x) = C+ ( 2 ) ( Z ) • (2.10) 

At the same time a fraction of the incoming wave is reflected as an electron 
(1, 0) and a hole (0, 1): 

= &+ ( J ) + «+(?) ■ (2-11) 

The corresponding wave numbers in the normal conductor and superconduc- 
tor are 
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\/ ^FN + (2m*/h^)E , 


(2.12a) 


kh = 


\Z^FN “ {2m*/n^)E , 


(2.12b) 


k(, — 


^Jkis + {2m,/n^){E^-Aiy/\ 


(2.13a) 


kh = 


J fcfs - {2m^l?y){E^ - Z\g)i/2 , 


(2.13b) 



respectively; fcpN = y (2m*/?i^)(/r — f/o) is the Fermi wave number in the 
normal conductor. All particles involved in the scattering process are depicted 
in the E(k) dispersion relation, shown in Fig. 2.4. From the continuity of the 




Fig. 2.4. Schematic illustration of the E{k) dispersion relation in the normal con- 
ductor (N) and in the superconductor (S). The filled circles denote electrons and 
electron-like particles, the open circles holes and hole-like particles. The arrows in- 
dicate the direction of the group velocity. The axis of the hole energy Ah points 
downwards 



wave function at the normal-conductor/superconductor interface (x = 0), 

yin(O) -h yrefl(O) = V'trans(O) , (2-14) 

and the boundary condition for the derivative, 

^(n(O-) + V';efl(0”) - V'Lns(0+) = 2fcFS^V’(0) , (2.15) 



the coefficients a+ ,..., d+ can be determined: 



a+ 



b+ 



UqVo 

•> 

7 

jup - + q) 

7 

Up [(1 -I- r)/2 - iZJ 



(2.16a) 

(2.16b) 



c+ 



7 



(2.16c) 
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^ 7 

Here, we define 7 = Uq(j> + 1) — VqP, q = jr + (1 
+ (r — l)^/4r. The ratio 



(2.16d) 



r^)/4r and p = 



fPN 

WFS 



(2.17) 



is the Fermi velocity mismatch between the normal-conductor and the su- 
perconductor. In case of |if| < Z\o the quasiparticle waves are exponentially 
damped since fee and fch contain an imaginary component. 

For a hole injected from the normal-conductor side, the reflection coef- 
ficients a_ ,..., d- can be calculated from the coefficients determined above 
by simply replacing Z by — Z, i.e. a_{Z) = a+{—Z). In case of injection of 
an electron-like particle from the superconductor into the normal-conductor, 
time reversal symmetry leads to the following transmission coefficients: 



4 = c±{ul-vl)/- 

d'± = d±{ul-vl)/ 



(2.18) 



The set of coefficients a±,..., d±, c(j_ and d(j_ will be used explicitly for the cal- 
culation of the supercurrent in a superconductor/two-dimensional-electron- 
gas/superconductor junction in Sect. 6.1.2. 

The electrical current I{V) through a single SN interface can be calculated 
from the probability currents [46]. In order to find an expression for I{V) it is 
sufficient to consider the contributions on the normal-conductor side, i.e. the 
Andreev-reflected contribution A{E) and the normal-reflected contribution 
B{E)-. 

I{V) = r \h{E+eV) - /o(A)][l + A{E) - B{E)] dE .(2.19) 

^ ^ J —00 

Here, fo{E) is the equilibrium Fermi distribution function, V is the voltage 
drop at the interface and W is the contact width. The Andreev reflection 
coefficient is defined as A{E) = a^a+, which is the probability current Jp 
of the Andreev reflected quasiparticles divided by ripN [16]. Similarly, the 
normal reflection coefficient B{E) is given by 4^+- 

A unique feature of the Andreev reffection is that the total current and 
thus the conductance are enhanced by this process. This is due to fact that 
the Andreev-reflected quasiparticle moves with the opposite group velocity 
to the incident quasiparticle but possesses the opposite charge at the same 
time. The contribution of the Andreev reflection can result in a so-called 
excess current /exc, which is the current added to the normal current /n = 
y/i?N at eV ^ Aq [16]. i?N is the normal-state resistance in the absence of 
superconductivity. 

For an isolated SN junction the Fermi velocity mismatch r can be exactly 
reproduced by shifting the Z parameter to a higher effective value, Zes = 
y/p =[Z‘^ -I- (1 — r)^/4r]^/^ [52]. If we define 7 = vq/uq, the Andreev and 
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normal reflection coefficients A{E) and B{E) used in (2.19) can be written 
as [53] 



A{E) 

B{E) 



|1 + Z2^(1-77^)P ’ 
\1 + {1-V^)Z!^\^ 



(2.20a) 

(2.20b) 



In Fig. 2.5, the coefficients A{E) and B{E) are plotted for Z^ff = 0 and 
^eff = 0.5. 





Fig. 2.5. (a) Andreev reflection coefficient A as a function of EjAo for ZeS = 0; (b) 
Andreev reflection coefficient A and normal reflection coefficient B for =0.5 



In the expression for the current through the interface (2.19), only nor- 
mal incidence of the quasiparticles was considered. The dependence of the 
Andreev reflection coefficient on the angle of incidence 9 was calculated by 
Hoonsawat and Tang ['' ']. Mortensen et al. [55] could show in an extended 
model, which also includes the Fermi velocity mismatch, that the Andreev re- 
flection process is forbidden when 9 is larger than a critical angle. Extensions 
of the BTK model concerning a finite-thickness barrier or a superconductive 
interlayer have been developed by Kupka [56] and by Sanchez-Cahizares and 
Sols [57]. 
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2.2.3 Inclusion of the Proximity Effect 



One key assumption of the BTK model is the step-like superconductive 
pair potential A{x) at the superconductor/normal-conductor interface, where 
Z\(x) is usually set to the value of the bulk superconductor. However, in many 
cases it can be necessary to assume a modified pair potential at the boundary 
with the normal-conducting layer in order to obtain a consistent interpreta- 
tion of the measured results [45, 58, 59]; i.e., if a superconducting material is 
deposited on a normal-conductor, the properties of the first atomic layers can 
be degraded owing to residual atoms or molecules on the surface. Hence, for a 
more realistic description of the superconductor/normal-conductor system, it 
can be worthwhile to assume an additional normal-conductive interlayer (N) 
in between. The normal-conducting electrode can, in our case be specified to 
be a semiconductor (Sm). Owing to the additional N layer, the semiconductor 
is in this case not in direct contact with the superconductor anymore. The 
complete layer sequence would thus be SN/Sm. However, if the N layer is suf- 
ficiently thin, the superconducting properties will be preserved to a certain 
degree. The underlying mechanism is the proximity effect [60]. Here, Cooper 
pairs from the superconductor can leak into the normal-conductor and diffuse 
a certain distance before they decay. This decay takes place on a length scale 
given by the coherence length ^n- 

For the discussion of the basic principle of the proximity effect, an SN 
sandwich in the dirty limit will be considered first. The dirty limit means 
that the coherence lengths £md are much larger than the corresponding 

elastic mean free paths and Zei,N in the superconductor and normal- 

conductor, respectively. The coherence lengths ^s.n in the dirty limit are 
given by 



Cs.N = \jKDs^^/2'Kk'B,Tc , (2-21) 

where I?s,n is the corresponding diffusion constant and Tc is the critical tem- 
perature of the superconductor. As mentioned above, the N layer will usually 
be composed of a degraded superconductor material. The metallic supercon- 
ducting films used for superconductor/semiconductor structures contain a 
large number of scattering centers. Therefore, the dirty-limit case is justified 
for the S as well as for the N layer. Independently, the semiconductor at- 
tached to this SN sandwich can be either in the clean or in the dirty limit. In 
the former case the coherence length is much smaller than the elastic mean 
free path. 

The large number of scattering centers in the SN sandwich makes it virtu- 
ally impossible to apply the Bogoliubov equation directly, in order to describe 
the system. A theoretical tool that allows one to handle the dirty-limit case 
in the S and N layers is provided by the Usadel equation [61], which follows 
from the Eilenberger equations [62] under momentum averaging: 



Cn,s^n,s(^) + i , rp 
tt/cbJc 



sm PN.sfxj 



tt/cbJc 
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Here, the system is described by a formally introduced complex angle 
(^n,s{E,x), where cos9^^s{E,x) is the normal angle-averaged Green’s func- 
tion G{E,x) and sin 0 ns (if, x) represents the anomalous angle- aver aged 
Green’s function F{E,x). The pair potential As(x) in the superconductor 
has to be determined from the self-consistency equation 



^s(a^) T r ^s(a^) 

tt/cbTc Tc Tc “ -KUlnksTc 

o;n>0 ^ 



sin 6^{E = icj„, x) 



= 0 . 



(2.23) 



Here, is the normalized Matsubara frequency given by (2n -|- 1)T/Tc. In 
the normal-conductor the pair potential Z\n( 3;) is assumed to be zero. The 
Green’s functions G{E,x) and E{E,x) contain a direct physical meaning. 
The real part of G{E, x) can be interpreted as the quasiparticle density of 
states (DOS) 

D0S(T;, x) = Re[G{E, x)] , (2.24) 

while Re[F(£l,x)] may be interpreted as the DOS of correlated particles. 

The proximity effect in an SN sandwich was treated theoretically by 
Kupriyanov and Lukichev [63]. They succeeded in deriving boundary con- 
ditions for the normal and anomalous Green’s functions G{E, x) and E{E, x) 
at the SN interface [64, 65]. Since the Usadel equations are not valid at 
distances of the order of /ei,s £md /ei,N from the interface, the boundary con- 
ditions, as given below, had to be derived from the more general Eilenberger 
equations [62]: 



7bCn0n(O’^) 

7eN0N(O+) 

0s(-oo) 

^N(dN) 


= sin [0 n(O+) - 0s (0 )] , 
= ^S^s(O-), 

= arctan[iZ\o(T)/£’] , 

= 0 . 


(2.25) 

(2.26) 

(2.27) 

(2.28) 


The parameter 






2 ^el.N 

3^' 




(2.29) 


in the first boundary condition (2.25) is a measure of the barrier strength at 
the interface. This is plausible, since the normal-transmission coefficient D at 
the SN interface in (2.29) is connected to the Z parameter in the BTK model 
by D~^ = 1 + Z'^ . In the expression for 7 b, < ... > denotes angle averaging. 
The pair-breaking parameter 


_ Qsis 




(2.30) 



in the second boundary condition (2.26), is a measure of the suppression of the 
order parameter in the superconductor due to the proximity effect; gN and £»s 
are the specific normal-state resistances of the normal and superconducting 
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material, respectively. For a thiir N layer of thickiress cZn with c^n/Cn 1, 7 
can be replaced by 



7m = 7 



Cn 



(2.31) 



The remaiiriirg two boundary coirditions (2.27) and (2.28) match the solution 
at —00 to a bulk superconductor at —00 and at the free interface of the N 
layer at cZn, respectively. 

The effect of the proximity effect oir the Cooper pair density F(x) in an 
SN sandwich is depicted schematically in Fig. 2.6. As cair be seeir here, F{x) 
decreases iir the viciirity of the interface on the superconductor side. Oir the 
normal-conductor side, F{x) is not zero but rather decreases monotonically. 
As shown in Fig. 2.6, 7 b determines the step iir the Cooper pair density F{x), 
while 7 is a measure of the suppression of F{x) in the superconductor near 
the interface. 



Foo 
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Fig. 2.6. Schematic illustration of the spatial dependence of the Cooper pair den- 
sity F{x) for a representative, fixed energy E — Eq ai an SN boundary 



As cair be seeir in Fig. 2.7, the density of states in an SN sandwich, 
calculated according to (2.24), differs significantly from the DOS of a bulk 
superconductor. At the SN interface, the maximum of the density of states 
around Z\s is widened and lowered if the coupling between the layers is en- 
hanced by decreasing 73. At the free surface of the N layer {x = -l-dN) a 
gap An is developed in the density of states owing to the coupling to the 
superconductor.^ This energy gap approaches the gap of the superconductor 
Aq if 7 b is decreased. No gap is expected if the N layer is very thick. This 
situation was considered by Belzig et al. [66]. Near the SN interface the DOS 
in the N layer has an induced maximum below the gap of the superconductor. 
This maximum is shifted towards zero if the distance from the interface is 

^ The gap An should not confused with the spatially dependent pair potential 
An(x). 
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increased. For very long distances the DOS approaches the Fermi-level DOS 
in the normal state. 




E/A, 



Fig. 2.7. Density of states in a thin N layer (dN <S ^n) of an SN sandwich nor- 
malized to the normal-state value at the free surface {solid lines) and at the SN 
boundary {dashed lines) for 7 m = 7dN/^N = 0.1 and for various values of 7 b. Data 
taken from [58] (© (1996) by the American Physical Society) 



In an SN sandwich with a thin layer and a small ratio 7m/7Bi the spatial 
gradients in S can be neglected, and /\s(a^) = ^o- Then the temperature 
dependence of the energy gap in the N layer can be determined from [58] 



ziN(r) 



^o(T) 

l + 7Bx/^^(T)-Zi^(r)/7rT, ■ 



(2.32) 



The two gaps zIn and Z\q merge if the temperature is increased towards Tc. 

Golubov and Kupriyanov [53, 67] succeeded in generalizing the BTK 
model by including the proximity effect. Instead of an SN bilayer, they as- 
sumed an SN/c/N' structure. Here c is a constriction, which ensures that 
the proximity effect in the SN sandwich does not extend to the N' layer. 
For the structures discussed later, N' is identical to the semiconductor layer 
(Sm). The constriction c is naturally given by the reduced dimensions of the 
conducting channel in the semiconductor. In contrast to the BTK model, the 
semiconductor is not in direct contact with the superconductor via the con- 
striction but rather with a normal-conducting layer (N), which is modified by 
the proximity effect. Owing to the proximity effect, the density of states at 
the interface with the semiconductor is altered, whereas in the BTK model, 
the quasiparticle density of states of the superconductor at the interface cor- 
responds to its bulk value DOS{E) = E/ — A^. 
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Owing to the modified density of states in the N layer, the Andreev reflec- 
tion coefficient A{E) and the normal-reflection coefficient B{E) of the BTK 
model, as given by (2.20b), have to be replaced by 



ME) = 


\E{E,x = 0-)\^ 


(2.33) 


|l-h2Z2-hG(A,x = 0-)|2 ’ 


B{E) = 


4^2(1 -hZ)2 


(2.34) 


|l-h2Z2-hG(A,x = 0-)|2 ■ 



In these expressions the Green’s functions G(A, x) and F{E^ x) have to be 
taken at the N/c boundary in the N region. The Z factor is determined by 
the barrier in the constriction between the N and the Sm layers. In Fig. 2.8, 
A{E), B{E) and the transmission probability T{E) are plotted as a function 
of energy. In accordance with the lower gap at the free interface of the N 




E/A, E/A, 

Fig. 2.8. Coefficients for Andreev reflection A{E), normal B{E) reflection and 
transmission T{E) for Z — 1 and 7 m = 0.1. In (a) ye = 0 and in (b) yn = 1 was 
assumed. Data taken from [58] (© (1996) by the American Physical Society) 



layer due to the proximity effect, the peak position in the Andreev reflection 
coefficient as a function of energy is found at E = Z\n. For large values of ye, 
a second peak near As is found in the Andreev reflection coefficient, e.g. for 
7b = 1 (Fig. 2.8b). This originates from the second maximum in the density 
of states (see Fig. 2.7). 

So far, it has been assumed that in the SN junctions the superconductor 
is decoupled sufficiently from the normal-conductor that no proximity effect 
has to be considered between these layers. The proximity layer appeared 
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exclusively on the side of the superconductor, where an additional nonsuper- 
conducting interlayer was assumed. However, if the coupling between the 
superconducting and normal-conducting electrodes is sufficiently strong, the 
proximity effect between these layers has to be considered. As a result the 
interface pair potential in the superconductor can be reduced compared with 
its bulk value Aq [68] . The Andreev reflection probability for an SN junction 
including the proximity effect was calculated by van Son et al. [69]. It was 
shown that, owing to the decrease of the pair potential in the superconductor 
at the interface, the Andreev reflection probability is decreased for E > Aq 
compared with the BTK result. 



3. Fabrication of Superconductor/ 
Two-Dimensional-Electron-Gas Structures 



Two-dimensional electron gases in semiconductors have some interesting 
properties, which differ significantly from those of the bulk material, e.g. the 
electron mobility of a two-dimensional electron gas can be larger by orders of 
magnitude at low temperatures compared with a bulk semiconductor. Here, a 
brief overview of different realizations of two-dimensional electron gases and 
their physical properties is given. After that, the principal fabrication steps 
required to fabricate a superconductor/semiconductor junction are outlined 
for a junction based on a surface inversion layer, as well as for a junction 
with a buried two-dimensional electron gas in a heterostructure. 



3.1 Two-Dimensional Electron Gases 

In a two-dimensional electron gas the electrons can move freely within a 
plane, while in the direction perpendicular to this plane they are confined by 
a potential well. Two-dimensional electron gases in semiconductors can be 
classified into two groups. First, a 2DEG can be formed at the surface of some 
semiconductors, e.g. InAs. Second, a 2DEG can be located at the interface 
between two different semiconductor materials, e.g. GaAs and AlGaAs. Since 
both kinds of structures are employed in superconductor/two-dimensional- 
electron-gas hybrid structures, the properties of both will be discussed briefly, 
below. Detailed information on 2DEGs can be found, for example, in [70, 71]. 

3.1.1 Two-Dimensional Electron Gases in Inversion Layers 

In some semiconducting materials two-dimensional electron gases are formed 
directly at the surface of the crystal. Here, a downward bending of the con- 
duction band at the surface leads to a potential well in which the electrons are 
trapped. This downward bending has its origin in the modified atomic bonds 
at the surface. Owing to the broken periodicity of the crystal, the electron 
structure at the surface differs considerably from the structure in the bulk 
semiconductor. The free valences left at the surface atoms are saturated ei- 
ther by a reconstruction of the surface structure or by adsorption of atoms or 
molecules. In order to describe these modified surface bonds, discrete energy 
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levels, so-called surface states, within the band gap of the semiconductor are 
introduced [72]. The downward bending of the valence band at the surface 
is a result of the balancing of the Fermi level between the charged surface 
states and the space-charge region in the semiconductor crystal. 

A surface two-dimensional electron gas, as described above, is found in 
p-doped InAs [73, 74]. Here, the filling of surface states leads to a depletion of 
holes near the surface and thus to the downward bending of the valence and 
conduction bands (Fig. 3.1) [75]. Directly at the surface, an inversion layer 
with accumulated electrons is formed. Inversion means that, in contrast to the 
bulk crystal, where holes are the majority carriers, the majority carriers are 
now electrons. Between the inversion layer and the bulk material, a depletion 
region exists. The electrons are spatially confined in the inversion layer and 
thus cannot move perpendicular to the surface. The width of the inversion 
layer is in the order of 10 nm and therefore comparable with the Fermi 
wavelength. As a result the electron energy is quantized in discrete levels 
(Fig. 3.1). Since free electron motion is only possible along the surface, a 
two-dimensional electron gas is formed. 




Fig. 3.1. Schematic band profile of p-doped InAs at the surface. Directly at the 
surface, an inversion layer is formed. The electrons are confined in a quantum well 
with discrete energy levels Ei and E 2 



The basic properties of a two-dimensional system can be derived from 
the Schrodinger equation. For a confinement potential eV{z) along the 2 ; 
direction, as shown in Fig. 3.1, the Schrodinger equation can be written as 

eV(z) 'ip{r) = Eifj{r) . (3.1) 



2m" 



/ 



dy"^ dz'^ 



Here, m* is the effective electron mass in the semiconductor, which can differ 
considerably from the free electron mass. Owing to the free motion in the xy 
plane, the following ansatz can be used for the wave function: 
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^(r) = (fi{z) exp{ikxX + ikyy) = (p{z) exp(ifc|| • r) . (3.2) 

The Schrodinger equation can thus be separated into an equation for the free 
motion along fe|| in the xy plane, 

^ exp{ikxX + ikyy) = exp{ik^x + ikyy) , (3.3) 

and an equation for the bound states within the potential well, 

Ei, where i = 1,2,..., are the discrete energy levels of the quantum well. As 
depicted in Fig. 3.2, the energy dispersion relation has the form of parabolas, 
representing the kinetic motion within the xy plane, shifted by the confine- 
ment energy Ei. The density of states of each subband Di{E) is constant, 

A(A) = ^ , (3.5) 

ttTi 

so that the total density of states increases stepwise each time the next energy 
level Ei is reached (Fig. 3.2). Owing to the small band gap of InAs, a non- 



2m* dy'^ 




Fig. 3.2. Energy dispersion relation of a two-dimensional electron gas. The parabo- 
las of the the kinetic energy are shifted by the energy of the qnantized levels in the 
quantum well. The corresponding density of states is shown on the right 



parabolic energy dispersion has to be taken into account at higher energies. 
In this case Di{E) is not constant anymore [76]. 

For the case when only one subband is occupied, the sheet electron con- 
centration ri 2 D can be calculated by using (3.5): 

«2D = / Di{E)dE = El) , (3.6) 

J El 7 tTi 

where y is the Fermi energy. Since y — Ei = ?i^fcpj^/2m*, the Fermi wave 
vector of the two-dimensional electron gas can be expressed as 
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fcpN = 



m fpN 

n 



= V27rn2D ■ 



(3.7) 



The confinement potential for a two-dimensional electron gas in InAs can 
be approximated by a triangular well [70]. Typical electron concentrations 
are in the order of ri 2 D = 10^^ cm“^ [74]. At low temperatures the electron 
mobility fie is limited by scattering on ionized impurities and by surface 
roughness scattering. A typical value at low temperature (4.2 K) is fig = 
10 000 cmVV s [74, 77]. 

The elastic mean free path /ei for carriers in the two-dimensional electron 
gas is given by 



lei = r'FNTel . 



(3.8) 



Here, Tei is the elastic scattering time, which can be extracted from the mo- 
bility by using Tei = m* fieje.- Taking the typical values for InAs of fig and 
ri 2 D results in an elastic mean free path of less than one micrometer. 

As we shall see in Sect. 3.2.1, low-resistance contacts to the two-dimen- 
sional channel can be prepared by evaporating a superconductor directly on 
the surface [8]. In order to improve the contact resistance the semiconductor 
surface is usually cleaned prior to the metal deposition, e.g. by Ar sputter- 
ing. The two-dimensional electron gas on the surface of p-InAs is difficult to 
suppress, with the consequence that adjacent structures are not sufficiently 
isolated. Isolation can be gained by growing an additional AlInAs layer on 
top of the p-InAs layer by molecular beam epitaxy [32]. Owing to the band 
offset between AlInAs and InAs, the band bending in InAs is reduced so that 
no electrons are accumulated. By simply removing the AlInAs layer by wet 
chemical etching, the two-dimensional electron gas in the inversion layer can 
be recovered. 



3.1.2 Two-Dimensional Electron Gases in Heterostructures 

In a semiconductor heteroj unction two different semiconducting materials 
with different band gaps are combined. The band profile for a typical het- 
erostructure is shown in Fig. 3.3. Here, GaAs is used as the lightly n-doped 
material with the smaller band gap, while the highly n-doped AlGaAs pos- 
sesses the wide band gap. At the interface the conduction band Eg and va- 
lence band change abruptly on an atomic length scale by AEg and AE^, 
respectively. Furthermore, the alignment of the Fermi level leads to a charge 
transfer between the two materials. As a result electrons are trapped at the 
interface within the lower-band-gap semiconductor. The electrons are sup- 
plied by the donors of the wide-gap material. Similarly to the inversion layer 
at the InAs surface, the electrons can move freely along the interface within 
the GaAs layer. The scattering of electrons on ionized impurities can be sup- 
pressed almost completely if the electrons in the well are separated from the 
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Fig. 3.3. Two-dimensional electron gas in a GaAs/AlGaAs heterostructure. The 
electrons are accumulated at the GaAs/AlGaAs interface in the low-gap material 



ionized donors by an imdoped spacer layer. By employing this so-called mod- 
ulation doping technique, the electron mobility is increased and thus very 
long electron mean free paths lei can be obtained [78]. 

Electron mobilities as high as 10^ cm^/V s have been reported for two- 
dimensional electron gases in GaAs/AlGaAs modulation-doped heterojunc- 
tions [79, 80]. In principle, these structures would be an ideal basis for ballistic 
superconductor/2DEG structures owing to the extremely large elastic mean 
free path of several micrometers. “Ballistic” means that the carriers can pass 
through the structure without any scattering. However, a serious obstacle 
concerning the coupling of a superconductor to the 2DEG in a GaAs/AlGaAs 
heterostructure is the high Schottky barrier at the interface. Owing to this 
barrier, the carrier transfer from the metallic superconductor into the 2DEG 
is considerably hindered. Nevertheless, superconductor/2DEG contacts based 
on GaAs/AlGaAs have been successfully prepared by using either annealed 
indium or tin contacts or by epitaxial growth of aluminum [o5, 81, 82]. 

With regard to a high interface transparency, heterostructures where the 
two-dimensional electron gas is located in an InAs layer are ideal candidates, 
since in this case no Schottky barrier is formed. In these structures, AlSb or 
Alo.2Gao.8Sb is used as a barrier material for the quantum well (Fig. 3.4a) [83, 
84, 85]. Typical electron concentrations and mobilities at low temperatures 
(< 50 K) are in the order of 1 x 10^^ cm“^ and 60 000 cm^/V s, respectively 
[84]. In order to make contact between the 2DEG and the superconductor, 
the heterostructure has to be etched down to the top of the InAs channel and 
covered subsequently by the superconductor layer. By this method a larger 
contact area and thus a low interface resistance are obtained [28, 36, 86, 87] . 
Alternatively, the layer system can be etched below the InAs channel so that 
the superconductor layer is contacted at the sidewall of the semiconductor 
mesa [23]. The latter method leads to a better geometrical definition of the 
contact at the cost of a larger interface resistance. 

Heterostructures with an InAs-based channel layer can also be realized 
by using Alo.4sIno.52As as a barrier material if the structure is grown by 
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Fig. 3.4. (a) Typical layer structure of an AlSb/IuAs lieterostructure (on a GaAs 
semi-insulating substrate). The lower barrier is formed by an AlSb/GaSb super- 
lattice terminated by AlSb. The top GaSb cap layer prevents oxidation of the 
AlSb layer [88]. (b) A two-dimensional electron gas in an AlInAs/In(Ga)As/InP 
heterostructure. The electrons are located in the InGaAs/InAs/InGaAs sandwich. 
The carriers are supplied by the AlInAs n+ dopant layer [89] 



molecular beam epitaxy (MBE), or by using an InP barrier if low-pressure 
metal-organic vapor phase epitaxy is employed. However, the large lattice 
mismatch between the channel and the barrier restricts the maximum al- 
lowed thickness of the InAs layer before strain causes dislocations in the 
crystal. In order to obtain a sufficiently low bound state so that electrons are 
accumulated in the quantum well, the quantum well width can be extended 
effectively by sandwiching the InAs layer between two Ino. 53 Gao. 47 As layers 
(Fig. 3.4b) [89, 90]. The latter is lattice matched to InP. Another feasible way 
is to use a strained InGaAs channel layer with a higher indium content, with 
respect to Ino. 53 Gao. 47 As [39, 40]. By this method the quantum well can be 
grown sufficiently thick to accumulate enough electrons. For a channel layer 
with an indium content of 77%, very high mobilities of up to 4.5 x 10® cm^/V 
s at a sheet carrier concentration of 6.0 x 10^^ cm“^ were obtained. The layer 
sequence of this type of heterostructure is shown in Fig. 4.1 in Sect. 4.1.1. 



3.2 Fabrication of Superconductor/ 
Two-Dimensional-Electron-Gas / 

Superconductor Junctions 

The different locations of the two-dimensional electron gas in a surface in- 
version layer and in a heterostructure require different methods to con- 
tact the 2DFG to a superconducting electrode. Only the basic princi- 
ples concerning the fabrication of superconductor/two-dimensional-electron- 
gas/superconductor junctions will be discussed in the following two sections. 
Of course, fabricating these elements is only a first step for the experiments 
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described in the subsequent chapters, and additional fabrication processes 
are required. 

3.2.1 Junctions Based on a Surface Inversion Layer 

In the case of a semiconductor structure with an inversion layer, the supercon- 
ductor can be deposited directly on the surface in order to make contact to the 
two-dimensional electron gas. Josephson coupling can be obtained if the two 
superconducting electrodes are linked by a small semiconducting bridge with 
a length of a few hundreds of nanometers. One possible way to obtain such 
small dimensions is to use electron beam lithography [8, 14, 91]. However, 
the use of electron beam lithography can be circumvented by using simple 
optical lithography in conjunction with some special preparation techniques 
in order to obtain a small separation of the superconducting electrodes. 

The first method is based on a step edge, which is used as a shadow mask 
for the definition of the superconducting contacts [30, 92]. The concept of a 
step edge junction was first introduced by Serfaty et al. [93]. A schematic illus- 
tration of this type of junction based on p-type InAs is depicted in Fig. 3.5a. 
In the first preparation step the upper Nb electrode is defined by optical 
lithography and lift-off. Prior to the deposition, the semiconductor surface is 
cleaned either by a water solution of NH3 or by Ar+ sputtering to remove 
oxides from the surface. After the first Nb layer is deposited, the step, which 
defines the channel length, is wet chemically etched (H202:H3P04:H20), us- 
ing the first Nb layer as an etch mask. Typical step heights are around 100 nm. 
Next, the second Nb layer is deposited. In order to avoid shorts between the 
electrodes, the sample is tilted during the deposition. The two-dimensional 
channel at the InAs surface follows the step from the top Nb electrode to the 
bottom electrode. In order to obtain a well-defined two-dimensional electron 
gas, a low-p-doped (< 10^^ cm“^) InAs layer is grown by molecular beam epi- 
taxy on an InAs substrate, which is usually highly p-doped. By this method, 
the two-dimensional channel is sufficiently insulated from the substrate. 

As an interesting feature, additional semiconducting layers embedded in 
the InAs can be inserted by using molecular beam epitaxy. These layers 
can, for example, serve as barriers in the channel along the step [20]. In 
other structures, a p-doped layer used to form the channel at the surface is 
sandwiched between two n-doped layers, which serve as a contact layer to the 
Nb electrodes [30]. 

Another very interesting method to fabricate Nb-InAs-Nb junctions was 
developed by Chrestin et al. [29, 59]. The structure is shown in Fig. 3.5b. 
After the first Nb electrode has been defined by sputter deposition and lift- 
off, the front edge of the Nb electrode is anodically oxidized. The insulating 
Nb20s layer serves as a spacer for the second Nb layer. The geometry for 
this second Nb layer is defined by optical lithography. Before the Nb film 
can be deposited, the oxidized InAs has to be removed by a selective wet 
chemical etch (H2S04:H20), which does not etch the Nb205 layer. Owing to 
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Fig. 3.5. Two alternative methods which can be used to prepare Nb-InAs-Nb 
junctions, (a) Step edge junction, where the Nb electrode separation is determined 
by a step etched into the InAs substrate, (a) Nb-InAs junction, where the channel 
length is given by an anodic oxidized Nb layer 



the overlapping Nb layers separated by the Nb20s layer, a relatively large 
capacitance is connected in parallel with the junction, which has a large 
influence on the current- volt age characteristics. 



3.2.2 Junctions Based on 
a Buried Two-Dimensional Electron Gas 

In order to prepare junctions with a buried two-dimensional electron gas in 
a heterostructure, electron beam lithography is mostly employed in order to 
obtain a superconductor electrode separation of a few hundred nanometers. 
Nevertheless, advanced self-aligned optical lithography processes can also be 
used to obtain electrode separations of less than half a micron [94]. 

In order to give an impression of how a typical junction with a buried 
two-dimensional electron gas is fabricated, the preparation steps required for 
a Nb-InGaAs/InP-Nb junction are described in detail. For the fabrication 
of these junctions, a two-step electron beam lithography process is employed. 
The junctions are based on the heterostructure introduced in Sect. 3.1.2. 
In the first step an InGaAs/InP wire structure is defined by electron beam 
lithography and subsequent reactive ion etching (RIE) with GH 4 /H 2 , as de- 
picted in Fig. 3.6. A Ti layer is used as an etching mask to fabricate wire 
structures having a width of a few hundreds of nanometers. 

After the etching procedure, the samples are mounted in an ultrahigh- 
vacuum chamber, where Ar+ sputtering is used to clean the semiconductor 
surface. Subsequently, a 100 nm thick Nb layer is evaporated in situ. De- 
positing the Nb layer directly on the semiconductor and thus avoiding any 
lithography step after the RIE process ensures that no residual organic con- 
tamination from the resist remains on the surface. 

The geometry of the Nb electrodes is defined by a second electron beam 
lithography step aligned with the first pattern. Again a Ti etching mask 
is employed for the reactive-ion-etching process (SFg) of the Nb layer. A 
scanning electron micrograph of a complete sample is shown in Fig. 3.7. It 
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Fig. 3.6. Preparation steps used for the fabrication of a Nb-InGaAs/InP-Nb junc- 
tion: (a) CH 4 /H 2 reactive ion etching of the semiconductor wire structure, (b) Ar^ 
sputter cleaning before Nb deposition, (c) Nb deposition, (d) geometrical definition 
of the Nb electrodes by electron beam lithography and SFg RIE 




Fig. 3.7. Scanning electron micrograph of a typical Nb-InGaAs/InP-Nb junction 
with a length L of 450 nm and a width IF of 2 | 0 ,m [165] 



can be seen that the InP surface is covered by In droplets. During the Ar+ 
sputtering, the phosphorus component is removed from the lattice so that the 
remaining surplus indium forms droplets on the surface [95, 96]. In order to 
prevent a conductive bypass channel on the surface, the samples are etched 
in an acetic acid solution [9 i, 97]. 

As an alternative to the fabrication methods described above, the width 
of the semiconducting channel can be defined by wet chemical etching [38, 98] 
or by wet chemical etching combined with Ar’*' milling [14]. For some types of 
junctions, the semiconductor is only etched down to the top of the conductive 
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layer [38, 98]. By this technique, the contact area between the superconductor 
and the semiconductor is enlarged. The definition of the channel by etching 
can be completely omitted if the superconducting electrodes are fabricated 
by alloying, as has been demonstrated for AlGaAs/GaAs-based junctions 
[35, 81]. 



4. Transport Studies 
on Single Superconductor/ 
Two-Dimensional-Electron-Gas Interfaces 



A detailed knowledge of single superconductor/two-dimensional-electron-gas 
structures is very important in order to understand more complex systems 
like the S/2DEG/S junctions treated in the following chapters. First, we shall 
focus on the effect of different surface treatments prior to the deposition of 
the Nb electrodes. Measurements of the differential resistance as a function of 
the voltage drop at the interface give direct information about the Andreev 
reflection probability and thus about the transparency of the interface. Only 
two-dimensional electron gases in heterojunctions are covered here. However, 
most of the findings also apply to inversion layer structures. Several alterna- 
tive fabrication methods for S/2DEG contacts which have been investigated 
and evaluated will be reported. By using a point contact to inject carriers 
into a superconductor/two-dimensional-electron-gas interface, the Andreev 
reflection probability can be measured, directly. In addition, point contacts 
can also be used to obtain information about the angular distribution of the 
Andreev reflection. 



4.1 Effect of the Surface Treatment 

The transparency of a superconductor/2DEG interface depends crucially 
on the surface treatment of the semiconductor prior to the deposition of 
the superconductor layer. Here, different cleaning methods for InGaAs/InP 
heterostructures will be compared. For the electrical characterization of 
the differential conductance, dl/dU was measured by using a lock-in tech- 
nique. This approach has been chosen earlier also for other superconduc- 
tor/semiconductor structures e.g. for Nb-Si [99, 100], Nb-Ino. 53 Gao. 47 As 
[33, 101, 102], Sn-GaAs/AlGaAs [103] and Nb-(Al,Ga)Sb systems [ 86 , 104]. 

4.1.1 Preparation of Single S/2DEG Contacts 

In order to study of the transparency of an S/2DEG interface prepared by 
different methods, it is sufficient to fabricate relatively large samples by op- 
tical lithography [45]. In contrast to the 2DEG at the surface of InAs, the 
2DEG in a heterostructure is covered by a nonconductive cap layer. There- 
fore, it is necessary to contact the two-dimensional channel from the side. For 
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the following studies InGaAs/InP structures were used. Of course, the results 
reported here can be transferred to the preparation of superconductor/two- 
dimensional-electron-gas structures based on a different type of heterostruc- 
ture. 

The Ino. 53 Gao. 47 As/Ino, 77 Gao. 23 As/InP heterostructure was grown by 
low-pressure metal-organic vapor phase epitaxy (LP-MOVPE) [39, 40]. In 
Fig. 4.1, the layer sequence together with a schematic conduction band pro- 
file is shown. The electrons of the 2DEG, which is located in the strained 
Ino. 77 Gao. 23 As layer, are supplied from a 10 nm thick InP dopant layer, which 
is separated by a 20 nm thick InP spacer layer. The dopant concentration was 
4.9 X 10^^ cm“^. On top of the strained Ino. 77 Gao. 23 As channel, a 150 nm thick 
Ino. 53 Gao. 47 As layer, lattice matched to InP, was grown. Since the quality of 
the InGaAs/InP interface is not as good as that of the inverse InP/InGaAs 
interface if grown by MOVPE, Ino. 53 Gao. 47 As instead of InP was used as a 
top barrier in order to obtain maximum electron mobilities. Shubnikov-de 
Haas and quantum Hall measurements at 1.4 K yield mobilities of about 
370 000 cm^/V s and a carrier concentration of 7 x 10^^ cm“^. An effective 
electron mass of m* = 0.04 We was extracted from temperature-dependent 
Shubnikov~de Haas measurements. 
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Fig. 4.1. Layer system of the Ino.53Gao.47As/Ino.77Gao.23As/InP heterostructure. 
The conduction band profile is drawn schematically 



In order to contact the 2DEG laterally to the superconductor, the het- 
erostructure is partially removed by a GH 4 /H 2 reactive-ion-etching process 
(Fig. 4.2). Subsequently, alloyed Ni/AuGe/Ni ohmic contacts are prepared 
in order to measure the voltage drop at the interface between the supercon- 
ductor and the 2DEG. The crucial step for the preparation of transparent 
Nb/2DEG interfaces is the cleaning procedure, which has to be carried out 
immediately before the evaporation of the Nb layer. Either wet chemical 
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etching, ion sputter cleaning or a plasma treatment can be performed for 
this purpose. Details of the surface cleaning procedure are described in con- 
junction with the corresponding measurements. The 100 nm thick Nb film is 
deposited by sputtering or by electron beam evaporation. The width of the 
contact is TT = 100 | 4 m. The Nb film has an overlap of a few microns with 
the Ino. 53 Gao. 47 As cap layer. However, the possible leakage current can be 
neglected, as shown by a comparative study of top and side contacts [105]. 



lOOpm 






Fig. 4.2. Schematic top view and cross section of the lateral contact from the 
Nb layer to the 2DEG. In order to obtain a homogeneous current distribution, the 
current flows from one Nb electrode to the opposite one 



4.1.2 Samples Cleaned by Wet Chemical Etching 

Wet chemical cleaning has an advantage over ion sputter methods because the 
semiconductor crystal is not degraded by the impact of ions. Especially for 
InP, this can be a serious problem [95, 96, 106]. For the wet chemical cleaning, 
the samples were etched in 1 HE : 20 H 2 O [45]. HE was chosen since oxides 
are expected on the surface owing to an O 2 treatment after the CH 4 /H 2 RIE 
and owing to exposure of the structures to air before the Nb deposition. As 
an alternative to HE, NH 4 OH [30] or HCI/CH 4 COOH/H 2 O 2 /H 2 [107] can 
be used to remove oxides from the surface before Nb deposition. For the 
HF-cleaned samples reported here, the Nb film was deposited by sputtering. 

The normalized differential conductance R^dl /AV of the sample cleaned 
by HE wet chemical etching is plotted as a function of eV j A q in Fig. 4.3. 
The normal-state resistance for this sample was 1 kQ and Z\q had a value of 
1.38 meV, corresponding to = 9.1 K. The experimental curve reveals two 
evident features. First, ctsn(O) = d//dV(v=o) is reduced to only 14% of its 
normal-state value cnn- Second, maxima are observed near eV = AzAq. From 
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the drop of the differential conductance at low bias voltage, it can be con- 
cluded that a barrier is present at the interface. This barrier could partially 
originate from the Fermi velocity mismatch between the superconductor and 
the 2DEG. Furthermore, adsorbates on the surface such as fluoride and water 
could give rise to an additional contribution to the barrier. 




eV/A„ 



Fig. 4.3. Normalized differential conductance as a function of eF/Zio (o). The 
BTK curve for Zbtk = 2 is plotted for comparison [dashed line). The ht assuming 
an SN/2DEG layer system with 71 = 1.6, 72 = 0.2, 7B2 ~ — 6, d^ = 0.06 

and /sN = 0.07 is plotted as a solid line [45] (© (1996) by the American Physical 
Society) 



The overall shape of the i?NdI /dF curve can be explained by the BTK 
model if a finite barrier [Z > 0) is assumed at the interface. Owing to normal 
reflection processes, a decrease of the differential conductance is expected for 
|eF| < Aq, while the singular density of states near Aq causes an enhance- 
ment of d//dF. A Z parameter of 2.0, which gives rise to the observed drop 
of ctsn(O), was determined from crsN(0)/crNN = 2(1 -I- Z'^)/(\ -|- 2Z^)^ [46]. 
As can be seen in Fig. 4.3, the curve calculated from the BTK model does 
not agree very well with the measured one. The peaks in the BTK model are 
found at a larger bias voltage and have a substantially higher peak value. 

Obviously, the BTK model does not describe the experimental findings 
sufficiently. The shift of the conduction maxima towards lower voltages can 
be explained if a normal-conducting film (N) between the Nb layer (S) and 
the 2DEG is assumed [45]. This normal-conducting film consists most likely 
of NbOa; compounds. Owing to the strong gettering of the Nb film, the first 
atomic layers in particular contain oxygen, which originates from residual gas 
atoms in the sputter chamber and also from water molecules which resided 
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on the surface after the wet chemical treatment. NbO^ {x ~ 1) has metallic 
properties but is not superconducting [108]. 



L =0 X 
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Fig. 4.4. (a) Schematic illustration of the Cooper pair density F{x) and (b) po- 
tential drop eip{x) in an SN/2DEG layer sequence. 71,2 describe the suppression of 
F{x) at the SN and N-2DEG interfaces, while 7 bi 2 S’l'® ^ measure of the disconti- 
nuity of F{x) at the corresponding interfaces 



Since it seems to be reasonable to include an additional normal-conducting 
layer, the extended model introduced in Sect. 2.2.3, which incorporates the 
proximity effect in the SN layer, can be applied. The slope of the Cooper pair 
density F{x), as well as the electrical potential eif{x), is depicted schemat- 
ically in Fig. 4.4. Owing to the proximity effect, F{x) decreases in the N 
layer. The lower value of F{x) at the N/2DEG interface in comparison with 
its value in the bulk superconductor is responsible for the shift of the con- 
ductance maxima towards lower voltages. No barrier was assumed between 
the S and N layers (7B1 = 0), because of the good electrical contact between 
the two metals. Owing to the comparable conductances of N and S, 71 must 
be in the order of one. Since the N layer is connected to a clean 2DEG, the 
definition of 72, as given by (2.30), cannot be applied directly. However, 72 
can be estimated from 72 « (3fcF.2DEG/XFN^ei,N)CN/^2DEGi as shown in [109]. 
?2DEG = (?i^/fcBTm*)-\/n2D/27r is the coherence length of the 2DEG in the 
clean limit. For a 2DEG as used here, the upper limit for 72 can be set at 
0.2 [15]. As can be seen in Fig. 4.3, the best fit to the measured data was 
obtained by assuming a barrier at the N/2DEG interface with 7B2 = 6. If the 
clean limit applies, this value is connected directly with the Z parameter of 
the BTK model by the relation 7B2 = .^^, so that Z w 2.5 [45]. It can thus 
be concluded that the barrier strength needed to obtain a good fit is larger 
than that estimated from the BTK model {Z = 2), with the consequence 
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that the amount of Andreev reflection is overestimated. Even if the shift of 
the maxima in the differential conductance can be explained by the model 
given above, inelastic scattering, represented by Isn, has to be incorporated 
additionally in order to explain the damping of the conductance peaks. For 
this structure an inelastic mean free path in Nb of 1.4 pm was taken, which 
was extrapolated from a value given for Nb at 10 K [110]. 

The experimental curve is fitted well by the model which considers a 
proximity layer (N). By using the BTK model and the extended model, a 
high barrier at the interface is found. In addition, the normal-state resistance 
i?N = 1 kQ is also very large, which can be explained by an inhomogeneous 
Nb-2DEG interface. 

4.1.3 Ar+ Sputter Cleaning 

As an alternative to wet chemical etching, the surface of the semiconductor 
can be cleaned by Ar+ sputtering. For the set of structures discussed first, 
RE Ar+ sputtering (800 V, 90 s) was performed in the same chamber where 
the 100 nm thick Nb film was deposited subsequently to sputtering [15]. The 
resulting normalized differential conductance as a function of bias voltage 
is given in Fig. 4.5. Two striking differences can be made out. First, the 
drop of Rjsidl /dV is significantly lower than for the wet-chemically-cleaned 
sample. Second, the conductance maxima are damped out more strongly. 
Even more important is the normal-state resistance, which is, with a value of 
only 8 Q, two orders of magnitude smaller than i?N for the wet-chemically- 
cleaned sample. 

The ratio of the differential conductances crsN(0)/crNN can be used to de- 
termine the Z parameter within the BTK model. For this sample Z = 0.8 
is obtained, which is significantly lower than the value calculated for the 
wet-chemically-cleaned sample. However, the general shape of the measured 
differential-resistance curve does not fit at all to the BTK curve for Z = 0.8. 
In order to obtain a more realistic description of the experimental situation, 
the degradation of the semiconductor surface by the Ar+ ions must be in- 
cluded in the model [45]. 

As depicted in Fig. 4.6, the previously presented model has to be ex- 
tended by an additional diffusive N' layer, which represents the degraded 
region in the 2DEG. The semiconductor has lost its ballistic properties ow- 
ing to crystal defects near the surface [95, 96, 106]. The inclusion of the 
nonballistic N' layer makes it necessary to assume a voltage drop not only at 
the metal/semiconductor interface but also within the N' layer itself. A the- 
oretical model, which accounts for a nonequilibrium proximity effect in the 
N' layer, was developed by Volkov et al. [111]. An extended model including 
this theory, which applies to the SN/N'2DEG layer structure, is described in 
detail in [45]. The major difference compared with the model used for the 
SN/2DEG structure is that now the current through the structure has to 
be calculated from the transmission probability T{E) through the N' layer 
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Fig. 4.5. Normalized differential conductance as a function of eV/Ao with Aq = 
1.38 meV (•). The curve calculated using an SN-N'2DEG model with 7mi = 0.15, 
7 b 2 , ^N' = 7 ^N' and 7^/ = Esn = 0.26 is plotted as the solid line. Fjqi and Esn are 
the parameters for the inelastic scattering in the corresponding layers [45] (© (1996) 
by the American Physical Society) 



and not of the coefficients for Andreev reflection A{E) and normal reflection 
B{E). This is necessary since here the voltage drop in the N' layer must be 
taken into account. The set of parameters which leads to a good fit to the 
experimental data is given in Fig. 4.5. Details of the simulation are given in 
[ lo] . Within this model, the parameter 7B2 = provides only a lower bound 
for the Z parameter. Here, a Z parameter of 2.3 is obtained, which is consid- 
erably larger than the corresponding value for the BTK model (2’btk = 0.8). 

For the Ar+-sputtered (800 V) sample it can be concluded that, although 
the BTK model suggests a lower barrier at the interface, the actual barrier 
is comparable to the barrier obtained for wet-chemically-cleaned structures. 
The underestimation of the Z parameter in the BTK model is due to the 
fact that the voltage drop in the degraded semiconductor layer is not in- 
cluded. Nevertheless, Ar+ sputter cleaning will be one of the preferred clean- 
ing methods for Nb-2DEG junctions, since the normal-state resistance of 8 Q 
is significantly lower than Rjq for wet-chemically-cleaned samples. 

In order to reduce the surface degradation, an Ar+ ion beam source with 
a lower acceleration voltage (300 V) can be employed [112]. In Fig. 4.7, the 
normalized differential conductance of samples cleaned for 1 h, 3 h and 6 h 
with a low ion current of 2 pA is plotted as a function of the voltage drop V. 
The most striking feature is that, with increasing sputtering time, the zero- 
bias conductance increases, while the maxima near EAq are damped out 
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Fig. 4.6. (a) Schematic illustration of the Cooper pair density F{x) and (b) the 
potential drop eip(x) in an SN/N'2DEG layer system. 71^2 is a measure of the 
suppression of F{x) at the SN and NN' interfaces. The barrier at the NN' interface 
is described by 7B2 ■ At the SN interface, no barrier is assumed (7B1 = 0) 



more strongly. Apart from that, the normal-state resistance also decreases 
with increasing sputtering time. For the samples sputtered for 1 h, 3 h and 6 h 
normal-state resistances of 60 Q, 13 Q and 8 Q, respectively, were measured. 
The set of samples sputtered for 1 h showed a large spread of i?N between 
32 Q and 500 Q [H2]. This is presumably due to an incomplete removal of 
residual atoms, which therefore leads to a large variation of the transmission 
probability. For the samples sputtered for 3 h and 6 h, this was not the 
case. Here, a variation in the range of 10% was measured. For the sample 
sputtered for 1 h, Zbtk has a value of 1.14, while values of 0.7 and 0.56 
were obtained for the structures sputtered for 3 h and 6 h, respectively. For 
increasing sputtering times the i?Nd//d% vs. V curves resemble more and 
more the damped curve measured for the 800 V Ar+-cleaned samples. At 
least partially, this tendency can again be attributed to the formation of a 
degraded semiconductor layer at the surface. 

An interesting question is to what extent the 2DEG is degraded by the 
Ar’*' sputter cleaning. This information is especially important for the Joseph- 
son junctions discussed in the following chapters, where the Nb electrodes are 
only a few hundreds of nanometers apart. In order to address the problem 
of degradation by Ar+ sputtering, wire structures were prepared. The un- 
derlying idea is that owing to the small width of the wire, the effect of Ar+ 
sputtering on the 2DEG should be directly visible in magnetotransport mea- 
surements. In order to separate this effect from other effects a comparative 
study was performed, where the voltage drop at the Ar+-ion-treated part of 
the wire was compared with the voltage drop at a nontreated section of the 
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Voltage (mV) 



Fig. 4.7. Normalized differential conductance as a function of bias voltage at T = 1 
K for samples cleaned by Ar+ (300 V) sputtering. The cleaning times were 1 h {solid 
line), 3 h (dash- dotted line) and 6 h (dashed line) 



wire. Selective sputtering was performed by opening a window in an electron 
beam lithography resist layer. The covered area of the semiconductor is not 
affected by the Ar+ sputtering. The 400 nm wide wires were prepared by 
electron beam lithography and CH4/H2 reactive ion etching. A schematic of 
the geometry and measurement setup can be found in Fig. 4.8b (inset). 

In Fig. 4.8a, the magnetoresistance Rxx of a wire cleaned with Ar+ for 
3 h is plotted together with R^x of the reference wire. The latter curve shows 
pronounced Shubnikov-de Haas oscillations for i? > 1 T. The decrease of 
Rxx{B) for H < 1 T can be explained by diffusive interface scattering in 
two-dimensional wire structures [113]. The wire cleaned for 3 h contains a 
conductive bypass channel. This is deduced from the lower zero-field resis- 
tance Rxx{B = 0) in comparison with the nonsputtered wire and from the 
overall increase of Rxx with B. The local maximum at i? « 0.3 T can be 
attributed to diffusive interface scattering [113]. 

From the magnetoresistance measurements, it can be concluded that the 
two-dimensional character of the transport in the InGaAs/InP wire is not 
destroyed by the Ar+ sputter cleaning, because the oscillations in Rxx still 
remain after Ar+ sputtering. However, the sputtered part of the wire pos- 
sesses a higher electron concentration, which can be attributed to an altered 
pinning of the Fermi level at the surface due to the sputtering. The bypass 
channel can be explained by conduction in the sputtered semiconductor sur- 
face. Especially on an InP surface, Ar+ sputtering has the effect that the 
phosphorus component is removed selectively, leaving a rough. In-rich sur- 
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Fig. 4.8. (a) Magnetoresistance Rxx of the wire section cleaned for 3 h by Ar+ 
ions and of the non-sputtered wire section, (b) Magnetoresistance after a wet- 
chemical-etching step using HF. The inset shows the geometry of the sample and 
the measurement configuration 



face, as shown in Fig. 4.9. On exposure to air, a conductive Iu 203 layer is 
formed, which is measured as a bypass channel. From the observation that the 
conductance of the bypass channel increases with increasing sputtering time, 
it can be inferred that the oxidation process is supported by the Ar+-induced 
roughness of the surface [95]. For InGaAs, the oxidation is probably much 
less pronounced, since here the surface is more stable against Ar+ sputtering 

[ 114 ]. 

The proof that oxidation of the surface is responsible for the conductive 
bypass was delivered by a subsequent HF etching of the wires. Fig. 4.8b shows 
the magnetoresistance of the sample after an HF dip, which was performed 
immediately before the wire structure was mounted into the cryostat again. 
Obviously, the bypass channel of the Ar+-sputtered wire section has vanished, 
because the overall resistance does not increase with B anymore. In addition, 
the zero-field resistance is now almost the same for the Ar^-sputtered and 
nonsputtered sections of the wire. For structures sputtered for 6 h at 300 V 
and sputtered for 90 s at 800 V, the bypass channel could not be removed 
completely by the HF dip. Possibly the degradation of the semiconductor sur- 
face is in this case deeper than the oxidized layer. Only the latter is removed 
by HF. 

As shown above, the study of the wire structures gives valuable informa- 
tion concerning the choice of the cleaning procedure. The measurements on 
Ar+-cleaned wires show that the 2DEG in the wire is effectively not degraded 
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Fig. 4.9. Scanning atomic force microscope picture of an InP surface sputtered for 
6 h with 300 V Ar+ ions. On the surface In droplets are found, which are formed 
after a selective removal of the phosphorus component 



by sputtering. This is a crucial result concerning the S/2DEG/S Josephson 
junctions. The InP surface is especially affected by the cleaning procedure, 
which can be deduced from the conductive oxide layer. Of course, in case of in- 
situ cleaning and Nb deposition, this oxidation does not occur. Nevertheless, 
the degraded InP affects the transmissivity through the Nb/semiconductor 
interface in any case. The oxidation of the surface can also be prevented if 
the sample is passivated directly after etching of the semiconductor mesa, by 
using (NH 4 ) 2 S [115]. Here, one or two monolayers of sulfur cover the surface. 
If a brief Ar’*' sputter cleaning is employed, this coverage can be removed 
before the Nb layer is deposited. 

4.1.4 Electron Cyclotron Resonance Plasma Cleaning 

Although Ar+ cleaning is a viable way to produce transparent supercon- 
ductor/semiconductor junctions, it is still worthwhile to search for more gen- 
tile cleaning procedures. A very promising alternative cleaning method is an 
electron cyclotron resonance plasma treatment of the semiconductor surface, 
which is based on the well-established cleaning method prior to the regrowth 
of semiconductor layers by molecular beam epitaxy [116, 117, 118]. In most 
cases a hydrogen plasma is used, which removes surface oxides, hydrocar- 
bons and other contaminants. However, it was shown for a phosphorus-based 
HI-V semiconductor, i.e. InP, that a hydrogen electron cyclotron resonance 
plasma can lead to a group-HI-enriched surface owing to the formation of 
P hydrids [118]. Therefore, a mixture of He/H or even pure He was used for 
the cleaning of the InGaAs/InP surface. The underlying idea is that owing 
to the lower mass of the helium or hydrogen atoms compared with Ar atoms, 
the momentum transfer at the semiconductor surface is considerably smaller. 
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After the CH4/H2 etching process for the definition of the semiconductor 
mesa, the samples were treated in the He plasma for 2 h and subsequently 
completely covered in situ by a 100 nm thick electron-beam-evaporated Nb 
layer. The final geometry of the Nb electrodes was defined by electron beam 
lithography and an SFg RIE process. The junctions had a width of 20 |4m. In 
the plasma source employed for the cleaning, the gas is excited to a plasma 
in a tunable microwave cavity by a microwave source at 2.45 GHz. 

InP is known to be the less stable material and thus the surface of a semi- 
insulating InP wafer was exposed to the plasma. The effect on the surface 
morphology can be investigated by using an atomic force microscope ( AFM) . 
For an He/H (50 : 50) mixture, the surface topology measured by an AFM 
reveals the formation of droplets with a size of the order of 50 nm (Fig. 
4.10a). This result shows the same effect as reported in studies on InP surfaces 
exposed to a pure H plasma [116, 118]. Here, the exposure of an InP surface to 
hydrogen led to the formation of volatile PH3 and to an indium-rich surface 
covered with In droplets. As can be seen in Fig. 4.10b, the formation of 
indium droplets is prevented if pure He gas is used for the plasma source. 
The InP surface remains smooth after the He plasma treatment [119]. 

The effect of the plasma cleaning on the interface transparency of the 
Nb-InGaAs/InP interface can be investigated by measuring the differential 
conductance as a function of the voltage drop. In Fig. 4.11 i?Nd//dP is 
plotted as a function of the voltage drop. The temperature was increased 
from 1.0 K to 8.0 K in steps of 1.0 K. As can be seen here, the normal-state 
resistance Rn has a value of about 21 Q at 1 K, which would correspond 
to a value of 4 Q for a 100 |4m wide junction. Obviously, the normal-state 
resistance is noticeably smaller than for moderately Ar’^'-cleaned samples. 
As shown in Fig. 4.11 (inset), the He/H treatment leads to a much larger 
normal-state resistance of about 115 Q, which is attributed to the surface 
deterioration (Fig. 4.10a). 

Below 7 K, the curves plotted in Fig 4.11 show a minimum at zero bias 
voltage and maxima at about ±Aq. The critical temperature of this Nb layer 
was 8 K, which results in a superconducting gap of 1.2 meV. A comparison 
with the measurements reveals that the maxima are found at voltages below 
Aq{T). This can be explained by the proximity effect at the interface. The 
increase of dl / dV at zero bias for increasing temperature can be attributed 
to the broadening of the distribution function. For temperatures near Tc the 
lowering of Aq with temperature leads to a further damping of the structure. 

From the normalized differential conductance RT^dl/dU at zero bias, a Z 
parameter of 0.85 was deduced by following the BTK model. It is remarkable 
that this Z parameter is larger than, for example, the Z parameter of the 3 h 
Ar ■'■-sputtered sample, although the interface resistance per micron is consid- 
erably smaller, 400 Q |4m“^ for the He-cleaned sample in contrast to 1300 Q 
|lm“'^ for the Ar^'^ sputter-cleaned sample. This result can be interpreted as 
showing that for the He-cleaned surface, the area where transmission through 
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Fig. 4.10. (a) Atomic force microscope picture of an InP surface treated in an 
He/H electron cyclotron resonance plasma. The effect of an He plasma is shown in 
(b) [119] 



the interface takes place is larger than for the Ar+-treated sample. The barrier 
at the interface seems to be higher at the interface for the He-cleaned sample, 
since the Z parameter is larger. Finally, it can be stated that the He plasma 
cleaning is a very promising method, since here a clean and homogeneous 
interface is prepared without deteriorating the semiconductor surface. 



4.2 Point Contact Spectroscopy 

Point contacts are a versatile tool for studying transport through a normal- 
conductor/superconductor interface [99]. When a point contact is employed. 
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Voltage (mV) 

Fig. 4.11. Normalized differential conductance as a function of bias voltage at 
various temperatures for a sample cleaned in an He plasma. The inset shows the 
differential conductance for an He/H plasma-cleaned sample [119] 



carriers are injected or collected via a very small opening. In conjunction with 
a superconductor/two-dimensional-electron-gas interface, it is even possible 
to directly separate normal-reflected electrons and Andreev-reflected holes. 

An experiment where a superconductor/two-dimensional-electron-gas in- 
terface was analyzed by a point contact was reported by Jakob et al. [120]. 
A scanning electron beam micrograph of the structure used for point contact 
spectroscopy is shown in Fig. 4.12a. Again, the two-dimensional electron gas 
is located in a strained InGaAs/InP heterostructure. The point contact ge- 
ometry was defined by reactive ion etching. Separated by a distance of L = 1 
pin, an Nb electrode is attached to the sidewalls of the semiconductor mesa. 
As depicted in Fig. 4.12, a current is driven through the point contact to the 
superconducting electrode. In the case of normal reflection, the electron is 
specularly reflected and thus does not reach the point contact opening. On 
the other hand, if Andreev reflection occurs, so that an electron is retrore- 
flected as a hole, the hole trajectory follows the path of the incident electron 
and passes through the point contact opening. 

The effect of the retroreflected holes passing through the point contact 
can be seen in Fig. 4.13, where the differential conductance Gnu of the point 
contact is shown as a function of the excess energy. As can be seen, at ex- 
cess energy zero the conductance is enhanced by 20% compared with Gnu 
at excess energies exceeding the superconducting gap of Nb. The enhanced 
conductance can be attributed to the additional contribution of the retrore- 
flected holes. The modification of the point contact conductance Gpc if a 
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Fig. 4.12. (a) Scanning electron micrograph of a point contact with an Nb electrode 
in front. The separation between the point contact and the Nb electrode is 1 |im. (b) 
A normal-reflected electron can only pass through the point contact opening if it is 
almost incident almost perpendicular to the interface. Owing to the retroreflection 
property, an Andreev-reflected hole reaches the point contact opening for all angles 
of reflection [120] 



superconducting electrode is placed in front can be expressed as 



Gnu = Gpc 



1 -I- A{E) exp ( — 



2L 

lin,h{E) 



- B{E)r]e exp ( - 



2L 

^ball,e(.G) J 

(4.1) 



The Andreev reflection A{E) leads to an increase of the conductance, which 
is reduced only by inelastic scattering, characterized by the scattering length 
^in,h(A’), whereas the normal reflection B{E) leads to a decrease. Only a very 
small fraction rje « 0.02 of the electrons, which are incident almost perpen- 
dicular to the two-dimensional-electron-gas/superconductor interface, can in 
principle reach the point contact on their way back [120]. This amount is fur- 
ther decreased by elastic scattering (/baii.e); hence the normal-reflection con- 
tribution can be neglected in (4.1). As a consequence (4.1) can be rewritten 
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Fig. 4.13. Differential conductance Gnu as a function of the excess energy. The in- 
set shows the normalized differential conductance of the Nb-InGaAs/InP interface 
[ 120 ] 



as an expression where the Andreev reflection probability can be determined 
directly from GNb(^^): 
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/ Gf^hjE) 
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(4.2) 



A comparison of the experimentally determined Andreev reflection prob- 
ability with the BTK model [46] showed that this model was not able to 
explain the experimental results. In particular, the conductance peaks at 
±Z\q in the BTK model are not found in the experiment. A better agreement 
is obtained if the proximity effect is taken into account, as in the model by 
Golubov et al. [121] (see also Sect. 2.2.3). Besides the energy dependence of 
the Andreev reflection probability, Andreev reflection can also be studied as 
a function of temperature by using the measurement setup depicted in Fig. 
4.12. It could be shown that the Andreev reflection probability A{T) at zero 
excess energy decreases with increasing temperature, which is in contradic- 
tion to the constant value of A expected from the BTK model [122]. The 
decrease of A(T) with temperature was explained here by inelastic scattering 
at the superconductor/two-dimensional-electron-gas interface, which leads to 
an energetic smearing of the quasiparticle states. 

The angular dependence of Andreev reflection can be studied if two very 
small superconducting electrodes (« 200 nm) are employed, separated by 
a small angle [123]. In order to adjust the angle of the trajectories which 
couple the two superconducting electrodes, a reflecting barrier induced by a 
gate electrode was used. 




5. Superconductor/ 
Two-Dimensional-Electron- Gas / 
Superconductor Junctions: 
Voltage-Carrying State 



If two S/2DEG interfaces are brought into close proximity, fascinating effects 
can be observed. In this chapter, only the voltage-carrying state of S/2DEG/S 
structures will be discussed. In this regime, subharmonic gap structures at 
voltages equal to 2AqIti, n = 1, 2, 3, ..., and excess currents, for example, can 
be measured; these are due to multiple Andreev reflections. The different the- 
oretical models explaining these effects will be introduced first. Subsequently, 
the experimental results will be presented and compared with these models. 
Effects related to the superconductive state of S/2DEG/S structures will be 
covered in the next chapter. 



5.1 Subharmonic Gap Structures: Theoretical Models 

Two alternative models will be introduced here in order to explain the sub- 
harmonic gap structures as well as the excess currents found in the current- 
voltage {I-V) characteristics. The first model is based on the motion of Gaus- 
sian wave packets within the normal conductor [124]. The second model, 
which will also be employed to interpret the experimental data, deals with 
ballistic particles in a Boltzmann equation approach [125]. For the latter 
model, the inclusion of the proximity effect will also be discussed briefly [58]. 



5.1.1 Kiimmel Gunsenheimer Nicolsky Model 

Kiimmel, Gunsenheimer and Nicolsky (KGN) [124] developed a model, where 
the I-V characteristics of a SNS junction are calculated by considering the 
trajectories of single electrons or holes. No barrier at the interface is included 
{Z = 0). The particles, described as Gaussian wave packets, are accelerated or 
decelerated by the electric field E in the normal conductor depending on their 
initial momentum. Between the Andreev reflection events at the interface the 
particles climb or descend the Andreev ladder as depicted in Fig. 5.1a. As long 
as the energy of the particles at the SN interface is within the superconducting 
gap, ideal Andreev reflection occurs. If the applied voltage is decreased, more 
and more Andreev reflections are necessary to gain or lose enough energy 
before the particle can enter the superconductor as a quasiparticle above 
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Zip or below —Zip. The increase of the particle contribution to the electrical 
current by multiple Andreev reflections is balanced by inelastic scattering 
processes. The latter limit the number of Andreev reflections. 



a) b) 




Fig. 5.1. (a) Electrons e and holes h accelerated by an electric field E = —exV/L. 
As long as the particle energy is within the gap ±Ap, Andreev reflections occur 
at the SN interfaces. Ek is the initial energy of the particles, (b) Energy range 
of electrons which can undergo one (A^) or two (Aj) Andreev reflections. The 
electrons start at a; = 0 



The total current through the SNS junction results from a competition 
between two opposite currents carried by electron and hole wave packets; 
the electrons and holes moving with the held make a positive contribution, 
while the corresponding particles moving against the held make a negative 
contribution. Calculating the time and spatial average of the Gaussian wave 
packets and using the symmetry properties of the Andreev reflection coeffi- 
cients leads Anally to the following expression for the Andreev current density 
[124]: 

Jar = -ex — 7^ ^ ftr(£’fc){/o(Afc)fce - [1 - 

k 

oo 

X ^ ^-(n+l/2)L/W. [Al{Ek) - A-{Ek)] . (5.1) 

n—1 

Here, fee and feh are the initial electron and hole wave vectors, respectively. 
17 n is the volume of the normal conductor and P^{Ek) is the probability to 
And the particle in the N layer. L is the distance between the superconducting 
electrodes, while kn is the inelastic mean free path. In the first sum all possible 
initial k vectors are considered. The second sum adds up the steps in the 
Andreev ladder. The Andreev coefficients A^{Ek) represent the probability 
that at least n Andreev reflections (see Fig. 5.1) are possible: 
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A^{E) = 



Here, A{E) is approximated by A{E) = 1 for \E\ < A and 0 elsewhere [124]. 
The (+) and (— ) signs represent particles moving with and against the field, 
respectively. A direct consequence of the form of A^ is the appearance of 
subharmonic gap structures at voltages 

Vn = 2Aolen, n= 1,2,3,... (5.3) 

At 14, a transition from n — 1 to n Andreev reflection processes takes place. 

In (5.1), only the additional current due to Andreev reflection is consid- 
ered. For the total current density, the normal-state current density has to be 
added. The excess current Axe is defined as the Andreev current accumulated 
at voltages far above 2Z\q. 

In order to calculate the current density Jarj information about the den- 
sity of Andreev states g{E) is needed. Information about g{E) can in principle 
be gained from the sum over k if one knows the dispersion of the Andreev 
states Eik'Px) as a function of the x component of the Fermi wave vector 
[124]. For an ideal SNS junction, an expression for the energy levels E(k-px) 
was given by Kulik [126], for example. Pereira et al. [127, 128] showed that 
for a fit of an experimental I-V characteristic, it can be sufficient to assume 
a simplified analytical Andreev density of states. 

One prediction of the KGN model is that for a large hn/L ratio a negative 
differential conductance (d/ /dV < 0) should be observable. The current in- 
crease with decreasing bias voltage is due to the increasing number of Andreev 
reflections. The only limiting process is the inelastic scattering. For a current- 
biased junction, the negative differential conductance will lead to a hysteresis 
in the I-V characteristics. In addition, the steep rise of the current density 
at low bias voltages, the characteristic foot of weak-link structures, can also 
be explained within the KGN model. Recently, a more rigorous treatment 
of inelastic scattering processes was provided by Gunsenheimer and Zaikin 
[129]. 
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5.1.2 Inclusion of Normal Reflection at the SN Interface 

In the Kiimmel-Gunsenheimer-Nicolsky model, no barrier at the interface 
was assumed, and therefore ideal Andreev reflection occurs for particle en- 
ergies within the superconductive gap. If a (5-shaped barrier is introduced 
at the interfaces, the particle wave packets are partially normal-reflected at 
the interface even for \E\ < Aq, as shown in Fig. 5.2. We shall now in- 
troduce a simple extension of the KGN model, where in (5.1) the sum over 
exp[—{n+l/2)L/lin][A:f{Ek) — A~{Ek)] will be replaced by a term which also 
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includes normal reflection processes [130]. Interference effects are neglected 
in this model. 

The probability to And an electron or a hole at x = 0 with energy E 
is deflned as Cf or , respectively. The (+) sign represents motion with 
the fleld, the (— ) sign motion against the fleld. As shown in Fig. 5.2, a hole 



Z6 Z5 




Fig. 5.2. Trajectories of Andreev and normal-reflected particles in an SNS struc- 
ture. and are the electron and hole probabilities, respectively, at x = 0; (-I-) 
denotes motion with the field, while (— ) denotes motion against the electric field 



with probability C^{E + eV), with energy E + eV at x = 0, which left the 
interface at x = T/2, moving with the fleld, can originate either from an 
Andreev-reflected electron with probability C^{E), of energy E, or from a 
normal-reflected hole with probability C^{E + eV), of energy E + eV, at 

X = 0: 

C+{E + eV) = e~^/^^''[A{E + eV/2)C+{E) 

+ B{E + eV/2)C^{E + eV)]. (5.4) 

A{E + eVj2) and B{E + eV/2) are the BTK coefficients for Andreev and 
normal reflection. A corresponding process can be found for an electron with 
probability C~{E), at energy E moving against the fleld: 

C-{E) = e-^/^^-'[B{E + eV/2)C+{E) 

+ A{E + eV/2)C^{E + eV)] . (5.5) 

The factor exp(— L//in) accounts for losses by inelastic scattering in the nor- 
mal conductor. The equations given above show that (7+ and C~ at energy 
E can be deduced from the values of and at the energy E + eV one 
step higher: 

e^/h.c+{E + eV/2) - B{E + eVl2)C^{E + eV) ^ 
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C-{E) = 



B{E + eV/2) 



C+{E + eV) 



A{E + eV/2) 

-L,u^ + eV/2) - B^{E + 6 ^/ 2 )] C^{E + eV) 

A{E + eV/2) 

(5.6b) 



+e 



One step further down, C^{E — eV) and C^{E — eV) can be written as 
functions of C^{E) and C~{E). 

Since electrons and holes have the same reflection and transmission prob- 
abilities at the SN interface, two series C^{E + neV) and C~{E + neV) can 
be defined, which describe the motion of a particle with and against the elec- 
tric held, respectively. Suppose that the initial particle has been an electron 
with energy Ek, then the even terms C^{Ek + 2neV) correspond to electrons, 
while the odd terms C^[Ek + {2n + l)eV)] belong to holes. 

The boundary conditions are given by the fact that for particle energies 
Ek + neV far beyond Aq, the Andreev reflection process can be neglected. 
Only normal reflection has to be considered. For Ek + neV ^ Aq, this means 



C~{Ek + neV) = B{Ek + neV + eV/2)C+{Ek + neV) , (5.7) 



while for the lower side {Ek + neV <§C — ^o)j one obtains 

C+{Ek + neV) = B{Ek + neV + eV/2)C-{Ek + neV) . (5.8) 



The total solution is given by a linear combination of two solutions 
satisfying the boundary conditions above Aq, and of a solution 
(Cq,Cq) with nonzero values only for E < Ek and with a step in Cq or 
Cq equal to one at Ek for an initial particle moving with or against the 
held, respectively. The linear combination is chosen in such a manner that 
the boundary condition for energies far beyond —Aq is satisfied: 






C+{E_) - B{E_ - eVl2)CQ{E_) 
^ ^ ^ Ct{E-) - B{E_ - eVl2)C^{E_) 



Ct(E) , 



(5.9) 



where = Ek + neV <C —Aq. The step at Ek is normalized to 1. The 
Andreev current is obtained by substituting X)™ ~ 

^n{Ek)\ in (5.1) by 



EE [C+{Ek + neV)A{Ek + neV + eV/2) 

p=w,a n 

- C~{Ek + neV)A{Ek + neV -eV/2)] . 

P = w,a denotes generation of a particle at Ek moving with (w) or against 
(a) the held. 



5.1.3 Octavio Tinkham Blonder Klapwijk Model 

The Octavio-Tinkham-Blonder-Klapwijk (OTBK) model [125] is based on 
a Boltzmann equation approach. As in the modified KGN model, normal 
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scattering at the SN interface is included by introducing a J-shaped barrier, 
whose strength is expressed by the Z parameter. In contrast to the previously 
discussed model, inelastic scattering processes are not considered here. The 
OTBK model can be regarded as a generalization of a preceding model with 
no barriers at the SN interface developed by Klapwijk et al. [131]. 

In order to calculate the current as a function of the applied voltage V 
between the superconducting electrodes, the electrons are divided into two 
subpopulations f^{E, x), f^{E, x) representing the two directions of motion. 
The nonequilibrium distributions at the left-hand SN interface (x = 0) and 
at the right-hand NS interface (x = L) are calculated from the boundary 
conditions 

f^{E,0) = A{E)[l-f^{-E,0)] + B{E)f^{E,0) + T{E)fo{E), 

(5.2010a) 

f^{E,L) = A{E)[l-f^i-E,L)] + BiE)f^iE,L) + T{E)foiE). 

(5.2010b) 

Here, A{E), B{E) and T{E) are the coefficients for Andreev reflection, nor- 
mal reflection and transmission through the interface, respectively. For the 
transmission from the superconductor, a Fermi distribution function fo{E) 
is assumed. Since no inelastic scattering within the normal conductor is as- 
sumed, the reference point for all distribution functions can be shifted to 
X = 0; for example an electron starting at x = L with energy E will arrive 
at X = 0 with energy E — eV : 

f^{E,L) = f^{E-eV,0) = U{E-eV) . (5.2011) 

The net current for a two-dimensional structure as a function of the ap- 
plied voltage is determined by the difference between the distribution func- 
tions f^{E): 

I = r [f^{E] - f^{E)] . (5.2012) 

Here, W is the width of the contact. The distribution functions themselves 
are calculated self-consistently [125, 132]. As an example, f^{E) and f^{E) 
are plotted in Fig. 5.3 for Z = 0.5, T/Tc = 0.6 and V = QAAoje. The 
appearance of subharmonic gap structures in the form of pronounced steps 
at Vn = eAol{en) are clearly seen in Fig. 5.4 for a structure with Z = \. 

The OTBK model for a two-dimensional structure can be extended by 
considering an angle-dependent scattering probability. Here, it must be taken 
into account that only the wave vector component normal to the interface 
is altered by the scattering at the barrier. The component parallel to the 
interface is preserved. As can be seen in (2.8), the Z parameter is defined 
with respect to the Fermi wave vector feps. In order to account for the fact 
that only the wave vector component normal to the interface is affected, the 
Z parameter has to be replaced by an angle-dependent factor Z{fE) = Zj cos 9 
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Fig. 5.3. Distribution functions and for Z = 0.5 and a Fermi distribution 
function at T/Tc = 0.6. The applied voltage is F = QAA^/e 



[54]. Here, 9 is the angle of incidence at the interface with reference to the 
axis normal to the interface. For particles with nonnormal incidence on the 
interface, an effectively larger Z parameter is taken into account now. The 
expression for the current including Z{9) is given by 



The I-V characteristics calculated from this expression with a given Z pa- 
rameter will correspond qualitatively to the I-V characteristics calculated 
from (5.2012) with a a slightly larger Z parameter. 

The OTBK model can be generalized by considering superconductors with 
different energy gaps on the two sides (SNS') [131, 133, 134] or by assuming 
different barriers {Z parameters) at the SN interfaces [135, 136]. Ineleastic 
scattering processes were included in the OTBK model by Heslinga et al. 



5.1.4 Inclusion of the Proximity Effect 

In Sect. 2.2.3 the effect of the proximity effect on the Andreev reflection pro- 
cess was discussed. Here, the semiconductor (Sm) is not in direct contact 
with the superconductor (S) but rather with a normal conductor (N), which 



I 





{f^[E,Z{9)] - f^[E,Z{9)]} . (5.2013) 



[137]. 
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Fig. 5.4. Normalized differential resistance (dV/d/)_RN as a function of the bias 
voltage eV/Z\o calculated from the OTBK model for Z — 0 and T = 0 {solid line). 
The dashed and dotted lines show (dV/d/)J?N if the proximity effect is included, 
for 7 b = 1 and 7 b = 5, respectively. The arrows indicate the first subharmonic gap 
structures of the 14 = 2At^/{en), Vn = {As — At^)/{en), V„ = (Zis + 4\N)/(en), and 
14 = 2As/{en) series (from left to right). The data was taken from [58] (© (1996) 
by the American Physical Society) 



is modified owing to the proximity effect. In order to suppress the proximity 
effect between the semiconductor and the N layer, the two layers are con- 
nected via a small constriction c. As a result, the density of states (DOS) of 
the N layer at the interface with the Sm layer possesses two peaks, one at 
the induced gap Z\n and one at 2ls = (Fig- 2.7), in contrast to the bulk 
DOS of a superconductor, which has only a single singularity at Aq [-53, 67]. 
Consequently, the Andreev reflection coefficient A{E) has also a two peak 
structure, as was shown in Fig. 2.8. 

Of course, the additional structures in the Andreev reflection coefficient 
A{E) have a direct consequence for the subharmonic gap structures. A more 
complex structure can be expected in the differential resistance dV/dl, which 
may be due either to the gap of the superconductor Z\s, to the induced gap 
Z\n or to a combination of the two, Z\s i Z\n- A generalization of the OTBK 
model in this respect was provided by Aminov et al. [-58]. Here, an SN/c/NS 
structure, instead of the S/c/S structure in the OTBK model, is considered. 
The quasiparticle trajectories which lead to the structures in the differential 
resistance are shown schematically in Fig. 5.5. 

The (dK/d/)/i?N vs. eVjAs curves calculated by applying the extended 
model of Aminov et al. [58] are plotted in Fig. 5.4 for a discontinuity pa- 
rameter between the S and N layers of 7 b = 1 and 7 b = 5. The suppression 
parameter 7m was set to 0.01, while the Z parameter between the N layer 
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Fig. 5.5. Schematic illustration of the quasiparticle trajectories in the constriction 
region suffering multiple Andreev reflection at T = 0 K, at bias voltages of (a) 
eF = {As - An)/3, (b) eF = {As + An)/3 (c) eF = 2Zls/3 and (d) eF = 2Z1n/3 
[58] (© (1996) by the American Physical Society) 



and the constriction was equal to one. Compared with the curve calculated 
from the OTBK model, also shown in Fig. 5.4, the curves calculated from 
the extended model show many additional features. These structures can be 
explained by referring to the diagram shown in Fig. 5.5. At F„ = 2Z\s/(en), 
the structures already known from the OTBK model are found (Fig. 5.5c), as 
indicated for n = 1 (arrow) at eV/As = 2 for 7 b = 5. This structure is very 
weak, since only a small peak at As remains in the density of states at the 
N/c boundary. The much more pronounced structures at eV/Ag « 1.6 and 
0.75 belong to the series Vn = (As + AN)/(en) and Vn = (As — AN)/(en) 
with n = 1, respectively (Figs. 5.5b, a). Note that, in contrast to the diagram, 
for the calculated curve An is smaller than As/2 for 7 b = 5 (Fig. 2.7). The 
most pronounced dip in the differential resistance is found for cF/An = 2 
at eF/As « 0.5. This structure corresponds to the enhanced differential con- 
ductance at cV/Aq = 2 in the OTBK model, since the induced DOS at the 
N/c boundary has its maximum at An. 
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5.2 Measurement of the Voltage-Carrying State 

5.2.1 Andreev Current in S/2DEG/S Structures 

The fact that the Andreev reflection process can indeed lead to an additional 
current contribution is visualized in Fig. 5.6. Here, the Andreev current /ar 
of a 450 nm long and 6 |im wide Nb-InGaAs/InP-Nb junction, prepared as 
described in Sect. 3.2.2, is plotted as a function of the voltage drop. The sheet 
electron concentration of the heterostructures used here was 7.05 x 10^^ cm“^ 
at 0.3 K, with a mobility of 270 000 cm^/Vs. An elastic mean free path lei of 
3.7 gm is deduced from these values. The Nb layer had a critical temperature 
of Tc = 8.9 K, related to a superconducting gap of Aq = 1.35 meV. 

/ar is calculated from the total current I by subtracting the contribution 
of the normal current: /ar = / — V/R^, see Fig. 5.6 (inset). /?n has a value 
of 50 Q. The measured /ar-E characteristics reveal a negative differential 
conductance at small voltages V < 0.1 mV. Since the measurement was 
performed in a current bias mode, a shunt resistor of 6.8 Q was connected in 
parallel with the junction in order to gain access to the negative-differential- 
conductance range. Details of the measurement setup are given in Sect. 6.2.3. 




V(mV) 



Fig. 5.6. Measured Andreev current as a function of V {solid line). Also plotted 
are the calculated Iar~V curves obtained from the KGN model for hn/L = 10 
{dashed line) and /„// = 2.5 {dash-dotted line) and from the modified KGN model 
for lin/L = 10, Z = 0.5 {dash-dot-dotted line). The inset shows the measured total 
current {solid line) and the ohmic characteristics {dashed line) extrapolated from 
V > 2Aq. The intersection with the current axis is the excess current /exc 
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It should be remarked that in this S/2DEG/S junction, a finite excess 
current of 10 |J,A is found. The excess current is taken here as the Andreev 
current Jar obtained for voltages far above 2Z\o- The fact that an excess 
current, in contrast to a deficit current [135], is observed means that in this 
structure the current gain due to the Andreev reflection is larger than the 
current loss due to the normal reflection process. 

The occurrence of a negative differential conductance could lead to the 
conclusion that this is an effect of current enhancement by multiple An- 
dreev reflections as predicted by the KGN model [124]. In order to compare 
the measurement with the results of the KGN model, two calculated curves 
based on this model are also plotted in Fig. 5.6. For a large inelastic mean 
free path {kn/L = 10) the theoretical curve shows the predicted negative 
differential conductance. However, the overall shape of this curve does not 
fit very well to the experimental results. First, the negative differential con- 
ductance is observed for 0.1 mV < V < 2.6 mV, whereas in the experiment 
a positive differential conductance is measured in this range. Second, the so- 
called “foot” appears for V <0.1 mV in the calculated curve. In contrast, 
no foot structure but, rather, a negative slope is found in the experimental 
curve. A further analysis shows that it is impossible to obtain theoretically 
an increase of Jar with increasing V in the voltage range where this was ob- 
served in the experiments and to maintain a negative differential conductance 
at low voltages. This is demonstrated by the curve for kn/L = 2.5 in Fig. 5.6. 
It can thus be concluded that the experimentally observed negative differ- 
ential conductance cannot be explained within the KGN model. In addition, 
as shown by Uhlisch et al. [138], the inelastic mean free path for this kind 
of heterostructure is larger than kn = 2.5 L, which supports the assumption 
that the experimentally obtained Andreev current must be restricted by an 
additional process. 

As can be seen in Fig. 5.6, the overall shape of the experimental Iar(I^) 
curve can be better reproduced by applying the modified KGN model. Here, 
(5-shaped barriers at the SN interfaces are assumed. A good fit to the exper- 
imental data is obtained by assuming a Z parameter of 0.5 with kn/L = 10. 
Nevertheless, the additional normal-reflection contribution also prevents a 
negative differential conductance, so that this experimentally observed fea- 
ture cannot be explained within this framework either. In Sect. 6.2.3 the 
discussion on the negative differential conductance will be continued by in- 
cluding the Josephson effect. 

5.2.2 Subharmonic Gap Structures 

Subharmonic gap structures can best be resolved by measuring the differential 
resistance dV/dl as a function of the voltage drop. In Nb-InGaAs/InP-Nb 
junctions, pronounced subharmonic gap structures are found in the dV/dI(V) 
measurements as shown in Fig. 5.7 for a 450 nm long, and 6 pm wide Nb- 
InGaAs/InP-Nb junction. 
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Fig. 5.7. Normalized differential resistance as a function of bV/Aq {solid line) at 
T = 1.5 K. The theoretical curves obtained from the OTBK model were calculated 
for Z = 0.55 (v) and Z = 0.6 (A) and in addition for Z = 0.55 (•), taking the 
angle-dependent Z parameter into account 



From the discussion in the preceding section it was concluded that the 
barriers at the interface mainly determine the magnitude of the Andreev 
current, while the influence of the inelastic scattering is less pronounced. 
This is the reason why, in the following discussion the OTBK model, which 
neglects inelastic scattering, is referred to. 

In Fig. 5.7, curves calculated using (5.2012) for Z = 0.55 and Z = 0.6 
are added. The features of the calculated subharmonic gap structures at 
eV = 2Aq, Aq and 2Z\o/3 agree reasonably well with the measured struc- 
tures. For Z = 0.6 the local maximum near Aq fits the measured maximum. 
However, the dV/dl modulations in the calculated curves are generally more 
pronounced than the experimentally obtained variations of the differential re- 
sistance. If the dV/dI~V characteristics are calculated from (5.2013), which 
considers the angle dependence of the Andreev reflection, the values of the 
differential resistance for eV smaller than « 1.6zio are larger than the corre- 
sponding values calculated from (5.2012). This is demonstrated in Fig. 5.7, 
where the characteristic for Z = 0.55 calculated from the extended OTBK 
model is plotted. Again, the dV/dl modulation is more pronounced than in 
the experiment. The less developed minima in the experimental curve might 
be due to scattering processes within the normal conductor or due to diffusive 
scattering at the interface. The OTBK model considers only the ballistic case 
with ideal Andreev and normal reflection at the interface. From studies of 
single S/2DEG junctions, it is known that the interface has a certain rough- 
ness due to the Ar+ sputter cleaning. This roughness can lead to diffusive 
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scattering at the interface. Furthermore, the scattering of carriers at residual 
donors in the channel can also lead to damping of the dV/dl modulations. 

An additional minimum at eV = 1.2 Aq, together with a broad maximum 
at 1.55Z\o, is found in the measured curve, which is not predicted by the 
OTBK model. A possible origin of the minimum could be a difference in 
the gaps in the superconductive electrodes [131, 133, 134] but since both Nb 
electrodes were prepared simultaneously, this effect could be excluded for our 
samples. Probably these additional structures are due to the influence of the 
proximity effect [58], where the maximum at 1.55 Z\q can be attributed to 
2 z1n, where An is the proximity-induced energy gap. 




Fig. 5.8. Differential resistance of a planar Nb-InAs junction as a function of 
voltage, with the temperature varied as a parameter [139] 



The occurrence of an excess current Jexc can be deduced directly from 
the fact that the major part of the normalized differential resistance for 
eV/Ao < 2 is smaller than 1. From the normalized excess current 
the Z-parameter can be estimated following the approach introduced by 
Flensberg et al. [132]. For Jexc = 10 flA and J?n = 50 Q, a Z parameter 
of 0.9 is estimated. This factor is larger than the Z parameter used for the 
modified KGN model. In the OTBK model Jexc is mainly determined from the 
minima of the oscillations, which are deeper than the experimentally observed 
minima. If the angle-dependent Andreev reflection is taken into account, a 
slightly smaller Z parameter results for a given excess current. 
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The differential resistance of a planar Nb-InAs junction is plotted in 
Fig. 5.8 as a function of the voltage drop, with the temperature varied as 
a parameter. It can be clearly seen that the subharmonic gap structures shift 
towards smaller voltages if the temperature is increased. The temperature 
dependence follows Ao{T) of the BCS theory. 

In Fig. 5.9 the subharmonic gap structures are plotted for samples with 
different geometries (samples A-F). The length of the junctions varied from 
300 nm up to 800 nm. Even for large separations of the Nb electrodes, an 
excess current is observed indicated by the lowering of dl /dV below 2 Z\q = 
2.7 meV. For larger electrode spacings the subharmonic gap structures are less 
pronounced. This can be attributed to an increasing scattering probability of 
the particles when they move through the 2DEG. The pronounced structures 
at low voltages observed for samples A and C are not due to multiple Andreev 
reffections but rather due to the Josephson effect, as explained in Sect. 6.2.4. 
The absence of these structures for the 800 nm and 600 nm long samples 
is due to a small magnetic field (< 1 mT), which suppresses the Josephson 
effect. 




Fig. 5.9. Normalized differential resistance for different junction geometries as a 
function of the voltage drop. The corresponding lengths and widths are given in 
the graph. The curves are shifted upwards in steps of 0.1 



5.2.3 Subharmonic Gap Structures and the Proximity Effect 

In some junctions the proximity effect has to be taken into account in order 
to interpret the subharmonic gap structures. Here, the SN/c/NS model in- 
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troduced in Sect. 5.1.4 can be applied instead of the OTBK model. A typical 
example of a sample showing evidence for a proximity layer is an Nb/p- 
InP/Nb junction, as depicted in Fig. 3.5b [59]. The proximity layer (N) is 
formed in the InAs inversion layer directly underneath the superconducting 
Nb electrodes, as illustrated in Fig. 5.10 (inset). Because of the strong cou- 
pling to the superconductor due to the absence of a barrier, a gap in the 
density of states is induced in the inversion layer (N) below the Nb electrode. 
The two N layers are connected via a short bridge, constituted from the two- 
dimensional electron gas in the surface inversion layer of the p-InAs. This 
bridge serves as the constriction c in the SN/c/NS model. 

As can be seen in Fig. 5.10, the measured subharmonic gap structures 
cannot be explained consistently by the OTBK model. In order to obtain a 
reasonable agreement with the OTBK model at low temperatures (1.8 K), 
a lower gap As of 1.2 meV compared with the value determined from the 
critical temperature has to be assumed. However, a look at the measured 
curve at a higher temperature of 6.6 K shows that the subharmonic gap 
structures calculated from the OTBK model appear at bias voltages that are 
much too low. The reason for this discrepancy is the gap value calculated 
from the BCS theory, which is too small at 6.6 K. 

A much better agreement between experiment and theory is obtained if a 
proximity layer with an induced gap An is assumed. In this case, it is possible 
to take a gap value of 1.45 meV in the superconductor (S), deduced from the 
measured Tc of the Nb layer. According to the SN/c/NS model, a reduced 
gap An is found at the boundary with the bridge c formed by the two- 
dimensional electron gas at the surface of the InAs. In the calculated curves 
the subharmonic gap structures at Vn = 2 An / (en) due to the induced gap An 
are much more pronounced than the remaining structures at Vn = 2As/(en) 
or Vn = (As ± AN)/(en) [59]. This is probably due to the fact that for the 
parameter set taken here, the peak at E = As in the induced density of 
states is only weak. According to theory, the induced gap An approaches the 
gap of the superconductor layer As as T ^ Tc [58]. This is the reason why 
the subharmonic gap structures in the SN/c/NS model are found at larger 
voltages than in the OTBK model and thus agree better with the experiment. 



5.3 Resonant Strnctnres and Interference Effects 

In some junctions the differential resistance reveals features which cannot 
explained by the subharmonic gap structures discussed above. A typical ex- 
ample is shown in Fig. 5.11a, where the differential resistance AV/AI of a 
p-InAs step edge junction (Fig. 3.5a) is plotted [140]. For voltages below the 
first subharmonic gap structure at 2 An/ 6 (upward arrow), distinct subgap 
peaks are observed, which cannot be explained either by the BTK or by the 
SN/c/NS model. These subgap peaks are much more pronounced than the 
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V (mV) 

Fig. 5.10. Comparison of the measured differential resistances {solid lines) of a 
Nb/p-InAs/Nb junction at four temperatures with the SN/c/NS model {dotted 
lines) and the OTBK model {dashed lines). The parameters for the SN/c/NS model 
are As{T = 0) = 1.45 meV, 7 ?n = 1-55 12, 7b = 0.65 and Z = 0.45. For the OTBK 
model, the following parameters were used: As{T = 0) = 1.20 meV, Rjq = 1.55 12 
and Z = 0.45. The curves for 1.8 K and 8.0 K are shifted by 0.7 and 0.5 12, 
respectively. The inset shows the realization of the SN/c/NS junction in a planar 
Nb/p-InAs/Nb junction [50] (© (1997) by the American Physical Society) 



subharmonic gap structures. If the temperature is increased, the peak posi- 
tions shift towards lower voltages. In Fig. 5.11b the position of the dominant 
peak (<0>) is plotted versus temperature. It can be seen that the peak position 
does not scale with the induced gap 2A-^/e but, rather, maintains a con- 
stant distance from 2A-^/e. Comparable resonances have also been observed 
in other types of junctions, e.g. junctions based on InAs/AlSb [28, 36]. 

A qualitative explanation of the sharp resonances was suggested by Bas- 
tian et al. [28]. According to their interpretation, the peak is caused by quasi- 
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Fig. 5.11. (a) Differential resistance AV/dl of a 2.5 |lm wide p-InAs step junction 
at different temperatures as a function of bias voltage. For clarity, the curves are 
shifted in steps of 5 12. The arrows at the bottom indicate the induced energy gap 
2A^/e. The arrows at the top point to a giant subgap peak, (b) Temperature 
dependence of 2A^/e (•) and the giant subgap peak (<C>) The dotted line is the 
calculated value of 2As/e (BCS theory). The dashed line is the calculated value 
of 2A^je, according to (2.32), with A^jA^ = 0.8 at T = 0 and an interface 
transparency parameter ye = 0.25 (SN/c/NS model) [140] 



particle interference. In Fig. 5.2 the electron and hole trajectories in an SNS 
junction are shown. Owing to the finite barrier at the interface {Z >0), there 
is a certain probability for normal reflection. By this normal reflection process, 
a quasiparticle can return to the same energy and position where it started 
from. For those voltage drops which satisfy the condition for constructive 
interference, a standing wave is formed. Since the particle is trapped in the 
cavity, the differential resistance is increased, as observed in the experiment. 
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The interpretation given above is supported by measurements on junc- 
tions with different separations of the superconducting electrodes. With an 
increasing width of the cavity, the number of spikes is increased as expected 
[28] . The shift of the peaks towards smaller voltages with increasing temper- 
ature can be explained by the fact that the effective width of the cavity is 
increased. The cavity is bounded by two induced proximity layers (for details 
see Sect. 5.2.3). With increasing temperature the width of these proximity 
layers decreases, leading to an effective widening of the cavity and thus to a 
shift of the spikes to lower voltages. 



6. Josephson Current in Superconductor/ 
Two-Dimensional-Electron-Gas Junctions 



During the discussion of the subharmonic gap structures, the phase of the 
superconducting electrodes was neglected. We shall see in this chapter that 
phase coherence is crucial for the occurrence of a Josephson supercurrent. In 
the first part of this chapter, the theoretical concepts will be discussed, start- 
ing with the simplest model of an ideal one-dimensional junction, without any 
barriers at the SN interface and with no Fermi velocity mismatch. Later on, 
more elaborate models will be introduced, which describe the experimental 
situation more appropriately. 

In the second part of the chapter, measurements of the superconductive 
state are discussed. After analyzing the current-voltage characteristics of var- 
ious types of Josephson junctions and comparing them with the theoretical 
models, the interference effect due to an applied magnetic field is discussed. 
The light sensitivity of the sheet resistance in a semiconductor heterostruc- 
ture can be employed to control the supercurrent in a Josephson junction. 
This will be demonstrated in the last section. 



6.1 Theoretical Description of the Supercurrent 
in an S/2DEG/S Junction 

For the description of a Josephson supercurrent within the framework of 
Andreev reflection, it is essential to maintain the phase coherence between the 
two electrodes. In the clean-limit case < le\) it was shown that a significant 
portion of this supercurrent is carried by discrete energy levels, which are 
formed within the superconductive gap [126]. The second part contributing to 
the total supercurrent is due to the continuous energy spectrum at \E\ > Aq. 
For the ideal case, an expression for the discrete energy- level spectrum with 
no barriers at the interface and no Fermi velocity mismatch was derived 
for a one-dimensional structure by Kulik [126], Ishii [141] and by Bardeen 
and Johnson [142]. The next step towards a more realistic description of 
the experimental situation was provided by Bagwell [143], who discussed the 
effect of a single impurity within the normal conductor modeled by a J- 
shaped barrier. A further generalization was achieved by Tang et al. [144], 
where J-shaped barriers were assumed at the SN interface. In addition, a 
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Fermi velocity mismatch was also included here. It is sometimes interesting, 
especially for modeling structures with regard to potential applications, to 
introduce an arbitrary potential shape, e.g. tunnel barriers, in the normal 
conductor, as discussed in [13, 145, 146]. 

In order to describe the transport in a realistic device, the one-dimensional 
model has to be extended to two or three dimensions. A detailed model for an 
Nb/n-type-InAs/Nb junction was provided by Schiissler and Kiimmel [76]. 
For a two-dimensional electron gas sandwiched between two superconducting 
electrodes, a model was developed for the ideal case with no barriers at the 
S/2DEG interface by Kresin [147]. Later this model was extended by Chrestin 
et al. [148] by including J-shaped barriers at the S/2DEG interfaces. 

6.1.1 Supercurrent in an Ideal One-Dimensional SNS Junction 

In order to understand the basic transport mechanisms in an SNS junction, 
it is instructive to discuss an ideal one-dimensional junction. “Ideal” means 
that neither barriers at the interface nor differences in the Fermi velocity be- 
tween the superconductor and the normal conductor are considered [126]. In 
contrast to a single SN interface, now the phase difference </> of the pair poten- 
tial between the two superconducting electrodes must be taken into account 
(Fig. 6.1). Similarly to the BTK model [46], a step-like superconducting pair 
potential A{x) is assumed, which is zero in the normal conductor: 

( Aoe-'-^/^ , a:<0 

A{x) = < 0 , 0 < a: < L . (6.1) 

[ Aoe‘‘^/2 , x>L 




A A 

0 N 0 



h ^ 



Fig. 6.1. Spatial dependence of the pair potential A{x) for an ideal one-dimensional 
SNS junction. The arrows represent the group velocities of an electron and a hole. 
The wave vectors of the two particles have the same direction 



The energy spectrum of an ideal SNS junction, as a solution of the the Bogoli- 
ubov equation (2.1), can be separated into a continuous Andreev spectrum 
for \E\ > Aq and discrete Andreev states for \E\ < Aq. The latter are due to 
the fact that states within the gap of the superconductor cannot exist. Here, 
the wave function attenuates exponentially in the superconductor. In a sense, 
the situation is comparable to a finite-depth quantum well for electrons, with 
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discrete states within and extended states outside the well. However, in con- 
trast to a quantum well, the Andreev bound states in an SNS junction can 
carry a supercurrent. 

Following the ansatz of Kulik [126], the wave functions ip± for excita- 
tion moving to the right (-I-) and to the left (— ) for \E\ < Aq can be expressed 

by 



A+ ( ) e+‘'=""^ 

' uo 



X < 0 



V'+(x) = < 






C 4 



voe 



-i(p/2 



^ e+ifeex + ^ , 0 < X < L (6.2) 



] g+ifee(x-L) 

xo 



X > L 



and by 



V’-(x) = < 



A_ 



B_ 



C- 



uoe 
xo 



, X < 0 



^ + B_(^ ^ e-‘'=‘>“ , 0 < X <L. (6.3) 



Uo 



, X > L 



Here, kg and fch are the wave vectors given by (2.12a) and (2.12b) for an elec- 
tron (e) and a hole (h), respectively, in the normal conductor. Note that the 
group velocity of the hole has the reverse direction to its phase velocity. The 
latter is oriented along fch- fee and fch are the wave vectors of an electron-like 
(e) and a hole-like (h) particle in the superconducting electrodes, as given 
in (2.13a) and (2.13b), respectively. For the ideal case, as considered here, 
the Fermi wave vectors in the normal conductor and the superconductor are 
identical; fcpN = ^fs = ^f- Since |if| < Aq, ke and fch possess an imaginary 
component, which is responsible for the aforementioned exponential attenu- 
ation within the superconductor. Continuity at x = 0 and x = L and use of 
the definition of uq and Vq, as given by (2.4a) and (2.4b), result in 

4 = ~ -E^ = _ ( 5 , 4 ) 

Cq ^0 '^0 

Here, +\4> corresponds to the right-moving particles and —i<j) to the corre- 
sponding left-moving particles. Since |£’| < Aq, the absolute value of Uq/xq 
is 1, thus we can define Uq/vq = cos2y -|- isin 27 = exp(i2y). Inserting this 
expression in (6.4), we arrive at 

I — g-i27g±i0gi(fce-feh)i 



(6.5) 
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For \E\ /i, (fee — fch) can be approximated by kpE/fi = E/{Ao^o), where 

Co = (J'/i^pAo) is the BCS coherence length. By equating the phase on 
the right-hand side to 27m, where n = 0,±1,±2,..., and using cosy = 
y^(l -I- cos 27)72 = EjAo, we find that the energy eigenvalues are given by 

The energy eigenvalue spectrum has a phase periodicity of 27 t. In the limit 
of a short junction {L Co), the left-hand side of (6.6) can be approximated 
by zero. In this case the two bound states Eq and EZi simply follow a cosine 
phase dependence: 

= +AoCos{cj)/2) , 

EZM = -^oCos((^/2) . (6.7) 

The states are degenerate at ^ = tt. 

For a finite length L comparable to Co, the eigenvalue spectrum is as 
shown in Fig. 6.2a. Now, the discrete states are also degenerate at <^ = 0. If 
the phase difference is increased, the degeneracy is lifted. A typical eigenvalue 
spectrum as a function of the phase difference for a long junction (L Co) 
is given in Fig. 6.2b. Similarly to the eigenvalues for electrons in a quantum 
well, the number of states is increased, while the energetic separation between 
adjacent levels is lowered. At low energies {E Aq) the energy dispersion 
is approximately linear [126]: 

En ('(') = [7r(2n -h 1) T </'] • (6.8) 

As we shall see in the next section, the states E^{(j)) carry a supercurrent 
in an SNS junction. In equilibrium the contribution of each of these chan- 
nels to the net supercurrent is determined by the Fermi distribution function 
fo{E). At zero temperature only the states with E < 0 contribute to the 
supercurrent, since they are fully occupied, whereas the states with E > 0 
are empty. Besides the bound states E^{cj)), part of the net supercurrent is 
also carried by continuous states with \E\ > Aq. Details of how the Joseph- 
son supercurrent in an SNS junction is calculated are given in the following 
section. 

6.1.2 Supercurrent 

in a General One-Dimensional SNS Junction 

In this section we shall treat the more general case of a one-dimensional 
SNS junction with barriers at the interfaces and a Fermi velocity mismatch 
between the superconductor and the normal conductor. The Andreev-level 
spectrum and the supercurrent can be calculated very elegantly by a trans- 
mission matrix method, as presented by Tang et al. [144]. An alternative 
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Fig. 6.2. Andreev levels of a short (T/^o = 2) (a) and long (T/^o = 20) (b) ideal 
SNS junction 



method to calculate the supercurrent would be to use the scattering-matrix 
approach [149, 150, 151]. 

Once again, a step- like superconducting pair potential Aq and a phase dif- 
ference (j) between the two superconducting electrodes are assumed (Fig. 6.3). 
The potential U{x) is given by the ^-shaped barriers and by the step due to 
the difference Uq between the Fermi energies between the normal conductor 
and the superconductor: 

U{x) = Uq [0{x) — 0{x — L)] + Z [<5(a:) -|- 5{x — L)] . (6.9) 
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Fig. 6.3. Spatial dependence of the pair potential A(x) and the potential U{x) in 
an SNS junction 
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6. Josephson Effect in Superconductor/2DEG Junctions 



For the calculation of the Josephson supercurrent, the transmission of in- 
coming electron-like particles and hole-like particles iFh through the junc- 
tion is described by a transfer matrix T(i/)) [144]; 



( 







( 6 . 10 ) 



T(())) can be decomposed into matrices representing the different steps dur- 
ing the transport. The transmission through the left (SN) and right (NS) 
boundaries is represented by the following two matrices, which contain the 
transmission coefficients given by (2.16d) and (2.18) for the generalized BTK 



model introduced in Sect. 


2.2.2: 










g-i0/4j/_ \ 




(6.11) 


I = 


1 + 
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e*0/4c'_ J 


, x = 0 , 


O = 


f e-‘'^/4c+ 


ei0/4g_ J 


, X = L . 


(6.12) 



After passing through the left SN interface, the electron and hole acquire 
phase shifts during the transfer through the normal conductor of kf.L and 
— k\^L, respectively, which can be expressed by the matrix P: 



P = 



0 




(6.13) 



An electron or hole can either be transferred directly through the N layer 
or be multiply reflected inside the normal conductor as in a resonant cavity. 
Each of the cycles can be described by a reflection at the NS interface S(</>/2), 
travel to the opposite interface P, a reflection at this interface S(— ^/2) and, 
finally, a transfer to the starting point P, as can be seen in Fig. 6.4. If this 
single cycle is represented by a matrix M(<^), the total transfer from the left 
to the right interface can be described by 



^M(0)" = [l-M(^)]-i. 

n—0 



(6.14) 



M((^) is composed of the matrices for the above-mentioned single steps, where 
the scattering matrices S for the Andreev reflection a± and the normal reflec- 
tion b± of an electron and a hole from the normal-conductor side are given 
by 



S(±^/2) = 



b+ e±‘'^/2a_ \ 
eTi-^/2a+ &_ ) • 



(6.15) 



Combining all components leads to 



M(((.) = PS(-()>/2)PS((/)/2) . 



(6.16) 



For the total transmission matrix T through the junction, the transmis- 
sion through the SN interface (6.11) and through the opposite NS interface 
(6.12) have to be taken into account: 
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Fig. 6.4. Transmission and reflection processes in a superconductor/normal con- 
ductor/superconductor structure 






C+(<^) D+icj,) 

D_{4>) c.{4>) 



0(<^/2)[l-M(<^)]-ipi(-<^/2). (6.17) 



The Josephson supercurrent /(</>) can be calculated by using the transmis- 
sion matrix T directly. However, for the practical computation of it is 
useful to separate the current into a component Idis, originating from the dis- 
crete Andreev levels inside the gap (|i5| < Zlg), and a component Icont, from 
the continuous spectrum outside the gap (|£1| > Aq). The discrete energy 
levels can be obtained from the poles of the transfer matrix T. By writing 






1 - 

det[l — M(</))] 



1 - 

T 



(6.18) 



where M* is the adjoint of M and F = det(l — M), it can be seen that the 
energy eigenvalues can be determined from the roots of F. Expanding F leads 
to 



F = 



( A 



COS <() - ^2p2 + 



+2{q^ — Z"^) cos [3 — 4Zgsin/3 



E 



(2p^ -I- 2p -I- 1) cos ( ~ 2(<7^ — Z^) cos/3 -I- AZqs\x\ (3 



+( 2 p+l)^Jl-(^) sinfaT. 



(6.19) 



where a = kp^LAoK^ — Uq) and (3 = 2fcpNT. In a and /3, (fee — fch) has 
been approximated by — Uq) and (fee -I- fch)/2 by fepNl P is defined 

as Z'^/r + {r — l)^/4r and q as Z'^/r -|- (1 — r^)/4r; r is the Fermi velocity 
mismatch given by (2.17). 
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At zero temperature (T = 0) only the discrete levels below the electro- 
chemical potential /r are occupied and contribute to the net supercurrent. The 
current /dis carried by the discrete levels, determined from (6.19), is obtained 
from 



.^dis — 



E 

n ,±,£'<0 



2e dA±(^) 

Ti d(f> 



( 6 . 20 ) 



The supercurrent /cont due to the continuous energy spectrum must be cal- 
culated directly from the transmission matrix T: 

2e 

leant = Y y i{E,(j))dE 

= I / <)>) - <^)] 

J-co l“0 ^ol 

- [Tt^niE, - Tt^niE, 0)] }dE 

n y_oo \ul - vl\ 

X [\ C + m ^ - |C+(-<^)p - + |C_(-<)>)P] d£(6.21) 

The continuous current is composed of a current carried by both electron-like 
and hole-like particles. In the above expression we have made use of the fact 
that the current direction can be inverted simply by reversing the sign of </>. 
The factors C±(^) are defined in (6.17). 

The general expressions (6.19) and (6.21) for the calculation of the total 
Josephson supercurrent are relatively complex. However, for some special 
cases expression (6.19) can be simplified significantly; i.e., for a clean junction 
with no barriers at the interface {Z = 0) and no Fermi velocity mismatch 
(t/o = 0), it reduces to (6.6). The general expression for bound states in a 
short junction (L ^ (/i — t/oj/fcpN^o), which can be obtained from (6.19), is 



A±(0) = ±Z\o 



/ cos^{(j)/2) + 2 p 2 + 2p- 2 q 2 + 2^2 

2p2 + 2p- 2g2 + 2Z‘^ + 1 



( 6 . 22 ) 



Here we have assigned E^ = Eq and E^ = E_^. Since 2p2 _|_ 2p — 2g2 is 
equal to 2^2, this equation reduces to 



Eb(4>) = ±^c 



cos2(<;i/2) -k4Z2 
4^2 -hi 



(6.23) 



which is the expression for the eigenvalues if a junction with J-shaped barriers 
but no Fermi velocity mismatch (r = 1) is considered. For the short-junction 
case the continuous current contribution is zero, so that only the current due 
to the discrete levels has to be considered for the calculation of the total 
supercurrent. In Fig. 6.5a the Andreev levels are plotted as a function of 
the phase difference </> for the case Z = 0, (6.7), as well as for Z = 0.5 and 
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Fig. 6.5. (a) Normalized Andreev levels and (b) Josephson current of a short 
junction (L ^ ^o) as a function of the phase difference ij> for Z = 0 {solid lines), 
Z = 0.5 {dashed lines) and Z — 1.0 {dash-dotted lines). The direction of the electron 
motion for the corresponding Andreev levels is denoted by arrows 



1.0; these were obtained by using (6.23). Owing to the additional normal 
reflection process for Z > 0, the degeneracy at <^ = tt in the quasiparticle 
spectrum is lifted, with the consequence that an energy gap appears. As 
can be seen in Fig. 6.5b, the supercurrent is reduced by the presence of a 
barrier at the SN interface. The maximum value of the supercurrent which 
can be obtained before the junction switches into the normal state is called 
the critical current. In Fig. 6.5b it can be seen that /(</)) gradually transfers 
from the clean-limit point contact behavior [24, 152] to a tunnel junction 
characteristic [I{(j)) oc sin</f] [153]. 

For a finite length of the normal-conducting region, more Andreev levels 
appear within the gap, as can be seen in Fig. 6.6a for a 450 nm long junction. 
The Andreev levels were determined by numerically calculating the roots of 
(6.19) for a given phase difference 4>. The degeneracy of the levels at (f> = 0, 
7T and 27 t for Z = 0 and r = 1 is lifted if Z > 0 or if r < 1.0. Since the E{(f>) 
dispersion is flattened for increasing Z and r, the Josephson current due to 
these discrete levels is also lowered according to (6.20). Here, r = 0.44 was 
chosen, as a typical value for an Nb-InGaAs/InP junction, as introduced in 
Sect. 5.2.1, with a typical carrier concentration of 7.1 x 10^^ cm“^ for the 
two-dimensional electron gas. In Fig. 6.6b, the discrete and the continuous 
contribution as well as the total supercurrent are plotted as a function of the 
phase difference ^ for Z = 0.5 and r = 0.44. Since now the junction has a 
finite length, the continuous Andreev-level spectrum also contributes to the 
net supercurrent. 
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Fig. 6.6. (a) Andreev levels for a 450 nm long contact with Z = 0.5, r = 1.0 
{solid lines), Z — 0.0, r = 0.44 {dash-dotted lines), and Z = 0.5, r = 0.44. {dashed 
lines), (b) Discrete, continuous and total Josephson supercurrent for a 450 nm long 
junction with Z = 0.5 and r = 0.44 



At finite temperature the occupation of the Andreev levels is not simply 
determined by a step-like distribution function but rather by a Fermi distri- 
bution function fo{E) broadened around jj,. As a consequence, Andreev levels 
above /x become occupied and levels below /x are partially emptied. In this 
case the supercurrent carried by the discrete bound states is given by 

= E 

n,± ^ 



The current contribution originating from the imbalance i{E, (p) in the elec- 
trical current per unit energy of the quasiparticles flowing in the continuum 
levels is given by 

-An 



I 

Jrnnt. — ^ 



7 

J A, 



i{E,(t>)fo{E)AE . 



(6.25) 



For a short junction with a two energy levels, as described by (6.23), the 
temperature dependence of the supercurrent can be calculated easily: 

2e Aq sin 4> 



ddis — 



^ 4v/(4Z2 -h l)(cos2 (j)/2 -h 4Z2) 
cAq sin (j) 



Tt 2y'(4Z2 -h l)(cos2 4,/2 -h 4Z2) 



[fo{E^) - fo{E+)] (6.26) 

tanh(ifg /2 A:bT) . 
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The current direction of the upper states at is opposite to that for the 
lower state at E^. As the temperature is increased E^ is filled more and 
more, while E^ is emptied at the same time. As a result the net supercurrent 
decreases with temperature. For junctions with a finite length, an exponential 
decrease of Ic with temperature, Ic oc exp(27rL), is derived. 



6.1.3 Supercurrent in a Junction 
with a Two-Dimensional Electron Gas 



The one-dimensional model, as introduced above, can be used very success- 
fully to elucidate the fundamental transport mechanisms in a superconduc- 
tor/normal-conductor Josephson junction. However, the experimentally in- 
vestigated junctions contain a normal conductor between the superconduct- 
ing electrodes which is extended in two dimensions. This implies that elec- 
trons and holes in the two-dimensional channel impinge at various angles 
onto the interface with the superconductor. Only the component of the wave 
vector normal to the interface is relevant for the phase-coherent transport 
mechanisms described in the preceding section. Thus, at a given Fermi wave 
number, a complete spectrum of wave vector components normal to the in- 
terface has to be taken into account. 

A model for a two-dimensional semiconducting channel between two su- 
perconducting electrodes was first introduced by Kresin [147]. An extension 
of this model, which includes 5-shaped barriers at the S/2DEG interface 
and a potential step due to the Fermi velocity mismatch (see Fig. 6.3), was 
developed by Chrestin et al. [148]. Here, the confining potential of the two- 
dimensional channel was modeled by a triangular well [70]. This shape ap- 
proximates the potential of an inversion layer at the surface of InAs or in 
a heteroj unction. If a single subband is occupied in this quantum well, the 
current through the junction can be expressed by 



cA:bT 1 sin (j) 

^ '“27ry_^ (1/2) cos (j) + f{ky,ujn) ^ 



(6.27) 



where 
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(6.28) 
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and 
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(6.30) 



Here, a is defined as m* /me and p = (fc|g — k^y^^ is the wave vector in the 
superconductor perpendicular to the interface. is given by y/oj/^ + Aq, 
where tUn = (2n + l)7rfcBT are the Matsubara frequencies. The wave vector 
k^ within the semiconductor is 







fc2±i 



2m* 

~W 






(6.31) 



where fcpN = y (2m*/?i^)(/x —Un — Eq) \s the Fermi wave vector in the semi- 
conductor. Eq is the energy level of the first subband in the 2DEG and Uq 
is the conduction band offset between the superconductor and the semicon- 
ductor. 

As we shall see in Sect. 6.2.1, the model of Chrestin et al. [148] can be 
employed to explain the temperature dependence of the critical current. In 
addition, this model is of interest if the characteristics of junctions with a 
gate on top of the two-dimensional electron gas are to be described. By ap- 
plying a gate voltage, the carrier density ri 2 D in the channel can be adjusted. 
This directly affects the Fermi wave vector as given by (3.7) in (6.31) and 
therefore also the critical current. If the carrier concentration is increased, 
the critical current is on the average increased, if superimposed oscillations 
due to resonances are neglected [148]. 

The coherent transport of an SNS structure including a 5 barrier at both 
interfaces was also discussed in the clean limit by Brinkman and Golubov 
[154]. The authors could show that for a thin normal layer, the junction is 
in the coherent regime, where the supercurrent is proportional to the angle- 
averaged transmission probability. This is expected for a thin interlayer (N), 
where the transmission resonances are broad. In contrast, for a wide inter- 
layer the supercurrent was found to be proportional to the square of the 
angle-averaged transmission probability. This can be explained by sequen- 
tial, incoherent tunneling through the two 5 barriers. The results given by 
Brinkman and Golubov [154] were obtained for a three-dimensional structure 
but they can also easily applied to the two-dimensional case. 



6.1.4 Supercurrent in a Junction with Proximity Layers 

In Sect. 5.1.4 it was shown that a proximity layer (N) between the super- 
conductor and a constriction (c) in a junction leads to additional features 
in the subharmonic gap structures. For the interpretation of the measured 
current-voltage characteristics of some particular superconducting junctions. 
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it is necessary to include a proximity layer between the superconductor and 
the 2DEG. Those models which are able to describe this situation usually 
consider the two-dimensional electron gas as a constriction. This is justified 
by the large difference in the number of states between the three-dimensional 
superconductor and the two-dimensional electron gas. 

A superconducting junction with two superconducting electrodes con- 
nected by a constriction (ScS) was first discussed by Kulik and Omel’yanchuk 
(KO model) [152, 155]. The more general case with two different supercon- 
ductors (ScS'), was considered by Zaitsev [50]. On the basis on this theory, 
Golubov and Kupriyanov succeeded in obtaining an expression for the su- 
percurrent in an SNcNS junction with a constriction shorter than [156]. 
The Green’s functions suitable for this problem have already been introduced 
in Sect. 2.2.3. Inserting these functions into the Zaitsev equation yields the 
following expression for the supercurrent: 



The function ^ is defined as ^{uJn) = Z\/[l -|- (Z\^ -|- w^j^/^ye/Tr/cBTc]. D is 
the normal transmission probability of the constriction and 7 b represents the 
influence of the SN interface transparency as defined by (2.29). 

In Fig. 6.7, the temperature dependence of elcRjsi / t^Tc is shown for dif- 
ferent values of the barrier strength at the SN boundary, 7 b, and of the 
suppression parameter 7 [156]. The temperature dependence is plotted in 
terms of the product IcRn, normalized to the critical temperature Tc. This 
IcRn product, also referred to as the characteristic voltage Vc, is a measure 
of the superconducting coupling strength. The maximum value of « 1.76 is 
obtained for the ideal case of a ballistic ScS junction (KO model) [152], this is 
twice the value obtained for an ideal superconducting tunnel junction [153]. 
As can be seen in Fig. 6.7, the critical current is lowered if a proximity layer 
is introduced between the constriction and the superconductor. Furthermore, 
the shape of the curve changes from convex to concave as the superconducting 
electrodes are increasingly decoupled. 



6.2 Measurements of the Superconductive State 

In the following part of the chapter, experimental results on superconduc- 
tor /two-dimensional-electron-gas Josephson junctions are discussed. First the 
I~V characteristics are analyzed in detail. Subsequently, further properties, 
e.g. interference effects due to a magnetic field, which depend on the phase 
coherence are covered. The last part deals with the light sensitivity of the 
semiconducting material. As we shall see, this effect can be employed to 
control the critical current in the junction. 




(6.32) 
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Fig. 6.7. Normalized critical current cIcRn /ttTc of an SNcNS junction as a function 
of T/Tc with 7 b (a) and 7 (b) varied as a parameter [156]. The characteristics of an 
ScS junction obtained from the Kulik-Omel’yanchuk model are also shown (KO) 
for comparison [152] 





6.2.1 Current Voltage Characteristics 

One of the basic methods to characterize a superconductor/two-dimensional- 
electron-gas junction is the measurement of the current-voltage character- 
istic. From these measurements, information on the junction properties can 
be obtained. Let us begin with a typical I-V characteristic of a Nb/p-InAs 
step edge junction, as shown in Fig. 6.8 [140]. It can be clearly seen that a 
Josephson supercurrent flows through the junction with a maximum value of 
7c = 220 pA. In contrast to Josephson tunnel junctions, the current- volt age 
characteristic is ohmic, which means linear, in the voltage-carrying state. At 
voltages above 2Z\o/e a normal-state resistance 7 ?n of 4.5 Q is measured, 
which results in a characteristic voltage 14 = Ici?N of 990 pV. Compared 
with the earlier results on Nb/p-InAs junctions [8, 30, 92] the IcRn prod- 
uct obtained here is larger by almost one order of magnitude and thus is 
relatively close to the maximum value (4/7r)eZ\o given by theory [152]. This 
improvement can be attributed to the optimized cleaning procedure prior to 
the deposition of the Nb electrode. An even larger /ci?N product of 1.35 meV 
at T = 2 K was obtained for an Nb/p-InAs junction fabricated by anodic 
oxidation, with a layout as shown in Fig. 3.5b [29]. The corresponding crit- 
ical current has a value of 870 pA. The very large IcRn product can be 
attributed to the small distance of only 20 nm between the Nb electrodes 
and the refined Ar+ sputter cleaning. The latter leads to an improved con- 
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Voltage (mV) 



Fig. 6.8. Current -voltage characteristic of a 30 | 0 ,m wide Nb/p-InAs step junction 
at a temperature of 4.2 K. The step height was « 140 nm. The geometry of the 
junction is depicted in Fig. .i..la [140] 



pling of the superconductor to the two-dimensional channel at the surface. 
In Fig. 6.9, the temperature dependence of the critical current is plotted for 
a 20 nm and a 40 nm long junction. In order to explain the decrease of Ic 
with temperature, the model for an SNcNS junction, as introduced in Sect. 
6.1.4, can be referred to. The reason for doing so is that the two-dimensional 
channel underneath the Nb electrodes can be regarded as a diffusive normal 
conductor N, which serves as the proximity layer within this model. The as- 
sumption that this layer is indeed diffusive can be justified by the generation 
of defects due to the Ar+ sputter cleaning prior to the deposition of the Nb 
electrodes. Owing to the proximity effect, a gap is induced in these N layers. 
The two-dimensional electron gas, serving as a bridge between the N layers 
covered by the Nb electrodes, is regarded here as the constriction. As can be 
seen in Fig. 6.9, the temperature dependence of the critical current can be 
very well explained by the SNcNS model, if a finite barrier between the su- 
perconducting electrodes and the 2DEG underneath (yp > 0) and a reduced 
transparency of the constriction (£><!) are assumed [29]. 

The current-voltage characteristic of a typical junction with a buried two- 
dimensional electron gas in a strained InGaAs/InP heterojunction is plotted 
in Fig. 6.10. The junction has a length of 450 nm and a width of 6 |4m (sample 
B in Table 6.1). At 0.3 K a critical current Ic of 3.8 pA is obtained. For bias 
voltages larger than 3 mV, a normal-state resistance i?N of 50 Q is measured, 
which results in a characteristic voltage Vc = IcRn of 190 pV. 

At low temperatures (< 1.3 K) the I-V characteristic shows a pronounced 
hysteresis. For sample B the return current /r, which is the current where the 
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Fig. 6.9. Critical current as a function of temperature for a 20 nm and a 40 nm long 
planar Nb/p-InAs junction. The experimental values are compared with theoretical 
results (solid line) for the values of Tc and _Rn given in the graph. The inset shows 
the I-V characteristic of a 20 nm long junction at 1.8 K (solid line) and the ohmic 
line corresponding to = 1.55 £2 (dotted line). The hysteresis can be explained 
by the capacitance caused by the overlapping parts of the electrodes [29] 



junction switches back from the normal-conducting to the superconducting 
state, is 2.4 pA (Fig. 6.10). In order to detect the hysteresis without any 
distortion, extensive filtering is required [d8, 157]. Here, tt and RC filters 
were used at 300 K and 4.2 K, respectively. Additional filtering was obtained 
from by Thermocoax Philips cables directly connected to the sample at low 
temperature [158]. 

In Table 6.1, Ic, Rn and Vc are given for a set of junctions with different 
geometries. Since the characteristic voltage is a measure of the supercon- 
ducting coupling strength, it is reduced if a barrier is present at the S/2DEG 
interface or if the separation of the superconducting electrodes is enlarged. 
The latter effect can be observed here. For the larger separations, 600 nm 
(sample E) and 800 nm (sample F), the characteristic voltages are 139 |lV 
and 112 pV, respectively, which are lower than the values obtained for the 450 
nm long junctions. However, it can be seen in Table 6.1 that Vc varies among 
the 450 nm long junctions. This is due to the fact that the average value of 
the transparency, as well as the transparency, along a particular S/2DEG in- 
terface, can vary. The characteristic voltages of the junctions discussed here 
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Fig. 6.10. Current -voltage characteristic of a 450 nm long, 6 |lm wide Nb- 
InGaAs/InP junction at a temperature of 0.3 K. The upward and downward 
branches are marked by arrows. The inset shows the I-V characteristic for larger 
voltages. The excess current 7exc is obtained by a linear extrapolation {dashed line) 
[157] 

Table 6.1. Characteristic data for a set of Nb-InCaAs/InP-Nb junctions with 
different geometries 



Sample 


Length 

(nm) 


Width 

(|lm) 


Nb electrode 

separation 

(nm) 


Ic 

(|iA) 


Rn 

(D) 


IcRn 

(HV) 


A 


300 


6 


300 


3.3 


43 


142 


B 


450 


6 


350 


3.8 


50 


190 


C 


450 


6 


450 


2.7 


69 


186 


D 


450 


16 


350 


13.0 


20 


260 


E 


600 


8 


400 


4.1 


34 


139 


F 


800 


8 


600 


3.4 


33 


112 



are of the same order of magnitude as 14 of comparable junctions with a 
heterostructure quantum well as the conductive channel [37, 38, 98]. 

In Fig. 6.11, the normalized critical current Ic/Ico and the normalized 
return current Ir/Iro of samples B and C are plotted as a function of tem- 
perature. Ico and Jro are the critical current and return current, respectively. 



82 6. Josephson Effect in Superconductor/2DEG Junctions 

at 0.3 K. Although the absolute values of Ic and A differ for the two sam- 
ples, the temperature dependence is almost identical. It can be seen that 
the critical current increases with decreasing temperature. However, at very 
low temperatures (< 0.6 K) tends to saturate. The I-V characteristic is 
hysteretic for temperatures below 1.3 K so that the junction returns to the 
resistive state at a lower current A than the critical current. It is interesting 
to notice that I^/Iro remains almost constant with temperature. 




Temperature (K) 



Fig. 6.11. Normalized critical current Ic/IcO and normalized return current /r//ro 
[157] as a function of temperature for samples B and C. The solid line shows the 
theoretical results obtained from the model of Chrestin et al. [148] 



In order to explain the experimentally measured temperature dependence 
of Ic, the model of Chrestin et al. [148], as introduced in Sect. 6.1.3, can be 
used. The calculated critical current as a function of temperature is plotted 
in Fig. 6.11. For temperatures below « 2.4 K, very good agreement with the 
experimentally determined critical current is obtained if a Z parameter of 
1.0 is assumed. For larger temperatures a minor deviation develops, which 
can be explained by uncertainties in the experimental determination of Ic 
due to the thermal smearing of the I-V characteristics. The Z parameter 
was estimated from the normalized excess current elexcR^i /‘^o, following the 
approach introduced by Flensberg et al. [132]. For the sheet electron con- 
centration, which determines the Fermi wave vector in the 2DEG, a value of 
7.05 X 10^^ cm“^ was assumed. 

From the saturation of R at low temperatures, the transport mechanism 
in the junction can be classified. In Sect. 6.1.2 the spectrum of the Andreev 
levels was calculated for an SNS structure which contains a barrier at the 
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interface and a potential step in the normal-conducting regioir. As a coir- 
sequence, the degeneracy of the Andreev levels at the chemical potential is 
lifted and a gap is formed instead. At very low temperatures only the An- 
dreev levels below the chemical potential are occupied, so that the critical 
current remains constant. Only if the temperature is raised to a value where 
Andreev levels above fi are thermally occupied will start to decrease with 
iircreasing temperature. The 1^ vs. T characteristic discussed here shows ex- 
actly the behavior described above. It cair thus be concluded that an eirergy 
gap is present in the Andreev level spectrum. A similar characteristic was 
also observed by Takayairagi and Akazaki [159]. Iir coirtrast, Mur et al. [38] 
found no saturation iir the Ic vs. T measurements but rather an exponential 
law, Ic(T) = lc(0) exp(— T/T q), that match the experimental values. They at- 
tributed their finding to transport through a few resonant transport channels 
induced by disorder with almost ideal transmission probability [160, 161]. 
Although their semiconductor material is similar to the InGaAs/InP het- 
erostructure discussed here, the technique of making contact between the Nb 
electrodes and the 2DEG differs substantially. Iir Mur et al.’s junctions the 
channel layer is contacted from the top, which might lead to contacts with a 
higher transmission probability owing to multiple Andreev reflections within 
the 2DEG channel layer. 



6.2.2 Effect of a Magnetic Field 



By means of Ic measurements as a function of an external magnetic field B, 
more detailed information about the supercurrent transport in the junction 
may be gained. Owing to the effect of the magnetic flux penetrating the 
junction, an interference pattern develops in the critical current. For the 
simplest case, with a short junction, a regular Fraunhofer-like interference 
pattern is expected for a homogeneous supercurrent distribution [162, 163]; 



Ic = 



IcO 



sin(7T^/^o) 

TT^/^O 



(6.33) 



Here, Ico is the critical current at zero magnetic field, <P is the effective mag- 
netic flux in the junction and = h/2e is the magnetic flux quantum. As 
can be seen in Fig. 6.12a, the assumptioir of such a constant supercurrent dis- 
tribution for the simulatioir fits the Ic vs. B interference pattern of a 450 nm 
long, 6 pm wide Nb-InGaAs/InP junction (sample B, see Table 6.1) almost 
ideally. 

For a 450 nm long, 16 pm wide junction (sample D), a step-shaped current 
density distribution (Fig. 6.12b, inset) has been assumed in order to explaiir 
the more complex measured interference pattern. As mentioned above, the 
variation in the current distribution is presumably due to an inhomogeneous 
interface. The experimental curve can be fitted reasonably by using the fol- 
lowing expression [162]; 



84 



6. Josephson Effect in Superconductor/2DEG Junctions 



B (mT) B (mT) 





Fig. 6.12. /c at 0.5 K as a function of B for samples B (a) and D (b) (o). The 
geometry of both samples is given in Table 6.1. The full line shows a fit as a function 
of using a current density distribution as depicted in the corresponding insets 






b — 2ab + 2a 
sin(7ra<?/<?o) 



b — 2ab sin[7r(l — 2a)(P/$o] 
[7t(1 - 2a)^/^o] 

cos[7t( 1 — a)<P/^o] • 



(6.34) 



Here, the outer parts of width aW carry the current density ico, whereas the 
inner part of width (1 — 2a)W carries the current density bico; W is the total 
width of the sample. For junction D, a larger supercurrent density at both 
ends of the junction, expressed by the parameters a = 0.2 and b = 0.62, was 
assumed (Table 6.2). Considering the simplicity of the underlying model, the 
agreement is relatively good. In Table 6.2 the fitting parameters of two other 
samples, C and E, are summarized, in addition. As can be seen here, the 
current density ico of junction D is considerably lower than for the 450 nm 
long sample. This can be explained by the smaller slope of the Andreev levels 
if the electrode separation is increased (see Sect. 6.1.2). According to (6.20), 
the Josephson current is determined partially by the slope of the Andreev 
levels as a function of (f). 

Because of the planar geometry of the junctions, the external magnetic 
field is focused into the semiconductor area owing to the shielding of the 
superconducting electrodes. Experimentally, the focusing factors (i?eff/i3)exp 
were obtained by matching the minima of the measured interference pattern 
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Table 6.2. Results of the measurements in a magnetic field; a and b are the htting 
parameters in (6.34) and a = 0 means that a homogeneous current distribution 
is present. (Beff/B)exp and {BeS / B)theo are the experimentally determined and 
theoretical focusing factors, respectively, of the magnetic field 



Sample 


a 


b 


*c0 

(pA) 


\ -LJ / exp 


(t) 

\ D J theo 


B 


0 


1 


80 


7.05 


8.55 


C 


0.23 


0.3 


62 


6.66 


8.55 


D 


0.20 


0.62 


82 


16.85 


16.44 


E 


0 


1 


44 


8.27 


8.21 



with the corresponding multiples of For the calculation of the effective 
magnetic field Beff, the London penetration depth (« 20 nm) was taken 
into account for the area penetrated by the magnetic field. A comparison of 
theoretically determined focusing factors (i?eff/.Bo)theo, calculated from an 
expression given by Gu et al. [164], shows good agreement with the experi- 
mentally obtained values (Table 6.2). 

The magnetic field dependence of Ic of Nb-InAs/AlSb junctions with a 
width W comparable to the length L has been investigated by Heida et al. 
[87]. These authors observed a large deviation from the periodicity (6.33) 
and interpreted their results in terms of mesoscopic phase effects. In con- 
trast to tunnel junctions or junctions with L ^ W , here all trajectories from 
one electrode (the left electrode in Fig. 6.13) to the opposite one have to be 
considered, rather than only the trajectories perpendicular to the interface 
{x\ — > X 2 in Fig. 6.13) as in the case of a short junction or a tunnel junc- 
tion. In a heuristic approach using a quasiclassical description, the following 







w 






Fig. 6.13. Josephson junction with a width W comparable to its length L. The 
supercurrent can flow from xi to any X 2 at the right electrode [87] (© (1998) by 
the American Physical Society) 
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expression was derived for the magnetic-field dependence of /c by averaging 
over all possible trajectories [87]; 



As can be inferred from this equation, a periodicity of 2^o is expected in this 
case, in contrast to the periodicity of junctions with L <C W. A period of 
2^q was indeed experimentally observed for junctions where L is comparable 
to W [87]. 

6.2.3 Hysteresis 

The I-V characteristic plotted in Fig. 6.10 reveals that for low temperatures 
a hysteresis occurs. Furthermore, the temperature-dependent measurements 
show that the return current A is almost constant for T < 1.0 K. The hystere- 
sis originates from a negative differential conductance at low bias voltages. 
If an ohmic shunt resistor i?g = 6.8 Q [Fig. 6.14 (inset)] is connected to the 
junction, the range of negative conductance can be resolved. Apart from this 
shunt resistor, the contact resistances (10 Q each) must also be taken into 
account, since i?s is not directly connected to the junction but rather is fixed 
to the bond pads. The value of the resistance i?s was chosen sufficiently low 
that a positive differential conductance was obtained for the complete circuit. 
The corrected I-V characteristics, measured on the shunted junction (sam- 
ple B) for various temperatures, are shown in Fig. 6.14 [165]. Here, a slightly 
larger critical current of 4.1 pA was measured, which can be attributed to 
changes due to the multiple thermal cycles the sample had to undergo be- 
fore this measurement. Up to a temperature of 1.3 K a negative differential 
conductance is observed. 

In contrast to Josephson tunnel junctions, the hysteresis cannot be ex- 
plained by a shunt capacitance within the “resistively capacitively shunted 
junction” (RCSJ) model [162, 166, 167]. An estimation of the junction capac- 
itance shows that it is two orders of magnitude smaller than the capacitance 
necessary to explain the hysteresis within the RCSJ model. In addition, 
is constant for temperatures T < 1.0, in contradiction to the predictions of 
the RCSJ model. For this type of junction the occurrence of the hysteresis 
can be explained by an effective heating of the electrons in the 2DEG, as 
pointed out by Heida [168]. This heating effect originates from a broadening 
of the electron distribution function in the 2DEG due to multiple Andreev 
reflections in the voltage-biased regime. The broadening of the distribution 
function can be seen, for example, in simulations using the OTBK model 
when compared with the corresponding equilibrium Fermi distribution func- 
tion, as described in detail in Sect. 5.1.3. Following the approach discussed 
in Sect. 6.1, the resulting supercurrent in an S/2DEG/S junction is deter- 
mined by the occupation of discrete and continuous Andreev levels. In the 
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Fig. 6.14. Current-voltage characteristics measured at various temperatures. The 
inset shows the circuit with the shunt resistor Rs and the contact resistances 7?c 



[165] 



case of the discrete levels the corresponding part of the net Josephson cur- 
rent is provided by levels carrying a supercurrent alternately in the forward 
and the reverse direction. A broadening of the distribution function leads 
to occupation of formerly unoccupied discrete states and thus the maximum 
obtainable supercurrent is reduced to A. In the nonshunted case the I~V 
characteristic switches back into the superconductive state at a current A, 
which corresponds to the critical current at this effective temperature. It is 
interesting to notice that the effective heating is relatively independent of the 
bath temperature, since the minimum in the I-V characteristic, as shown in 
Fig. 6.14, is found at the same current for various temperatures. 

In principle a negative differential conductance can also result from the 
conductance enhancement due to multiple Andreev reflections, as described 
within the KGN model [124]. For a decreasing bias voltage, more and more 
Andreev reflection processes occur. Since each Andreev reflection event in- 
creases the conductance, an increase of the total current and thus a negative 
differential conductance may be predicted for decreasing voltage. However, 
in order to obtain a noticeable conductance enhancement the interface trans- 
parency needs to be sufficiently large. This is in contradiction to our find- 
ing that a relatively large Z parameter of about 1.0 had to be assumed in 
order to describe the temperature dependence of Ic- A further indication 
that the measured negative differential conductance cannot be explained ad- 
equately within the KGN model is given by measurements in a magnetic held 
(Fig. 6.15). A small magnetic held of 0.19 mT is sufficient to suppress the neg- 
ative differential conductance completely. This is not expected in the KGN 
model, since here the transport is described by the nonphase-coherent motion 
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of Gaussian wave packets, which should be insensitive to a small magnetic 
field. 




Fig. 6.15. Current and differential conductance of sample B as a function of voltage 
for B — 0 {solid lines) and B = 0.19 niT {dashed lines) [165]. The inset shows the 
Shapiro steps of an 800 nm long junction (sample F) for 12.7 GHz radiation at 0.7 K. 
The radiation power was increased step by step from an attenuation of —15.7 dB 
to -11.3 dB 



6.2.4 Self-Detection of Josephson Radiation 

A closer inspection of the differential conductance plotted in Fig. 6.15 shows 
that for i? = 0, modulations are superimposed on the AI /AV curve. From the 
absence of these modulations at B = 0.19 mT, it can be directly concluded 
that these features are not due to subharmonic gap structures. A more likely 
explanation is that these oscillations are due to the AC Josephson effect. This 
conclusion is supported by the measurement depicted in Fig. 6.16a. 

At zero magnetic field the differential conductance shows peaks at 140 pV 
and at 230 pV, indicated by the concentric contour lines. If the magnetic field 
is increased these maxima vanish at first and then reappear at a field of 0.12 
mT. In particular, at a fixed voltage of 140 pV the differential conductance 
oscillates periodically with the magnetic field. The plot of AV/AI as a func- 
tion of B and I (Fig. 6.16b) reveals that the structures are connected to 
the Josephson effect since they are only observable if the Ic{B) Fraunhofer 
pattern has a local maximum. The superconducting state corresponds in this 
plot to the black area, where AV/AI = 0. 

A possible physical origin of the structures in the AI/AV curves is 
a self-detection of Josephson radiation. Similarly to the coupling of an 



B (mT) 



6.2 Measurements of the Superconductive State 89 



a) dl/dV b) dV/dl 




Voltage (mV) Current (yttA) 

Fig. 6.16. (a) Differential conductance of sample B as a function of the applied 
voltage and the magnetic field [165]. The contour lines trace constant values of 
dl /dV , while the absolute value is represented by a gray scale. In (b) the differential 
resistance dV/dl is plotted as a function of current and magnetic field. The black 
area, which represents a supercurrent, is identical to the Fraunhofer interference 
pattern 



Nb/2DEG/Nb junction to an external microwave radiation, as shown by 
an example in Fig. 6.15 (inset), the junction can also couple resonantly to 
its own AC Josephson radiation [168]. From the dl/dV plot for i? = 0 in 
Fig. 6.15, a peak separation of about 90 pV can be extracted, which results 
in a frequency of / = 2eV jh = 43.5 GHz. 

6.2.5 Light Exposure 

A unique property of superconductor/semiconductor junctions is that /c and 
i?N can be controlled by light exposure. The first measurements using this 
concept go back to the work of Giaever [169] and Andreozzi et al. [170], who 
studied junctions based on a CdS light-sensitive interlayer. Another normal- 
conductive material, which has also proved its suitability for light-sensitive 
weak- link junctions, is Pbi_a,Sna;Te [171]. In many IIHV semiconductor het- 



90 



6. Josephson Effect in Superconductor/2DEG Junctions 



erostructures a persistent photoconductance is observed at low temperatures 
[172, 173]. Persistent photoconductance means that the photogenerated car- 
riers remain in the channel layer even after the light source is switched off. 

The increase of Ic relative to the initial value IcO of an Nb-InGaAs/InP 
junction after light exposure is shown in Fig. 6.17a. Here, a red (665 nm) 
GaAs light-emitting diode (LED) operating at room temperature was used 
as a light source, with a plastic fiber guiding the light down to the sample 
[174]. In order to increase the carrier concentration in the 2DEG stepwise, the 
LED was operated in a single-pulse mode with a pulse duration of 5 ms. The 
initial critical current Ico of the 300 nm long, 6 pm wide junction was 3.7 pA. 
First, the normalized critical current Ic/Ico increases very steeply for up to 15 
light pulses. However, further exposure leads to a slower increase of Ic/Ico to 
a final relative increase of 12% after 90 accumulated light pulses. The critical 
current can be reset to its initial value by increasing the temperature above 
100 K. 




1.0 



0 . 9 ^ 
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Fig. 6.17. (a) Relative increase of the critical current 7c//co and the normal-state 
resistance 7 ?n/Ano of a 300 nm long, 6 pm wide Nb-InGaAs/InP junction as a 
function of the number of light pulses, (b) Gritical current 7c and normal-state 
resistance as a function of the LED bias current [174] 



At the same time as the critical current increases, the normal-state resis- 
tance i?N decreases with light exposure, as can be seen in Fig. 6.17a. Again 
the largest effect is measured for the first 15 pulses, while further light pulses 
lead to a slower decrease of 7 ?n- A relative decrease of about 16% was ob- 
served after 90 pulses. 
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In order to relate the increase of /c and the decrease of i?N to the proper- 
ties of the two-dimensional channel layer, Shubnikov-de Haas measurements 
on Hall bars were performed under pulsed illumination. Since for these mea- 
surements at high magnetic fields a different sample holder had to be used, 
the dose per light pulse was different in this case. However, Shubnikov-de 
Haas effect measurements allowed us to determine the sheet electron concen- 
tration 7T.2D and, at the same time, the sheet resistance Q 2 t> at zero magnetic 
field as a function of the number of light pulses (Fig. 6.18). The sheet elec- 
tron concentration increases from an initial value of 7.35 x 10^^ cm“^ to 
a value of 8.45 x 10^^ cm“^ after four light pulses. A further illumination 
does not change the sheet carrier concentration in the 2DEG. As can be seen 
in Fig. 6.18 (inset), a conductive bypass channel appears, indicated by the 
parabolic increase of the curve measured after seven light pulses. For less 
than four light pulses, no bypass channel was observed. 




number of light pulses 



Fig. 6.18. Sheet electron concentration ti 2 d (•) and specific sheet resistance Q 2 d 
(v) as a function of the number of light pulses. The measurement was performed 
at 0.3 K. In the inset the Shnbnikov-de Haas oscillations are plotted before illumi- 
nation {lower curve) and after seven light pulses {upper curve) [174] 



The photogenerated carriers are due to an excitation across the band gap 
of InP, as found by Kane et al. [175], from the wavelength dependence of 
the persistent photoconductance. Electrons and holes are separated by the 
built-in electric field of the heteroj unction. According to Kane et al. [175], 
the holes are trapped in the InP buffer layer. Up to a certain light dose the 
electrons are collected in the quantum well, whereas for even higher doses the 
electrons are found in the InP dopant layer. The latter effect is responsible 
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for the bypass channel. The sheet resistance q^d plotted in Fig. 6.18 first 
shows a steeper decrease, followed by a saturation in accordance with the 
measurement of i?N in Fig. 6.17a. By matching the measured drop of Q 2 t> to 
the drop of the normal-state resistance i?N of the junctions, it is now possible 
to gain information about the corresponding sheet electron concentration in 
the junction. 

As mentioned above, the intensity of the light pulses was different for the 
measurements on the Hall bars and on the Josephson junctions. However, it 
is possible to match the measured drop of the normal-state resistance i?N 
with that of the specific resistance g 20 by rescaling the intensity of the light 
pulses. By this method, information about the electron concentration ri 2 D in 
the junction can be gained. 

Referring to a simple picture, the increase of the critical current with 
increasing carrier concentration can be understood from the increase of the 
coherence length = hvp / . If the number of carriers in the channel 
is increased, the Fermi velocity wp = h^/2Tm2T>/'m* is increased, leading to 
a larger value of ^n- For a more detailed analysis, the model of Chrestin et 
al. [148] can be used (Sect. 6.1.3). However, owing to interference effects, the 
calculated critical current has an oscillating component as ri 2 D is increased 
[27]. No clear indication for this oscillation is seen in the measurements, pre- 
sumably because of fluctuations of the effective Nb electrode separation along 
the junction. Neglecting the oscillatory part of Ic in the model, an increase 
of about 21% was estimated for a 300 nm wide junction when the carrier 
concentration was increased from 7.35 x 10^^ cm“^ to 8.45 x 10^^ cm“^. For 
the calculations, a Z parameter of 1.0 was assumed in accordance with the 
simulations presented in Sect. 6.2.1. A possible explanation for this discrep- 
ancy in the increase of the critical current might be that a part of the 2DEG 
channel is covered by the superconducting electrodes so that the interface 
region is shadowed. The separation of the Nb electrodes is smaller than the 
distance L between the Nb/2DEG interfaces. Therefore, fewer photogener- 
ated carriers are present in this covered area. While the steep increase of 
the critical current can be explained by additional carriers in the channel 
layer, the bypass channel observed in the Shubnikov~de Haas measurements 
for more than four pulses is presumably responsible for the slower increase 

of Ic/IcO- 

The critical current can also be controlled by a continuous illumination, 
as shown in Fig. 6.17. Starting at an initial value of Ic = 3.8 pA, the critical 
current is increased by 24% to a final value of 4.7 pA. Similarly to the ex- 
periments with the pulsed light source, the normal-state resistance decreases. 
Since a conductive bypass channel is already observed after a few light pulses, 
the carriers excited under continuous light exposure will probably also be lo- 
cated in a bypass channel. The mechanism leading to a change of the critical 
current under continuous light exposure can in principle be used as a light- 
sensitive switch [174]. 
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In an InGaAs/InP heterostructure, the carrier concentration can only be 
increased, not decreased, by illumination. If a GaSb/AlSb/InAs heterostruc- 
ture is used as a two-dimensional electron gas instead, the critical current 
can be switched back and forth between two stable states by using two dif- 
ferent photon wavelengths. The underlying mechanism can be explained by 
referring to the band profile shown in Fig. 6.19. The heterostructure used for 




Fig. 6.19. Band profile of a GaSb/AlSb/InAs heterojunction. CB(J') and CB(A) 
are the conduction band edges at the F and X valleys, respectively. VB(T') is the 
valence band edge at the F point. If the layer system is illuminated by a 632 nm 
source, electron-hole pairs are generated in the GaSb layer. The holes are dragged 
to the InAs well by the built-in electric field. By illuminating with a 1.3 |tm source, 
deep donors located in the AlSb barrier are activated [176] 



this study is similar to the layer system depicted in Fig. 3.4a. In this case, the 
InAs channel layer has a thickness of 12.5 nm, while the upper AlSb barrier 
is only 50 nm thick. 

As depicted in Fig. 6.19, illumination with infrared light with a wave- 
length of about 1.3 |lm ionizes deep donors inside the AlSb barrier, which 
are transferred to the InAs quantum well. The carrier concentration in the 
quantum well is increased by this mechanism, similarly to the InGaAs/InP 
heterostructure discussed above. If the sample is illuminated by light of a 
wavelength shorter than 950 nm, electron-hole pairs are generated in the 
GaSb cap layer. As can be seen in Fig. 6.19, owing to the built-in electrical 
potential, the holes are dragged to the InAs well, where they recombine with 
part of the electrons accumulated in the quantum well [177]. As a result, the 
carrier concentration in the quantum well is decreased by this mechanism. 

The effect of illumination with two different wavelengths on a Josephson 
junction can be seen in Fig. 6.20. During illumination with a wavelength of 
632 nm, the critical current Ic decreases from an initially higher value to 
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Fig. 6.20. Current -voltage characteristics for an illumination sequence with 
632 nm and 1.3 |tm wavelength light sources and intermediate stages where the 
light source is off. (a) I-V characteristics for a large bias current range, (b) Details 
showing the characteristics close to zero voltage. Both the normal-state resistance, 
given by the slope in (a) , and the critical current are changed persistently by light 
illumination [176] 



7.4 |J,A. Simultaneously, the normal-state resistance 7 ?n increases. As men- 
tioned above, this effect can be explained by a decrease of the carrier con- 
centration in the InAs well due to hole transfer from the GaSb layer. After 
the light source is switched off, the normal-state resistance decrease slightly 
owing to some carriers trapped in the InAs layer, as can be seen in Fig. 6.20 
(the second curve in each part of the figure) [176]. At the same time Ic in- 
creases by almost a factor of three. When the 1.3 |im light source is switched 
on, is further decreased owing to an excitation of deep donors in the AlSb 
barrier. The carrier concentration in the well is increased by this effect. How- 
ever, because of the detrimental effects of Cooper pair breaking and heating, 
the critical current remains approximately at the previous value. Only if the 
1.3 |im source is finally switched off does Ic increase (the fourth curve in each 
part of the figure). The sequence as described above can be employed as a 
switch that cycles between two stable states [176]. 



7. Gate-Controlled Superconductor/ 
Two-Dimensional-Electron-Gas Junctions 



Compared with metals, semiconductors have the great advantage that the 
electron concentration can be changed very easily by means of a gate elec- 
trode. As we discussed in the preceding chapter, where light exposure was 
employed, a change of the electron concentration has an immediate effect 
on the superconducting properties, i.e. the critical current. The first struc- 
ture discussed below, which utilizes this effect, is the Josephson field effect 
transistor (JoFET). This is a three-terminal device; the supercurrent of the 
superconductor/semiconductor Josephson junction can be switched on and 
off by applying a gate voltage. Besides questions related to the JoFET, which 
might be used in future electronic circuits, some more fundamental questions 
related to the wave nature of the particles can also be addressed by study- 
ing superconductor/semiconductor structures controlled by a gate. Typical 
examples, which will be discussed below, are interference effects due to Fabry- 
Perot resonances, fluctuation and localization effects, and the quantization 
of the critical current. 



7.1 The Josephson Field Effect Transistor 

The possibility of controlling the electron concentration in a semiconductor by 
means of the field effect is one of the foundations of modern microelectronics. 
In a field effect transistor, the electron concentration and thus the resistance 
of the channel between the source and drain contacts are controlled by a 
gate electrode. Biasing the gate with a positive voltage leads to a higher 
electron concentration and, consequently, to a lower drain-source resistance. 
If a negative gate voltage is applied, the electron concentration is reduced, 
leading to an increase of the drain-source resistance. 

The concept of a field effect transistor can be transferred to superconduct- 
ing junctions. Here, the critical current in a superconductor/semiconductor 
junction is controlled by placing a gate electrode on top of the semiconductor 
channel, as depicted in Fig. 7.1a. This so-called Josephson field effect tran- 
sistor was first proposed by Clark et al. [7]. The idea behind the JoFET is 
the control of the critical current by changing the coherence length in 
the semiconductor link. As discussed in Sect. 6.2.5, and the critical cur- 
rent /c are increased with increasing electron concentration. A JoFET can 



Thomas Schapers: Superconductor/Semiconductor Junctions, STMP 174, 95—110 (2001) 
(c) Springer- Verlag Berlin Heidelberg 2001 



96 7. Gate-Controlled Superconductor/2DEG Junctions 

be employed to switch between two output voltages between the drain and 
source contacts. Let us define the “on” state as the state, where the link be- 
tween the source and drain is superconducting. The drain-source current is 
adjusted to a value below the critical current Ic of the junction (Fig. 7.1b). In 
the “ofP state a sufficiently negative gate voltage is applied that the electron 
concentration is reduced such that the critical current is suppressed. Now, a 
nonzero voltage V°^ appears between the drain and source contacts, which 
is determined by the load resistance i?L and the drain-source resistance . 




Fig. 7.1. (a) Schematic illustration of Josephson field effect transistor. The super- 
current flows through the channel between the source and drain contacts. The gate 
is separated by an insulator in order to suppress leakage currents, (b) Schematic 
current -voltage characteristics of a JoFET. In the “on” state the transistor is in 
the superconducting state. The dashed line represents the one state that occurs 
without superconductivity 



A few years after the proposal by Clark et al. [7] JoFET operation was 
reported by Nishino et al. [9, 178] and by Takayanagi and Kawakami [8]. The 
first device was based on an inversion layer formed by an Si 02 /Si structure. 
For the second device, the native inversion layer of InAs was used as a channel. 
Owing to the lack of a Schottky barrier for InAs, the gate electrode was 
insulated by an Si 02 layer from the channel in order to prevent gate leakage 
currents (Fig. 7.1a). In addition to these structures based on dielectric gate 
insulation, control of a Josephson supercurrent was gained in an InGaAs- 
based junction field effect transistor using a pn junction in order to insulate 
the gate from the channel layer [179]. 

The first JoFET with a two-dimensional electron gas in a semicon- 
ductor heterojimction was realized by employing a GaAs/AlGaAs layer 
system [10]. However, as outlined above, the large Schottky barrier sup- 
presses the superconducting coupling and thus reduces the device perfor- 
mance. In contrast, InAs-based heterostructures are the material of choice 
for a transparent coupling to a superconductor, owing to the inversion 
layer at the surface, and are consequently well suited as a basis for a 
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JoFET device. The first JoFET based on this type of heterostructure, 
realized with an lup 53 Gao. 47 As/InAs/Iuo. 53 Gao. 47 As 2DEG channel (Fig. 
3. jb), was fabricated by inserting an Si02 layer in order to insulate the A1 
gate electrode [180]. Unfortunately, the high density of interface states at 
the Si02/semiconductor boundary results in a relatively large gate voltage 
(> 10 V) being necessary to suppress the supercurrent between the drain 
and source contacts. Luckily, the Alin As top layer by itself provides a suf- 
ficiently high Schottky barrier. Therefore, the gate electrode can be placed 
directly on top of the AlInAs cap layer, providing that the gate electrode 
can be aligned between the source and drain contacts, which are only a few 
hundreds of nanometers apart. A cross section of this most advanced type of 
JoFET, fabricated by Akazaki et al. [11], is shown in Fig. 7.2. 



Gate 

AlInAs banier 
InGaAs 
InAs 
InGaAs 
AlInAs spacer 
AlInAs n+ 

AlInAs buffer 



/ InP s.i. substrate / 

Fig. 7.2. Cross section of Josephson field effect transistor based on an InAs 
inserted-channel InAlAs/InGaAs heterojunction. The gate length is typically 
100 nm, while the Nb electrode separation is a few hundreds of nanometers [11] 





Figure 7.3 shows the current-voltage characteristics of an InAs inserted- 
channel InAlAs/InGaAs heterojunction Schottky gate JoFET, as depicted 
in Fig. 7.2 [11]. At zero gate voltage a critical current of about 5 pA and a 
normal-state resistance of 16 Q are obtained for a 40 pm wide junction at 
a temperature of 1 K. Biasing the gate with a negative voltage results in a 
decrease of 1^ and, at the same time, an increase of i?N- A gate voltage of 
about 1 V is sufficient to suppress the Josephson supercurrent completely. 

The decrease of the critical current with increasing negative gate bias Vg, 
accompanied by a corresponding increase of i?N, is plotted in Fig. 7.4. Here, a 
steep increase of Ic with Vg is observed around a gate bias of —0.7 V, whereas 
Ic tends to saturate towards zero bias. 

In order to describe the dependence of the critical current on the gate volt- 
age, the expression derived by van Dover et al. [181] is frequently used. For 
a superconductor/normal-conductor/superconductor junction in the limit 
L, the critical current is given by 



Ic = R 



-1 

N 



7^'Aq gN^e ^^2 j2 R ^-L/£n 
2ekBTc 



(7.1) 
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Fig. 7.3. Current -voltage characteristics of an InAs inserted-channel In- 
AlAs/InGaAs heterojunction JoFET for successively increasing gate bias from 0 
to -1.1 V [11] 




(V) 



< 

d. 



Fig. 7.4. Critical current 7c and normal resistance 7 ?n as a function of applied 
gate bias at T = 1 K. The slope of the 7c vs. Fg curve is the superconducting 
transconductance [11] 



Here, gs and pN are the resistivities in the superconductor and normal con- 
ductor, respectively, and / is the ratio of the order parameter at the interface 
to its bulk value. The factor by which the order parameter changes on cross- 
ing the interface is described by A. Strictly speaking, (7.1) is only valid near 
Tc since it is derived from the Ginzburg-Landau theory; however, in most 
practical cases it also applies at lower temperatures [182]. If no Schottky 
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barrier is present at the interface, as for InAs, is given by gsm'^/ g^nie, 
which simplifies (7.1) [182]. However, if a Schottky barrier has to be taken 
into account, as for Si or GaAs, A is proportional to the tunneling probability 
through the interface barrier [182]. 

In (7.1) i?N and depend on the carrier concentration ri 2 D- For the 
dirty-limit case, the coherence length is given by 



yd _ 
SN — 



hv 



N 






2Trk^T 



2m*ek'BT 



V^2D , 



(7.2) 



while the coherence length in the clean limit can be expressed as 

fiV-p 

^ 2i:kpT ^ ' 



(7.3) 



In both cases the coherence length is proportional to ^/n^. Since ri 2 D cx Ig — 
Vt, where Vt is the threshold voltage, the coherence length is proportional to 
y/Vg — Vt. Kleinsasser et al. [182] analyzed the gate voltage dependence of Ic 
of an Si/Si02 JoFET fabricated by Nishino et al. [9]. A very good agreement 
with the experimental data was obtained by using (7.1), with A determined 
by the Schottky barrier transparency. For a JoFET with an InAs inserted 
channel, a transition between the clean and the dirty limit was observed 
when the carrier concentration in the channel was decreased by applying a 
negative gate voltage [180]. As mentioned above, (7.1) is strictly valid only in 
a limited temperature range close to Tc. For long contacts (L > ^n), Kresin 
[147] developed a model for the critical current as a function of ri 2 D in the 
clean and in the dirty limit. Within this model, the different effective masses 
and the potential step between the superconductor and the semiconductor 
were considered, whereas the barrier at the interface was omitted. For the 
clean-limit case, a barrier was taken into account in the model presented by 
Chrestin et al. (see Sect. 6.1.3) [148]. Volkov and Takayanagi [183] analyzed 
the dependence of Ic on ri 2 D for a JoFET in the dirty limit, where a gate 
length shorter than the channel length was taken into account. A comparison 
with experimental data on an InAs-based JoFET [184] confirmed the validity 
of the model. 

One parameter which characterizes the device performance of a JoFET is 
the superconducting transconductance g^, defined as [11] 



= dJc/dVg , 



(7.4) 



where Jc is the critical current for a width of 1 mm. From the Ic vs. Vg 
curve in Fig. 7.4, a maximum of 350 |iS/mm is found at a gate bias of 
—0.7 V, which is considerably larger than the g®, values obtained for JoFET 
structures with an Si 02 insulation layer [11]. 

With regard to voltage gain, one of the crucial parameters concerning any 
possible application in electronic circuits, the JoFET has some important 
limitations. This has to do with the fact that a gate bias of the order of at 
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least one volt has to be applied, while the output voltage between the source 
and drain contacts is of the order of A^je. Of course the output voltage swing 
can be increased, but then the superconducting properties of the JoFET are 
not taken advantage of, so that the device is operated more as a conventional 
field effect transistor [179]. The voltage gain is defined as the voltage 
ratio between the output voltage and the gate voltage swing. For an inverter 
circuit with a load resistance i?L, the voltage gain is given by [185] 

I I (i?L + I r 

Here I™ and refer to the superconducting “on” state, while R^ and 
Vg°® are the normal resistance and the gate voltage in the resistive state 
(“ofF state) of the junction. Owing to the superconducting “on” state of 
the JoFET, is zero (Fig. 7.1). For the transistor shown in Fig. 7.2, a 
voltage gain of the order of one could be obtained since relatively small gate 
voltages (E ~ 1 V) were sufficient to deplete the drain-source channel [11]. 
Nevertheless, a relatively large load resistance i?L had to be chosen in order 
to obtain G > 1. This has a negative effect on the maximum operation 
frequency, which is determined by i?L and the gate capacitance Gg. For a 
real application of a JoFET in an electronic circuit, one of the major tasks 
is to increase by increasing Ic [11]. 



7.2 Interference Effects 

Besides its potential use as an electronic device, a JoFET structure is also 
very interesting for studying interference effects. In the clean-limit regime 
Fabry-Perot-like interference patterns can be found in the critical current, 
while in the dirty limit localization effects and mesoscopic fluctuations are 
observed in the critical current. 

7.2.1 Critical-Current Oscillations 

If the superconductor/2DEG interfaces of a JoFET are nonideal so that nor- 
mal reflection can occur, a JoFET can serve as a resonator similar to a 
Fabry-Perot resonator in optics. Fabry-Perot-like oscillations in a Josephson 
junction were first discussed by Gudkov et al. [12] in conjunction with the 
interpretation of three-dimensional Nb/a-Si junctions. For a junction with a 
two-dimensional electron gas as a normal conductor, oscillations due to inter- 
ference effects were discussed by Ghrestin et al. [148], while critical-current 
oscillations in a one-dimensional junction were investigated theoretically by 
Rittenhouse and Graybeal [13] and by Tang et al. [144]. 

In order to discuss the origin of the interference effects in a supercon- 
ductor/2DEG junction, (6.27) and (6.28), derived by Ghrestin et al. [148], 
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can be referred to. Here, the first two terms in f(ky,u>n) can be attributed 
to an electron Andreev-reflected as a hole at the superconductor/2DEG in- 
terface and vice versa. The remaining two terms in (6.28), leading to the 
critical-current oscillations, arise when normal reflection is involved in ad- 
dition. One of the processes contributing to the oscillations is sketched in 
Fig. 7.5, where two Andreev reflections and two normal reflections occur. 
If the carrier concentration is changed, thus changing the electron and hole 
wave vectors, resonant states within the S/2DEG/S cavity arise, leading to 
a maximum in the critical current. 





e 




AE 




h 


)ar 




c 


h 






\NR 


h 


)nr 


S 


2DEG 


S 



Fig. 7.5. One of the processes contributing to the oscillations of the critical current. 
This process contains two Andreev reflections (AR) and two normal reflections (NR) 



The oscillations of the critical current as a function of the carrier concen- 
tration, calculated from (6.27), are shown in Fig. 7.6 for a 2DEG in an InAs 
inversion layer contacted to Nb electrodes. Two features can be made out. 
First, the oscillation amplitude increases with increasing barrier strength Z 
at the Nb/InAs interface. For Z = 0 only weak oscillations are found, while 
for finite Z parameters more pronounced oscillations occur owing to the in- 
creased normal-reflection probability. Second, for Z = 0 the oscillations al- 
most disappear at an electron concentration of n « 1.2 x 10^^ cm“^. This 
is due to the fact that at this particular electron concentration, the Fermi 
velocities of the superconductor and the 2DEG are equal, so that even the 
usual normal reflection due to the Fermi velocity mismatch is suppressed. 

The local maxima of the critical current as a function of the carrier density 
ri 2 D can be approximated by [148] 

?T-2D,max = 

where L is the length of the junction and the integer N equals half the number 
of quasiparticle states with resonant values of ky, defined by the maxima of 

cos (^2^kp^ — k'^L^; ky is the wave number perpendicular to the current 
direction. 
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Fig. 7.6. Critical current as a function of the carrier concentration Us for various 
barrier strengths Z [148] 



Initial experiments showing critical-current oscillations as a function of 
the carrier concentration were reported by Takayanagi et al. [27]. However, 
in these JoFET structures the Nb electrodes were deposited on top of the 
InAs channel of the heterostructure shown in Fig. 2, !a, so that the Nb/2DEG 
interface was not defined sharply. Consequently, the oscillation amplitude of 
Ic was relatively small. In addition, the junction length L deduced from (7.6) 
did not fit well to the geometrical Nb electrode separation. A much stronger 
oscillation amplitude was found for the improved JoFET structure shown 
in Fig. 7.2 [186]. Here, the Nb electrodes are contacted to the 2DEG at the 
mesa edge, so that the geometry of the resonator is well defined. The measured 
oscillations of Ic as a function of the carrier concentration underneath the gate 
electrode are shown in Fig. 7.7. Since the gate does not cover the complete 
2DEG area between the Nb electrodes, the electron and hole wave vectors are 
only altered in the area underneath the gate. The Ic peak positions calculated 
taking account of this effect, marked by arrows, show a reasonable agreement 
with the experimental values. 

7.2.2 Mesoscopic Fluctuations and Localization Effects 

Localization effects and universal conductance fluctuations (UGFs) are two 
remarkable phenomena in mesoscopic physics, which directly demonstrate 
the wave nature of the electrons in the conductor [149]. In addition to such 
effects in purely metallic or semiconducting systems, both effects have been 
predicted and experimentally observed in the Josephson supercurrent of su- 
perconductor/2DEG junctions. 

With regard to universal conductance fluctuations, Al’tshuler [187] and 
Lee and Stone [188] showed, for a disordered wire, that the mesoscopic flue- 
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Fig. 7.7. Critical current as a function of the carrier concentration underneath the 
gate electrode. The arrows indicate the calculated peak positions [186] 



tuations in the conductance G at zero temperature have a root mean square 
value of rmsG = yj (G^) — (G)^ ~ e^jh. Since this value is independent 
of the wire length and the elastic mean free path, these fluctuations are 
called “universal” . Universal conductance fluctuations manifest themselves 
as sample-to-sample fluctuations. However, owing to the ergodic property, 
changing the Fermi Energy E-p or the magnetic field B is equivalent to av- 
eraging G over an ensemble of samples. For a given sample, the fluctuations 
as a function of Ep or B are sample-specific and also reproducible, provided 
that the sample is kept at low temperature. One can say that each sample 
possesses its own “fingerprint” . 

From the theoretical point of view, fluctuations in superconductor/2DEG 
Josephson junctions have been treated in two limits. Al’tshuler and Spivak 

[189] discussed a junction which was long compared with the elastic mean 

free path and with the superconducting coherence length in the dirty limit 
(i 3> ^ei)- Their result for the fluctuations in the critical current rmsic = 

(/^) — (/c)^, expressed for a two-dimensional system, can be written as 

[190] 

eV fW 

rms/c«2.2^^^ (7.7) 

for L <^W and T = KD/2kp,L^. In order for one to be able to observe 

these mesoscopic fluctuations, the correlation energy E^ = Titt'^'D/L'^ must 
be larger than kpT so that thermal smearing is prevented. As can be seen 
in (7.7), the fluctuations in Ic depend on the sample length L as well as on 
the elastic mean free path Zei, since V depends on Consequently, these 
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fluctuations cannot be regarded as “universal” in the spirit of the universal 
conductance fluctuations. 

Universal fluctuations of the critical current were predicted by Beenakker 
[191] for the case when the junction length is much shorter than the super- 
conducting coherence length. In this short-junction regime, 

rmsic ~ eA^jh (7.8) 

for T <?; T(,. This result was obtained by relating the Josephson supercurrent 
to the scattering matrix of the normal region of an SNS junction. As for the 
theory of UCFs, the statistical properties of the scattering matrix were used 
to obtain (7.8). 

In order for one to observe fluctuations in A,, the transport in the junction 
must be diffusive so that l^i ^ L is fulfilled. Consequently, the first critical- 
current fluctuations were found in a junction based on an inversion-layer two- 
dimensional electron gas in p-InAs [190], where the ionized impurities and 
the rough interface lead to a sufficiently large number of elastic scattering 
events. For the superconducting electrodes, Nb was employed. 

With regard to UCFs, an external magnetic held is often used in or- 
der to trace the fluctuations in G. Since in spatially extended superconduc- 
tor/semiconductor junctions the supercurrent is very sensitive to B, it is in 
this case more appropriate to change the Fermi energy E-p by using a gate 
electrode on top of the semiconductor bridge. As can be seen in Fig. 7.8, 
the critical current fluctuates if the gate voltage is changed. Furthermore, Ic 
follows precisely the fluctuations of the differential normal conductance Gn. 
A detailed analysis of the experimental data showed that the sample satisfied 
the conditions of the theory of Al’tshuler and Spivak [189], since and lei 
were much smaller than the separation of the superconducting electrodes. 

In the diffusive transport regime the elastic mean free path lei is much 
shorter than the sample length, and thus the electrons undergo many scat- 
tering events while passing through the sample. At some locations multiple 
scattering leads to the situation that the electrons are moving in closed loops. 
In these closed loops the electrons can move in two opposite, time-reversed 
directions. Since the phase accumulations of the electron waves are identical 
in both directions, a constructive interference occurs, leading to a localization 
of the electron within the loop. In transport measurements this localization 
results in an additional contribution to the resistance. If a magnetic held is 
applied, the time-reversed paths acquire opposite additional phase shifts due 
to the vector potential A. This destroys the constructive interference and 
reduces the sample resistance. The same effect is obtained by increasing the 
temperature. 

Similarly to the decrease of the normal conductance of a metal due to 
localization, a reduction of the critical current in a superconductor/dirty- 
two-dimensional-metal Josephson junction was predicted by Fukuyama and 
Maekawa [192, 193]. In their theory the diffusion constant T> is modified by 
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Fig. 7.8. (a) Critical current and differential normal conductance Gn as a 
function of the gate voltage, (b) Fluctuations of Gn as a function of B for a fixed 
gate voltage [190] (© (1995) by the American Physical Society) 



localization. As can be seen in (7.2), this directly affects the coherence length 

and the critical current Ic (7.1). 

Verifying localization effects in superconducting junctions faces the same 
problem as does the detection of fluctuations in Ic, a magnetic field cannot 
be employed, owing to the strong dependence of Ic on B. Therefore, Ic was 
measured as a function of temperature (Fig. 7.9) [184]. Here, the inversion 
layer in p-type InAs was used as the normal-conducting bridge between the 
superconducting Nb electrodes. By applying an appropriate gate voltage, the 
2DEG channel can be adjusted to the weak (higher electron concentration) 
or strong (lower electron concentration) localization regime. The localization 
manifests itself in a saturation of Ic at low temperatures. 

As can be seen in Fig. 6.11 for a junction in the clean limit, a saturation 
of Ic at low temperature does not always mean that localization effects occur. 
For junction with L and a constant i?Nj the critical current usually 

shows a saturation at T < 0.2Tc [1]. It was shown that for the junction used 
for the measurements depicted in Fig. 7.9, this condition was not fulfilled. 
This proves that the saturation must be due to localization. 



7.3 Superconducting Qnantnm Point Contacts 

In a two-dimensional electron gas, the free motion of the carriers is restricted 
to a plane. One can proceed one step further and allow transport in only 
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Fig. 7.9. (a) Critical current of an Nb/p-type InAs/Nb junction as a function of 
temperature. The solid and dotted lines show the calculated 7c vs. T curves with 
and without localization effects. In (b) the temperature dependence of a reference 
junction in the clean limit with an InAs inserted channel (see Fig. 'b) is shown 

for comparison [184] (© (1995) by the American Physical Society) 



one direction. This can be obtained by using a split-gate geometry, where 
a small opening is defined by two gates, separated by a few hundreds of 
nanometers. If a sufficiently large negative gate voltage is applied, the two- 
dimensional electron gas underneath the gate is depleted, leaving only a small 
opening for the electrons to pass through. Since the Fermi wavelength Ap 
of a 2DEG in a semiconductor heterostructure is comparable to the width 
of the point contact defined by the split gate, the transverse momentum 
is quantized. In ballistic structures, where the channel is shorter than 
quantized conductance in units of 2e^ /h can be observed [194, 19-5]. 

The following question arises: is it possible to observe a quantization of 
the supercurrent if two superconducting electrodes are connected by a point 
contact? This problem was addressed theoretically by Furusaki et al. [2.5, 196, 
197] and by Beenakker and van Houten [’’j;, 198]. Here, the latter approach, 
which applies to a short adiabatic constriction, will be discussed briefly. 

The geometry of the constriction is shown in Fig. 7.10. It widens from a 
minimum width of 5’min to 5’max, with Smax at both sides. The length 

L is assumed to be much shorter than the coherence length = Tt-vf/ttAq. 
The number of modes propagating through the constriction is given by 
N ~ S’lnin/Ap. Adiabatic propagation, and thus no scattering between the 
modes, is assumed, which requires L ^ Ap. As depicted in Fig. 7.10, the pair 
potential at both ends of the constriction is approximated by Aq exp(— i</>o/2) 
for X < — T/2 and by Z\o exp(i</fo/2) for x > L/2. The Bogoliubov equa- 
tion (2.1) for the constriction is solved by expanding the eigenfunction 
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Fig. 7.10. Schematic illustration of a constriction between two superconducting 
reservoirs (after [24]) 
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into transverse modes ipn{x, y), which are the eigenfunctions of the Hamilto- 
nian — {h^/2me)d‘^/dy'^ + U{x, y). The corresponding eigenvalue for the mode 
n at the position x is En{x). In order find the solutions for the propagation 
along the constriction, the WKB method [199] is used. Owing to this approach 
only the discrete solutions jifl < Aq contribute to the Josephson current [24]. 
In the short-junction limit (L ^ 0 ) the normalized wave function is given 
by 

{Aq- 

X exp ^±i J k{x')dx''^ , (7.10) 

where ko{x) = {2me[/r — £'„(x)]/?i^}^/^. For Z\q ^ fi — En{x), the wave vector 
k{x) in the integral can be approximated by 

t(x) = (1 - , (7-11) 

where 4>{x) is the superconducting phase within the junction. In contrast to 
the one-dimensional case discussed in Sect. 6.1.1, the wave vector does change 
along the constriction owing to the quantization in the transverse direction. 
Since ip±{x) must decrease exponentially for jxj < T/2, the wave function 
ip+{x) corresponds to 0 < <(>0 < tt, whereas ip-{x) is appropriate to the 
solution for — 7T < <()o < 0. In the limit T/^o ^ 0, the Andreev bound-state 
energy is given by 

arccos(A/Ao) = ±^<('0 , (7.12) 

which corresponds to the solution of a one-dimensional junction in the short 
limit (6.7). It should be stressed that, for a constriction in the short limit 
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as discussed here, the Andreev bound-state energy is identical for all propa- 
gating modes n. From the normalized wave function (7.10), the supercurrent 
per mode can be calculated, which is 

In = ^Aosm(j)o/2 (7.13) 



at T =0. Since all modes contribute by the same amount to the supercurrent 
the critical current for N occupied modes is given by 



Ic 




(7.14) 



If the number of channels N is increased, the critical current increases step- 
wise. The steps have a constant height of eA^/h, which is determined solely 
by fundamental constants and the bulk pair potential Aq. 

The result that the critical current of a superconducting point contact is 
independent of the parameters of the junction holds only for a short junction 
in the adiabatic approximation. Furusaki et al. [25] treated the more general 
case beyond this limitation for various geometries of the constriction. 

For a long constriction in the adiabatic approximation, L 3> Furusaki 
et al. arrived at the following expression for the critical current at zero tem- 
perature: 







(7.15) 



Here the velocity of the nth mode is defined as 




where = {2mefJ./h^ — [mi / . S{x) is the width of the constriction 
at position x. In contrast to a constriction in the short limit, the critical 
current for L ^ depends on the junction parameters such as L and 
For the simplest adiabatic constriction with a constant width S{x) = 5o, the 
velocity of the nth mode is given by Vn = [2/r/me — {mrfi/ ■ This 
implies that an increase of the width Sq of the constriction would lead to a 
smooth increase of Thus no sharp steps, as in the short limit, are expected. 
In addition, the step height depends on the width of the constriction. In 
contrast to the adiabatic approximation, where the junction length is much 
larger than Ap, in constrictions with sudden changes of width, mode mixing 
and reflection can occur. As a result the critical current as a function of the 
constriction width can show additional features, i.e. resonant structures [25]. 

A superconducting quantum point contact can be realized by placing a 
split gate between two closely spaced superconducting electrodes (Fig. 7.11). 
In order to operate the structure in the ballistic regime, a high- mobility two- 
dimensional electron gas in a semiconductor heterostructure has to be chosen. 
If a sufficiently large negative bias voltage is applied to the split gate, the 
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Fig. 7.11. Schematic drawing of a split-gate superconducting point contact 



electron gas underneath the gate can be depleted. The opening width of the 
point contact can be adjusted by the gate voltage. 

The quantization of the critical current was demonstrated by Takayanagi 
et al. [2G]. Figure 7.12 shows the conductance and Ic as a function of the volt- 
age applied to the split gate. The conductance G of the constriction measured 
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Fig. 7.12. Conductance of the constriction and the critical current as a function 
of the gate voltage [26] (© (1995) by the American Physical Society) 



at a small magnetic field shows clear steps. From theory it is predicted that 
for an ideal point contact, facing two superconducting interfaces, the conduc- 
tance should be quantized in units of 4e^//i owing to Andreev reflection [200]. 
Although a quantization in steps of 4e^/h was not obtained, the conductance 
steps are larger than 2e^ jh. The latter step size is expected for a completely 
normal-conducting quantum point contact [2G]. It can thus be concluded, at 
least, that Andreev reflection processes are involved in the carrier transport 
through the constriction. 
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The steps in the critical current Ic follow the steps of the conductance 
G, as can be seen in Fig. 7.12. The step height in 1^ is considerably smaller 
than predicted by the theory of Furusaki et al. [25, 196]. This might be due 
to the fact that in the theoretical model no normal reflection was included. 
As we already know from the Josephson junctions discussed in Sect. 6.1.2, 
a finite barrier at the interface or a mismatch of the Fermi velocities can 
lead to a decrease of the critical current. One discrepancy, which has not 
been explained up to now, is the reduced number of steps in compared 
with the number of steps in G (Fig. 7.12). This result implies that some one- 
dimensional channels do not contribute to the superconducting transport. 



8. Nonequilibrium Josephson Current 



Control of a supercurrent in weak- link Josephson junctions can be realized 
by two alternative methods. The first approach, which was described in the 
preceding chapter, is to control the carrier concentration in the semicon- 
ductor with a gate electrode. A typical example of this kind of structure is 
the Josephson field effect transistor. In this chapter we shall discuss how a 
Josephson supercurrent can be controlled by injecting nonequilibrium carri- 
ers. As explained in Sect. 6.1, the net supercurrent is partially determined 
by the occupation of Andreev bound states located energetically within the 
superconductive gap. If carriers are injected, the occupation of these Andreev 
levels is altered. This directly affects the net supercurrent in the junction. In 
the first part of this chapter the underlying theoretical models will be dis- 
cussed in detail, while in the second part the corresponding experiments are 
described. 



8.1 Supercurrent Control by Nonequilibrium Carriers: 
Theoretical Concepts 

The concept of coupling a normal reservoir to a ballistic SNS junction in 
order to control the Josephson supercurrent by applying a bias voltage to the 
normal electrode was first proposed by van Wees et al. [15]. These authors 
could show that the coupling to the normal reservoir has two effects. First, the 
Andreev resonances are broadened owing to phase-breaking processes. This 
leads to a decrease of the critical current even for zero bias voltage. Second, 
if an appropriate voltage is applied to the normal reservoir, the occupation 
of the Andreev levels is altered, which directly affects the supercurrent in the 
junction. Since the current direction changes sign for energetically neighbor- 
ing Andreev bound states, even a reversal of the supercurrent was predicted 
for certain ranges of bias voltages. The injection of hot carriers into ballistic 
junctions of various lengths was investigated by Chang et al. [17]. By taking 
the scattering phases at the connection point between the junction and the 
injector into account, Samuelsson et al. [18, 201, 202] were able to show that 
for a one-dimensional structure, an additional nonequilibrium current contri- 
bution has to be considered. The mechanisms described here can in principle 
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be used for a superconducting device where the supercurrent is controlled 
by an injection current. From a fundamental point of view, a three-terminal 
structure could even be useful as a spectrometer in order to resolve single An- 
dreev levels [202, 203]. The finite phase difference between the two electrodes 
of a Josephson junction under supercurrent flow was studied theoretically 
by Nakano and Takayanagi [20-l] in an interferometer setup. Here, the phase 
difference was employed to control the conductance of a normal-conducting 
Y-shaped beam splitter, which was attached to both superconducting elec- 
trodes. 

For the diffusive case, Volkov [IG] was able to show that in this regime a 
reversal of the supercurrent can also be expected. Furthermore, the results 
of Wilhelm et al. [205] and Yip [19] have demonstrated that the spectral 
supercurrent i{E) oscillates as a function of energy. This directly implies 
that the supercurrent can be reversed if a certain bias voltage is exceeded. 

8.1.1 Carrier Injection in a Three- Terminal Structure 

In order to discuss the effect of carrier injection on the Josephson supercur- 
rent in detail, we follow the approach of Samuelsson et al. [18, 202]. Here, 
a three-terminal junction, as depicted in Fig. 8.1, is considered, where the 
two superconducting electrodes are connected by a superconducting ring. By 
threading a magnetic field through the loop, the phase difference (j) can be 
adjusted independently. Additional carriers are injected from a normal reser- 
voir biased at a voltage V. The superconducting loop is at ground potential. 

superconducting 



Fig. 8.1. Schematic illustration of a three-terminal junction, where the phase dif- 
ference can be controlled by a magnetic field threading the superconducting loop. 
The normal conductor is connected to the grounded loop via a voltage source biased 
at V. The sketch on the right shows a detail of the junction area with the directions 
of the current flow in each lead (after [202], © (2000) by the American Physical 
Society) 
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The transport is considered to be one-dimensional, and the normal conductor 
is connected halfway between the superconducting electrodes (see detail in 
Fig. 8.1). The ingoing and outgoing waves ifin and iFout at the connection 
point are related by a scattering matrix S [15, 18, 202]; 

tf'out = s , (8.1) 

where 

/ VI - 2£r a/s \ 

S = Ve r d \ . (8.2) 

\ Ve d r J 

Here, r and d are the reflection and transmission amplitudes for the scattering 
between leads 2 and 3; we define D = \d\^ and R = jrp. (0 < e < 0.5) 
is the scattering amplitude from the injection lead 1 to lead 2 or 3, where 
D+R = 1 — e and Re(rd*) = e/2. Moreover, Im(r<i*) = RD — e^/4, where 
(T = ±1 depending on the phase of the scatterer. Maximum coupling of the 
normal reservoir is achieved for e = 0.5, while the reservoir is completely 
decoupled for e = 0. Normal reflection at the SN interface is neglected. 

When the dwell time of the Andreev resonances in the junction is smaller 
than the inelastic scattering time, the quasiparticle distribution in the nor- 
mal region is determined by the Fermi distribution function of the normal 
reservoir. Under this condition, the currents in the leads j = 2 and 3 resulting 
from the injected particles can be written as 

/ -l-oo 

dU(VVe + + i^k) , (8.3) 

-OO 

where are the current densities of injected electrons and holes and /e,h = 
fo{E ±eV) are the corresponding Fermi distribution functions in the normal 
reservoir. The last term in (8.3) is the contribution of the quasiparticles 
injected from the superconductor. 

In order to obtain the currents in the different leads, the scattering states 
in the three normal regions have to be calculated by employing the Landauer- 
Biittiker scattering approach [206, 207]. The plane waves in the different 
regions of the normal conductor, as well as in the superconducting electrodes, 
are similar to (2.10) and (2.11) in Sect. 2.2.2, but now we must take both 
directions of wave propagation into account. The matrix S, as defined in 
(8.2), is used to describe the scattering of the plane waves at the junction 
point. The resulting current densities in leads 2 and 3 are given by [202] 

= — ^^{2D sin 4> sin 20 ± (4 ct \/ RD — j A sin sin^ 0) 

+ 2esin^ 0(1 -I- cos(/)} , (8.4a) 

^ — {2Z?sin 0sin 29 ± (4ct \/ RD — e^/qsin i/sin^ 0) 

— 2£sin^ 0(1 -I- cos(/)} , 



(8.4b) 
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where the superscripts e and h are related to the + and - signs, respectively, 
and where 0 is defined as 

9 = arccos(i5/zio) — EL/{Tivp) . (8.5) 

r is given by 

/^ = [(1 — e) cos 20 — R — Dcostjif' + sin^ 29 . (8.6) 

In order to separate the different contributions resulting from the coupling 
of the normal conductor, the currents in leads 2 and 3, as given in (8.3), can 
be rewritten as [202] 

(/e - /h) + , (8.7) 

where and i~ = ~ In equilibrium (F = 0) the current 

through the junction is simply given by 

= JdE{i+ + i^)fo, (8.8) 

where 

= i 2 = = — ^ — D sin cj) sin 29 . (8.9) 

The resulting current density for a short junction is shown in Fig. 8.2. Owing 
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Fig. 8.2. Charge current density i^/{eh) as a function of the energy EjAo for a 
junction in the short limit (L — 0). The coupling strength e is varied as a param- 
eter from 0.1 to 0.3. The transmission and reflection probabilities were calculated 
according to D = (Vl — 2e — 1)^/4 and R = {y'l — 2e + 1)^/4 
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to the coupling to the normal reservoir (e > 0), the charge current density 
is now extended over a broad energy range. This is in contrast to the case 
where the junction is decoupled from the reservoir. For the latter case, only 
a sharp current density peak at 



= TZio^l - T»sin2(,/)/2) (8.10) 

is found. The broadening of the current density spectrum for finite values 
of e results in a rounding of the current-phase relation /(</>) and thus to a 
decrease of the critical current [15, 17]. 

Part of the nonequilibrium current (V yf 0) can be associated with the 
nonequilibrium population of the resonant Andreev states that exist: 

A = JdE (^^(/e + /h - 2fo)j . (8.11) 

The current contributions Aq and A are both related to the current density 
which is equal in both leads. However, the difference current densities 
ij are not equal in leads 2 and 3, as can be seen by comparing (8.4a) and 
(8.4b). Nevertheless, the second terms of the two expressions are identical 
and, in addition, this contribution survives even in the limit of e — > 0. The 
corresponding current density contribution of i~ = defined as the 

anomalous current density Za, is given by 

Za = — RD — £2/4sin(()sin^ 9^ , (8.12) 

resulting in an anomalous current 

7a = |d£;(|(/e-/h)) . (8.13) 



The anomalous current is due to the important modifications of the Andreev 
states caused by the open normal lead. This situation is depicted in Fig. 8.3. 
Owing to the coupling to the normal reservoir, the electron and hole current 
densities z® and z^ differ. This leads to a finite difference current z“, which 
has the same sign for both resonances. 

The contributions to Z 2 and , which are different in leads 2 and 3, 
originate from the current density Zinj injected into lead 1: 



4ee^ . 

Zinj = Z 3 — Z 2 = — — sin 6^(1 + cos cp) . 
The corresponding injection current in lead 1 is 



h = /inj = dE\ ^(/e - /h) 



(8.14) 



(8.15) 



For a short junction (L = 0), only two resonant Andreev states exist for 
\E\ < Z\o and the continuous current is zero (see, e.g.. Fig. 6.5). At zero 
temperature and for e 1 only the lower one of these levels is populated. 
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Fig. 8.3. The charge current density for two resonant Andreev levels for injected 
electrons {dotted line) and holes {dashed line)', {solid line) and i~ {dash- 
dotted line) are the sum and difference current densities, respectively. The difference 
current has the same sign for both resonances. The inset shows the corresponding 
current paths for an injected electron {solid line) and a hole {dashed line) at reso- 
nance [202] (© (2000) by the American Physical Society) 



If a voltage is applied {V > 0), the electron and hole populations are shifted 
upwards and downwards, respectively, in energy. For a voltage that reaches 
the upper Andreev level ), this level becomes populated and the current 
in the junction jumps abruptly, as shown in Fig. 8.4. The contribution of Igq 
is totally compensated by the regular current A, so that Jeq -I- A is zero for 
eV > Ag . This cancellation can be understood as a result of the occupation 
of two resonant Andreev levels with opposite current directions (see Fig. 8.2). 
The anomalous current /a jumps by the amount [202] 



eAo V^I?sin(^| sin((^/2)| 
^ 2fi 1 — Dsin(^/2) 



(8.16) 



and is thus the only remaining current for eV > Ag (Fig. 8.4). 

In a long junction with A 3> ^Oj several resonant Andreev levels are found 
within the gap (|£’| < Aq), as can be seen, for example, in Fig. 6.2. Each of 
these Andreev levels carries a supercurrent in the opposite direction to its 
neighboring states. As a consequence, most of the current is compensated in 
equilibrium. If an increasing voltage is applied to the normal lead the current 
leq + A reverses sign each time a further Andreev level is reached by eV, as 
shown in Fig. 8.5. In contrast to the current reversal due to the oscillating 
contribution of A, the anomalous currents /a add up with the same sign for 
all levels, as can be seen in Fig. 8.6. The corresponding current-voltage char- 
acteristic of the anomalous current contribution has a staircase-like shape. 
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Fig. 8.4. The sum of the equilibrium and regular currents /eq + 7r {dash-dotted 
line) and the anomalous current 7^ {dashed line) as a function of the applied voltage 
eV/Ao for a short junction with L = 0, for D = 0.8, (j) = 37t/4, e = 0.01 and a = — 1 
at T = 0. 7eq + It jumps to zero if the applied voltage reaches the upper Andreev 
level. The remaining total current 7eq+7r+7a {solid line) is in this regime exclusively 
carried by the anomalous current [202] (© (2000) by the American Physical Society) 




In principle, this anomalous nonequilibrium Josephson current is expected to 
be observable even when the equilibrium supercurrent or the regular current 
is exponentially small at finite temperature. However, since the anomalous 
current is a quantum interference term, this contribution depends very sensi- 
tive on the scattering phases. Thus, for junctions which are extended in two 
dimensions, the scattering phases for different modes will be different and 
on average will roughly be zero [208]. Under these conditions. It would be 
the only observable nonequilibrium current contribution at sufficiently low 
temperatures. 

8.1.2 Carrier Injection in a Four- Terminal Junction 

In the treatment above, a three-terminal configuration was assumed. A four- 
terminal junction, with two normal electrodes on each side of the junction, 
has the advantage that a symmetric injection condition can be obtained by 
biasing one normal electrode at U/2 and the opposite one at — U/2. The ex- 
tension to a four-terminal structure in the fully ballistic regime was discussed 
theoretically by Than and Bagwell [209] and by Samuelsson et al. [202]. Here, 
we shall consider an alternative concept of a four-terminal junction, where the 
normal-conducting part of the superconductor is found in the center part of a 
diffusive normal-conducting wire. This configuration was chosen by Wilhelm 
et al. [205] in order to model a diffusive junction with two long diffusive wires 
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Fig. 8.5. The equilibrium current Jeq plus the regular nonequilibrium current T as 
a function of the applied voltage cV/Aq at T = 0 (soZid line), ksT = 0.04Aq (dash- 
dotted line) and ksT = 0.07Ao (dashed line). The following parameters were used: 
L = 10^0, (/) = 7 t/ 2, D = 0.8 and e = 0.05. For eV > Aq the current corresponds 
to the continuous current [202] (© (2000) by the American Physical Society) 

connected to normal reservoirs. Here, we shall discuss the case of diffusive 
normal-conducting wires attached to a junction in the ballistic regime [210]. 
This condition is fulfilled if the length Tw/2 of each piece of wire is much 
longer than the elastic mean free path but the junction length L, which is 
equal to the wire width, is smaller than lei- The junction itself is regarded as 
one-dimensional . 

If we assume that the two ends of the normal reservoirs are connected by 
a sufficiently long wire, with a total length smaller than the inelastic scat- 
tering length lin, then, according to Pothier et al. [211, 212], the distribution 
function along the wire can be expressed as 



f(E) = (1 - y/L^)fo{E - eV/2) + (y/L^)fo(E + eV/2) . (8.17) 



Here, fo(E ± eV/2) are the distribution functions of the normal reservoirs 
shifted by ±eV/2 with respect to the chemical potential of the supercon- 
ductive electrodes. Let us assume a short one-dimensional junction located 
in the middle of the wire (y = Lw/2), then the distribution function at the 
position of the junction has a double-step shape, as depicted in Fig. 8.7. The 
supercurrent in such a short junction, with d barriers at the SN interfaces, 
can be obtained by inserting the modified occupation of E^ into (6.26) [210]: 



7dis(^/^) 



eA 




^ 4y^[cos^((j)/2) + 4Z2](1 -b 4Z2) 




X 



(8.18) 
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Fig. 8.6. The anomalous current as a function of the applied voltage eV/Ao for 
{curve a) (fi = tt/ 4 and {curve b) (f> = 37t/4 for a junction with the following 
parameters: L = lO^o, 4' ~ '^/2, D = 0.8 and e = 0.05. At zero temperature {solid 
line) steps are found, which are averaged out if T is increased to 0.1 Ao/ltB {dashed 
line). The upper inset shows the critical anomalous current at eV = Aq as a function 
of the transparency D for e = 0.1. The lower inset shows the anomalous current 
Ia{eV — Ao,kB S> hv-p/L) as a function of the phase difference for D = 0.1, 0.5, 
and 0.9 [202] (© (2000) by the American Physical Society) 




Fig. 8.7. Decrease of the critical current due to nonequilibrium carriers. The partial 
filling of the Andreev levels for energies — eF/2 < E < eV j2 leads to a cancellation 
of the supercurrents carried by the two Andreev levels. As a consequence the critical 
current is reduced from 7co to E 



The effect of a nonequilibrium distribution function on a short junction 
at T = 0 with no barriers {Z = 0) is depicted in Fig. 8.7. Owing to the 
voltage V applied between the ends of the normal-conducting wire, the An- 
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dreev levels in the energy interval [— ey/2,eV/2] are only halffilled because 
of the modified distribution function f{E). As a consequence the supercur- 
rent Idis{4>) switches to zero when the phase difference (j) lies in the interval 
[cj)[, 2 tt — (j)i] (see Fig. 8.7). Here is the phase where the upper step of f{E) 
matches the upper Andreev level. Consequently, the critical current, which 
was originally found at <f> = n with a value of IcO = sAqITi, is now reduced 
to a lower value Ic- This effect can be used as a control mechanism, since by 
changing the applied voltage V on the normal conductor, the critical current 
can be adjusted. 




Fig. 8.8. Critical current as a function of the applied voltage for various values 
of the Z parameter at T = 0 {solid lines) and T — 0.2Tc {dashed lines). The inset 
shows the normalized critical current Ic{T, V = 0)/7co vs. T for the corresponding 
Z parameters [210] 



In Fig. 8.8 the critical current, calculated from (8.18), is plotted as a 
function of the applied voltage. In case of no barrier at the interface {Z = 0) 
the critical current decreases immediately at T = 0, as well as at T = 0.2 Tc 
if a voltage V is applied. This behavior can be understood by considering 
the Andreev-level spectrum sketched in Fig. 8.7. Even a small voltage leads 
to a rearrangement of the Andreev-level filling in the vicinity of a phase 
difference tt. The situation changes completely if a barrier is present at the 
interface {Z > 0), since now a gap appears in the Andreev -level spectrum 
(see Fig. 6.5). Here, a finite bias voltage is necessary to cause the upper 
Andreev level to be occupied. Consequently, up to bias voltages eV < Eg^p, 
no influence on the critical current is expected at T = 0; a plateau develops. 
Increasing the temperature leads to a shortening of this plateau owing to the 
thermal smearing of the nonequilibrium distribution function. 



8.2 Supercurrent Control by Nonequilibrium Carriers: Experiments 121 

8.2 Supercurrent Control by Nonequilibrium Carriers: 
Experiments 

Concerning the experimental realization of Josephson junctions which can 
be controlled by injected carriers, two type of structures will be introduced. 
First, an InAs step edge junction is discussed, where the carriers are injected 
from a 5-doped layer into the channel along the step. In the remaining part of 
this section it will be shown how carrier injection can be utilized in junctions 
with a buried two-dimensional electron gas in a heterostructure. 

8.2.1 Carrier Injection in an InAs Step Junction 

An InAs step edge junction, as depicted in Fig. 3.5a, has the advantage that 
the channel length is simply defined by the depth of the step. Furthermore, 
the characteristics of the junction can be modified by using an epitaxially 
grown layer system. In Fig. 8.9 such a type of InAs step edge junction is 
shown, where an inserted p-type Adoped layer is introduced, in order to 
control the two-dimensional channel along the step. The conductive 5-doped 
layer is contacted by an additional An electrode. The underlying idea of this 
kind of structure is that a barrier is formed in the two-dimensional channel 
owing to the pn junction at the intersection between the p-type 5-doped 
layer and the surface channel. This barrier determines the coupling e to the 
junction. By applying an appropriate voltage to the control electrode, the 
barrier height can be controlled by changing the space charge area. Owing to 
the relatively high transparency of the barrier formed by the pn junction, hot- 
carrier injection is a major contribution to the characteristics of the structure. 

In order to obtain isolation between the different terminals of the struc- 
ture, an InGaAs barrier combined with an InGaAs cap layer was grown in 
addition [32]. Only on those areas where this bilayer was removed by wet 
chemical etching was the InAs surface two-dimensional electron gas recov- 
ered. This was done before the two Nb electrodes were deposited, as shown 
in Fig. 8.9. 

A typical set of current-voltage characteristics of an InAs step junction 
with a 5-doped layer as an injector is shown in Fig. 8.10a. For zero voltage 
applied to the control electrode, the junction has a critical current Ic of 
320 pA and a normal resistance i?N of 1.6 Q at 4.2 K, thus yielding an Ici?N 
product of 510 pV. When a voltage is applied to the control electrode, the 
width of the superconducting state in the characteristic is reduced. This can 
be seen better in Fig. 8.10b, where the difference 2/c = Ic+-Ic- between the 
positive value ^c+ and the negative value Ic- of the critical current is plotted. 
For negative control voltages the double critical current 2/c is reduced from 
an initial value of 640 pA to a value below 500 pA at a control voltage of 
—500 mV. A maximal superconducting transconductance 5 ^, as defined in 
(7.4), of 12 mS/mm was obtained for the control voltage interval [—200 mV, 
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Fig. 8.9. Schematic illustration of an luAs step junction. The supercurrent flows 
in the surface two-dimensional electron gas along the step between the two Nb 
electrodes. The junction has a length of 120 nm and a width of 20 |tm. The InAlAs 
and luGaAs layers serve as a barrier and a cap layer, respectively. Control of the 
junction is obtained by injecting carriers via the p-type 5-doped layer. The Si02 
layer is used to prevent leakage currents from the An electrodes 



— 100 mV], which is larger than the values obtained for JoFET structures 
[11]. For positive voltages applied to the control electrode, 2/c first increases 
slightly and then decreases for Vg > 200 mV. 





Fig. 8.10. (a) Current -voltage characteristics of an InAs step edge junction at 
4.2 K controlled by a current injected from a 5-doped layer. From left to right, 
voltages Vg of —300, —200, ..., 200 and 300 mV were applied between the control 
electrode and the upper Nb electrode, (b) The difference between the positive value 
(7c+) and the negative value (7c-) of the critical current (27c = 7c+ — 7c-) as a 
function of Vg. The inset shows the positive value (7c+) and negative value (7c_) of 
27c [20] 
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The observed shift of the current-voltage characteristic (Fig. 8.10) when 
a control voltage is applied indicates that a considerable amount of carriers 
are injected from the (5-doped layer into the channel. The shift is due to the 
fact that this current is superimposed on the current flowing between the two 
superconducting Nb electrodes. The large injection current indicates that the 
barrier formed by the pn junction is relatively transparent. 

The asymmetry shown in Fig. 8.10b can be explained by the pn junction 
at the intersection between the p-type J-doped layer and the n-type surface 
two-dimensional channel. As mentioned above, owing to this pn junction a 
barrier arises between these two degenerate sections [77] . If a negative control 
voltage is applied, the barrier height is increased and the carrier concentration 
is reduced owing to pinning at the surface. The resulting higher resistance of 
the channel leads to the observed reduction of the critical current. If a positive 
voltage is applied, the barrier in the channel is lowered, leading to an increase 
of Ic- Comparison with the theoretical models in the preceding section shows 
that in this structure it is possible to control the coupling strength e with the 
applied voltage. The reduction of Ic for higher positive control voltages can 
be attributed to heating effects. Owing to high transparency of the barrier, 
a considerable current flows from the J-doped layer into the two-dimensional 
channel. At zero bias a low dV/dl value of 140 Q was observed. Thus the large 
current injected from the J-doped layer into the channel leads to a heating 
and thus to a broadening of the Fermi distribution function. As explained in 
the preceding section, this leads to a rearrangement of the occupation of the 
Andreev bound states in the junction. The net supercurrent decreases owing 
to the occupation of Andreev states carrying a supercurrent in the opposite 
direction. 

When the p-type (5-doped layer was replaced by an n-type (5-doped layer, 
a reduction of the critical current was observed as well, but now almost 
symmetrically with respect to the polarity of the control voltage or current 
[32]. In this case no barrier is formed between the two-dimensional channel 
at the step edge and the (5-doped layer. The change in the supercurrent is 
in this case solely due to carrier injection effects, without any change of the 
coupling of the injector to the channel. The resulting heating of the carriers 
by hot-carrier injection results in a reduction of the critical current. 

8.2.2 Multiterminal Josephson Junctions 

In order to feed an additional current through the buried 2DEG located be- 
tween the superconducting electrodes, ohmic contacts have to be fabricated 
in addition. A top view of a typical junction with normal-conducting 2DEG 
leads is shown in Fig. 8. II, while a detail of the junction is depicted in 
Fig. 8.12. For the InGaAs/InP heterostructures, an Ni/Au:Ge/Ni layer sys- 
tem was used. The four ohmic contact areas are connected to the junction via 
two-dimensional-electron-gas leads. As an alternative to structures based on 
an InGaAs/InP layer system, an Alo. 2 Gao. 8 Sb/InAs heterostructure grown 
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by molecular beam epitaxy can also be used as a 2DEG [85]. The ohmic con- 
tacts for the corresponding Nb-AlGaSb/InAs-Nb junctions were fabricated 
by using a Pd/ An layer system [23]. 




Fig. 8.11. Top view of an Nb-IuGaAs/InP-Nb junction with additional ohmic 
contacts to the normal-conductive area of the junction 




Fig. 8.12. Scanning electron micrograph of an Nb-InGaAs/InP juirction. The 
normal-conductive area of the junction is connected on both sides to ohmic con- 
tacts so that a current /b can be injected into the junction. The separation of the 
superconducting electrodes shown is 400 nm, while the width of the Nb electrodes 
is 6 |lm. The arrows indicate the flow of the bias current 7bc and the injection 
current Is 
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The fact that an injection current /b can be utilized to control a Josephson 
junction is demonstrated in Fig. 8.13. Referring to the measurement setup 
shown in Fig. 8.14, the control current is fed from terminal B into the Nb elec- 
trode C, which serves as a common ground for this circuit. For small control 
currents (/b < 50 nA) a strong effect on the junction I-V characteristic is ob- 
served, while for larger values of Ib the effect on the characteristic is smaller. 
A closer inspection of the curves in Fig. 8.13 shows that the I-V character- 
istics are slightly shifted upwards for reverse control currents (/b < 0) and 
downwards for the forward control currents. Therefore the switching of the 
junction into the resistive state in the forward direction ^c+ and in the reverse 
direction 1^- was plotted as two separate lines in Fig. 8.15 for the following 
analysis. 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 



v,c(mv) 

Fig. 8.13. Current -voltage characteristics for different injection currents, varied as 
a parameter. The curves for reverse (Ib < 0) and forward (7b > 0) control currents 
are shifted in steps of 50 |lV [22] 



The critical current 7c = (7c -i- — 7c_)/2 and the offset current FoSset = 
(7c+ -I- 7c-)/2 can be determined from 7c_ and 7c+ (Fig. 8.15). The offset 
current is due to the fact that in the superconducting state of the junction 
both Nb electrodes are at the same electrical potential. Therefore, part of 
the injection current flows first into the opposite electrode E, as depicted in 
Fig. 8.14. Subsequently, this current component is transferred into electrode 
C as a supercurrent. Owing to this current contribution, the threshold of 
Iec for switching into the resistive state is shifted. In fact, in this particular 
junction the current contribution which flows first into E is even larger than 
the contribution which directly flows into C, since the /offset (f^e) curve has 
a slope of —0.7, instead of the value of —0.5 for symmetric injection. The 
asymmetric injection is due to different interface transparencies on the two 
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Fig. 8.14. Schematic illustration of the setup for a three-terminal measurement. 
Contact C is used as a common ground. The arrows show the current paths in the 
junction. The broken arrow represents the injection current which flows first iirto 
the opposite electrode E 




-0.4 -0.2 0.0 0.2 0.4 

IbCI^A) 

Fig. 8.15. Upper critical current 7c+ (A) and lower critical current 7c_ (v) as 
a function of the injection current 7b- The critical current 7c (•) and the offset 
current 7ofjset (*) were determined from 7c± [22] 



sides of the junction. For this sample, the resistance from B to E was higher 
by 20 Q than the resistance from B to C. 

As noted above, for small injection currents (Jb A 50 nA) a very strong 
effect on the critical current 7c is observed, e.g. 7 b = 50 nA is sufficient 
to reduce the initial critical current of 1.48 |lA by 43%. As can be seen in 
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Fig. 8.16. |Z\7c/Z\/b| as a function of the injection current. The inset shows Ic (v) 
and Jr (A) vs. temperature [22] 



Fig. 8.16, a ratio of critical current to control current \AIc/AIb \ as large as 20 
was observed at low values of |/b|- For larger injection currents the decrease 
of Ic is slower. 

The model describing the effect of nonequilibrium carriers on the Joseph- 
son effect introduced in Sect. 8.1 can be utilized to explain the measured 
data. The condition that the junction should be in the short limit is fulfilled, 
since only two Andreev bound states are present below /i (see Fig. 6.6). Nev- 
ertheless, some differences have to be considered. First, in the experiment a 
current Jb was injected into the junction, while in the theoretical model a 
voltage was applied to the normal leads. Second, the junction discussed in 
the model was one-dimensional, whereas the measured junction has a width 
of 6 pm. Concerning the first statement, the injection of a current into the 
normal-conducting region of the junction is of course accompanied a volt- 
age drop between the semiconductor and the superconductor. However, the 
voltage drop is not constant but, rather, decreases from the entrance point 
of the junction to the opposite side. Typical values for the voltage drop at 
the Nb/2DEG interface were 0.18 mV at the injection point and 5 pV at the 
other end for an injection current of Jb = 0.1 pA. Consequently, in a more 
elaborate model an averaging along the junction would be desirable. 

Comparing the theoretical results plotted in Fig. 8.8 with the experimen- 
tal Ic vs. Jb values shows that the plateau predicted in the model for a finite 
Z parameter is not observed in the experiment. Nevertheless, since the tem- 
perature dependence of the critical current tends to saturate at low tempera- 
tures in accordance with the finding discussed in Sec. 6.2.1, a barrier must be 
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present at the interface and thus a finite Z parameter can be assumed. One 
possible reason for the absence of the plateau in the Ic{Ib) measurements is 
the two-dimensional character of the junction, which implies that the gap in 
the Andreev-level spectrum is lowered owing to wave vectors with a compo- 
nent parallel to the Nb/2DEG interface. In addition, Ilhan and Bagwell [209] 
showed that even if the supercurrent switches abruptly as a function of the 
injection voltage, the supercurrent changes only slowly as a function of the 
injection current. 

As demonstrated by Morpurgo et al. [21], an injection current can lead to 
a heating of the normal-conductive region of the junction up to an effective 
temperature T* [211]. In this model, the broadening of the distribution func- 
tion f{E) due to electron-electron scattering processes leads to a decrease of 
the net supercurrent. However, as explained below, in contrast to the results 
obtained on Nb/Au junctions [21], it is not possible to explain the sharp 
decrease of the critical current in the Nb/2DEG structures solely by heating. 

For the following discussion we shall refer to experiments on Nb-AlGaSb/ 
InAs-Nb junctions [23, 107]. During these measurements the injection cur- 
rent Ib.b' was driven through the 2DEG from the ohmic contact B on one 
side of the junction to the contact B' on the opposite side (Fig. 8.14). In 
accordance with the results in Fig. 8.15, the critical current first decreases 
strongly for smaller injection currents, while a slower decrease is observed for 
larger injections currents (Fig. 8.17). From the voltage drop Vjnj measured 
between the two ends of the control line, the effective temperature T* was 
determined by using the following expression [211]: 

T*{y) = + (j//Aw)(l - 2//Lw)(3eV7r2fc|)Vi2. . (8.19) 

For the analysis, the effective temperature T*{y = L^j2) at the center of the 
junction, with T = 0.6 K, was taken as an average value. In Fig. 8.17 (inset) 
the critical current is plotted as a function of T and T*. Obviously, the Ic{T*) 
curve beginning at T = 0.6 K decreases considerably more steeply than the 
/c(T) curve; only at higher temperatures T, T* do both curves merge. 

From the difference between the two curves at low temperatures it can 
be concluded that, at least at low injection currents, different mechanisms 
are responsible for the steep decrease of Ic with Ib and with T* . In order to 
obtain a broadening of the distribution function, which can be interpreted 
as a higher effective temperature T*, electron-electron scattering processes 
are necessary. By using the formula of Giuliani and Quinn [213], the inelastic 
scattering length kn,ee can be estimated for typical injection energies. For an 
excess energy of 2 meV and a Fermi energy of 30 meV, ^in,ee ~ 7pm. This re- 
sult shows that even for injection energies where the critical current is almost 
completely suppressed, the inelastic scattering length is larger than the width 
of the junction. This result leads us to the conclusion that the distribution re- 
sembles more a step-like distribution function, as in the model introduced in 
Sect. 8.1.2, than a broadened distribution function. This statement holds for 
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Fig. 8.17. Critical current vs. injection current fed from terminal B to ter- 
minal B' for an AlGaSb/InAs-based junction. The inset shows 7c as a function of 
the bath temperature T and effective temperature T* [210] 



the AlGaSb/InAs-based junction as well as for the Nb-InGaAs/InP-Nb junc- 
tion. However, for larger injection currents heating effects due to electron- 
electron scattering will at least partially contribute to the reduction of 7c, 
which can be concluded from the fact that both curves in Fig. 8.17 (inset) 
merge for larger temperatures T, T*. 

Morpurgo et al. [214, 215, 21G] used a junction connected to a super- 
conducting ring, similar to the configuration shown in Fig. 8.1, in order to 
detect the discrete Andreev levels in the junction formed by the two super- 
conducting terminals of the ring and the 2DEG in between. For an increasing 
magnetic field, the differential resistance dV/dl of this structure showed an 
oscillation pattern with a period of the flux quantum threaded through the 
ring. If now an injection voltage was applied, a change of the dV/dI{B) os- 
cillation amplitude, accompanied by a phase shift, was observed. By using a 
microscopic model, based on an analysis of the interference processes occur- 
ring within the junction, Morpurgo et al. could explain the modification of 
the oscillation pattern by the presence of discrete Andreev bound states. 

In contrast to junctions with a semiconducting channel, a reversal of the 
supercurrent, as discussed in Sect. 8.1, was observed by Baselmans et al. 
[217, 218] in a purely metallic four-terminal Nb/Au junction. The observed 
effects could be explained consistently by the model in the diffusive regime 
of Wilhelm et al. [205]. 

By employing superconducting aluminum electrodes as injectors, Kutchin- 
sky et al. [219] obtained an enhancement of the critical current in an Al/n- 
type GaAs/Al Josephson junction when the injection voltage reached Aq/c. 
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This effect can be explained by the peak of the density of states in the super- 
conducting injectors at Aq. If this singularity in the density of states of the 
injectors matches the chemical potential of the junction, the critical current 
is enhanced. 

8.2.3 Local Suppression of the Josephson Current 

Depending on the width of the structure, carrier injection from one side or 
the other into the junction is not always sufficient to suppress the Josephson 
supercurrent completely. This can be seen in Fig. 8.18, where the critical 
current of the junction is plotted for different modes of carrier injection. In 
the case of a carrier injection from either contact A or contact B into SI, a 
finite critical current is found even for large injection currents. In contrast, 
if the injection current is fed from contact A to B a complete suppession is 
possible. 




Injection current (pA) 

Fig. 8.18. Critical current as a function of the injection current for three differ- 
ent conhgurations for carrier injection. Contacts A and B are normal-conducting 
terminals. The critical current was determined for the junction formed by the su- 
perconducting electrodes SI and S2 with the two-dimensional channel in between. 
The results of injection from either A or B into SI, as well as for a control current 
passing from A directly into B, are plotted 



The partial suppression can be explained by the fact that for injection 
from A or B into SI only part of the normal-conducting area between the 
superconducting electrodes is affected. This is depicted in the schematic illus- 
tration attached to the curves in Fig. 8.18. As can be seen here, the injection 
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current affects only the junction area in the vicinity of the injection point, 
whereas the area on the opposite side is more or less unaffected. Owing to the 
S/2DEG interface, which is not perfectly homogeneous, the efficiency of the 
supercurrent suppression differs according to whether the current is injected 
from A or B. In the case of injection from A through the normal region of 
the junction into B a complete suppression is observed, which shows that in 
this case the whole 2DEG area of the junction is affected. 

More detailed information about how the injection current influences the 
supercurrent distribution can be gained with the help of a magnetic field. In 
Sect. 6.2.2 it was demonstrated that by analyzing the interference pattern of 
the critical current as a function of an external magnetic field, information 
about the supercurrent distribution in the junction area can be extracted. 
This method can now be employed in order to trace the region where the 
injection current affects the supercurrent flow. 

In Fig. 8.19 the Ic{'P/^o) interference pattern of an AlGaSb/InAs-based 
junction, which will serve as a reference, is plotted [23]. Here, no current is 







Fig. 8.19. /c(^/^o) pattern of an AlGaSb/InAs-based junction without injection 
current {solid line) and a fit using the distribution in the inset {dash-dotted line) 
[..;.] (© (1999) by the American Physical Society) 



injected from a normal lead. By using a generalized form of (6.34), the current 
distribution in the junction can be simulated. As can be seen in Fig. 8.19, the 
supercurrent distribution is largest at the left side, while almost no current is 
flowing in the middle. The higher current densities at the edges of the contacts 
can be explained by the fact that this junction can be regarded as being in the 
intermediate range between a short and a long Josephson junction {W ~ Aj, 
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where Aj is the Josephson penetration depth). The asymmetry is presumably 
due to inhomogeneities at the superconductor/2DEG interface. 

The effect of an injection current on the supercurrent distribution in the 
junction is visualized in Fig. 8.20. The upper graph shows the Ic vs. 




Fig. 8.20. (a) /c(^/^o) for an injection current of 500 nA from the left-hand 
side (solid line) and a simulation (dash-dotted line) for the modified supercurrent 
distribution shown in the inset. The side of the current injection is marked by 
an arrow, (b) Corresponding measurement and simulation for injection from the 
opposite side [23] (© (1999) by the American Physical Society) 



characteristic for the case where a current of 500 nA is injected from the left. 
Here, a three-terminal configuration as depicted in Fig. 8.14 was used. The 
experimental interference pattern was fitted as described above by assum- 
ing a step-like supercurrent distribution. Now, the best fit is obtained for a 
distribution where the supercurrent density is reduced on the left, as shown 
in Fig. 8.20 (inset). The interference pattern calculated using this modified 
current distribution matches well to the experimental curve. Since now the 
current distribution is more balanced, the interference pattern resembles the 
shape corresponding to a superconducting quantum interference device. 

Analogously to the case of injection from the left-hand side, an injec- 
tion current from the right-hand side leads to a decrease of the supercurrent 
density on this side, as shown in Fig. 8.20. Owing to the lowering of the 
supercurrent in the area where the current is injected, the dominant contri- 
bution to the supercurrent is now found on the opposite side. As a result, the 
interference pattern has approximately the shape corresponding to a shorter 
junction. 



9. Conclusion 



In this book the properties of superconductor/III-V semiconductor hybrid 
structures with a two-dimensional electron gas as a conducting channel have 
been reviewed. A crucial parameter for all of these structures is the trans- 
parency of the superconductor/ two-dimensional-electron-gas interface. Opti- 
mizing the transparency is of course mainly a technological challenge but, 
as was shown, a detailed understanding of the transport processes is neces- 
sary in order to assess the experimental results obtained. At present, argon 
sputtering is the best-established interface-cleaning technique and is used in 
many laboratories. However, there is still room for further improvement. In 
this respect He plasma cleaning seems to be a promising method, since here 
the semiconductor crystal is not degraded by the sputtering process. 

Comparing Josephson junctions that contain an inversion layer with junc- 
tions based on a buried two-dimensional electron gas in a heterojunction 
shows that the former type of structure has some advantages concerning a 
large characteristic voltage, which is an important parameter with respect 
to possible applications. However, if more fundamental questions about the 
transport mechanisms are to be addressed, the latter structures based on het- 
erojunctions are almost ideal model systems, since here the transport can be 
regarded as ballistic. In this regime the transport processes can be described 
by relatively simple models. 

Control of the critical current by light exposure is a typical example of 
how specific properties of semiconductors, in this case persistent photocon- 
ductance, can be employed. Another rather unique property is the large Fermi 
wavelength of the semiconductor, which is of the order of the feature size ac- 
cessible by electron beam lithography. Owing to this large Fermi wavelength, 
remarkable quantum effects, such as the quantization of the critical current, 
can be observed very clearly. 

In principle, a Josephson element with a semiconductor bridge between 
two superconducting electrodes can serve as a basis for a three-terminal elec- 
tronic device, since the electron concentration of the semiconductor can be 
controlled easily by a gate. However, the performance of the Josephson field 
effect transistors obtained so far shows that, owing to many technological 
problems, it is questionable whether this concept will find its way into real 
electronic circuits. 
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9. Conclusion 



Control of the Josephson current can also be obtained by injecting 
nonequilibrium carriers via an additional ohmic contact into the normal- 
conductive area of the junction. This effect can be explained theoretically 
within the framework of a nonequilibrium occupation of phase-coherent An- 
dreev bound states. Currently these kinds of multiterminal structures are 
being intensively studied, theoretically as well as experimentally. Here, one 
aspect is their potential application as a transistor-like device. 

Regarding future research on superconductor/two-dimensional-electron- 
gas structures, it will be very interesting to study the Andreev reflection 
process in large magnetic fields. Owing to the Landau quantization so-called 
edge states are formed in the two-dimensional electron gas, which couple 
to the superconductor. In order to preserve the superconductive state in a 
magnetic field, different electrode materials, e.g. NbN, have to be employed. 
Initial steps in this direction have been taken already [220, 221]. Another field 
of interest is the coupling of superconducting electrodes to a semiconductor 
quantum dot. Concerning theory, many effects have been predicted already 
[222, 223, 224, 225], but experiments on hybrid quantum dot structures are 
still lacking. 
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99 

superconductive ground state, 6 
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supercurrent, 65 

supercurrent distribution, 83, 131 
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surface states, 22 
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